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PREFACE 


This book has been written at the suggestion of some students 
who have been unable to find a book on the subject suitable 
for their requirements, and also with a view to completing 
the series of text-books on Mechanics which I have produced 
in recent years. 

The book is only an introduction to a subject which, in the 
eighteenth and nineteenth centuries, occupied the interest 
and attention of many eminent mathematicians. Nowadays 
the theory of potential is usually studied in connection with 
electricity and magnetism, which is given priority as a subject 
of study because of its wider developments and importance 
in relation to everyday affairs, so that there is a tendency to 
overlook the theory of the attraction of gravitating solids or 
regard it as a subject to be ‘crowded out’. But the subject 
is of importance in connection with dynamical astronomy 
and the Figure of the Earth, and Todhunter’s history of the 
subject ‘from Newton to Laplace’ ran to nearly a thousand 
pages. 

In the early chapters I have not forgotten the needs of 
students reading for a pass degree, though the book as a whole 
is intended for an honours course. A good many of the 
theorems in pure mathematics specially required for the 
subject are collected in the first chapter, but students who 
prefer to do so may begin to read at Chapter n and only refer 
back to Chapter i as the need arises. The later chapters are 
on Green’s theorem. Harmonic functions and the Attraction 
of ellipsoids. 

I have indicated the source of the examples when taken 
from examination papers, the abbreviations M. T., C,, and P, 
denoting the Mathematical Tripos, a College or Intercollegiate 
Examination, and the Preliminary Examination in Mathe¬ 
matics in the University of Cambridge, and I gratefully 
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acknowledge the kind permission of the Senate of the Uni¬ 
versity of London to use questions from the examinations of 
that University, the source of which is also indicated. 

A line in each set of examples indicates a rough division 
of those which are easier from those which are more difficult 
of solution. 

The proofs have been read by my friend Dr S. Verblunsky 
and again I tender to him my thanks. The book has been 
greatly improved by his valuable criticism and useful sugges¬ 
tions. I also wish to thank the staff of the University Press 
for their careful composition and assistance in eliminating 
errors. 

A. S. R. 

Cambridge 
May 1940 
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Chapter I 

PRELIMINARY MATHEMATICS 

1-1. We propose in this chapter to give a brief account of 
some mathematical ideas and propositions which the reader 
will find of use in later chapters. 

1*2. Surface and volume integrals. Inthetheory of attrac¬ 
tions it is necessary to make considerable use of surface and 
volume integrals, and the applications of the theory frequently 
involve the evaluation of such integrals. The reader is therefore 
advised to acquire some familiarity with methods of integra¬ 
tion from a textbook on the subject. We propose here merely to 
explain what is implied when such symbols as 

jfix,y,z)dS and jf{x,y,z)dv 

are used to denote integration over a surface or through 
a volume, and then to prove some important propositions 
connecting surface and volume integrals. 

A definite integral of a function of one variable, say J / (x) dx, 

may be defined thus: let the interval from a to 6 on the aj-axis 
be divided into any number of sub-intervals Sj, Sg,... and 
let/,, denote the value of f(x) at some point on 8 ,.. Then, pro- 

n 

vided that the limit as oo of S exists and is indepen- 

r»l 

dent of the method of division into sub-intervals and of the 
choice of the point on 8,. at which the value of f{x) is taken, 
this limit is the definite integral of f{x) from a to 6, It can be 
shewn that the conditions for the existence of the limit are 
satisfied iff(x) is a continuous function. 

In the same way we may define y, z) dS over a given 

surface; let the given surface be divided into any number of 


RNA 
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small parts 81 , 82 ,... 8 „ ti-nd let/,, denote the value off {x, y, z) 

n 

at some point on 8 ^, then the limit as -> oo of S 8 ^, provided 

rtml 

the limit exists under the same conditions as aforesaid, is 


defined to be the integral J given surface. 

Any difficulty as to the precise meaning to be attributed 
to 'area of a curved surface’ may be avoided thus: after 
choosing the point on each sub-division 8 ^ of the surface at 
which the value of f{x,y,z) is taken, project this element of 
surface on to the tangent plane at the chosen point, and take 
the plane projection of the element as the measure of 8 ^ in 
forming the sum. 

The integral jf(x,y,z)dv through a given volume may be 
defined in a similar way. 

The evaluation of a surface integral of course involves 
a double integration, and that of a volume integral a triple 
integration, and when it is desired to exhibit the forms chosen 
for the elements of area or elements of volume such symbols 


asjj/(a;, y, z) dxdy JJj" f dxdydz are used; but when 


there is no doubt as to what is implied a single integration 
symbol f is employed as above for the representation of an 


integral along a fine or over a surface or throughout a volume. 


1*3. Green's Theorem- Ifn, v^w are functions of x^ y, z which 
have continuous derivatives with respect to x, y, z respectively 
throughout a singly-connected region bounded by a closed surface 
8 free from singularities and I, m, n denote the direction cosines 
of the normal drawn outwards from a point on an element dS 
of the surface, then 

///(£ ^ 5 ^ S +mv + nw) dS, 

where the surface integral is taken over the boundary and the 
volume integral throughout the region enclosed. 
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Consider the integral 

section dxdy parallel to the 
axis of z. Let it intersect the 
surface 8 in elements dS-^, 
dS^ and let normals to these 
elements have direction co¬ 
sines Zi, Zgj ^ 2 * 

Since the elements dS-^, dS^ 
project into an area dxdy on 
the ajy-plane, taking account 
of signs we have, from the 
figure, 



dxdydz. Take a prism of cross 



dxdy= —n^dSy^ — n^dS^ .(1). 


Hence, by integration with regard to z, we have for this 
particular prism 

J ^ dxdydz = {wdxdy).^ — {wdxdy )^, 

the limiting values corresponding to the intersections of the 
prism with the surface, so that 

The integration through the whole bounded region is now 
completed by summing for all similar parallel prisms which 
intersect the surface 8, so that 


=jjnwd8. 


The integrals relating to the functions u and v may be 
treated similarly by considering prisms parallel to the x and y 
axes, so that 
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This is Green’s Theorem* connecting a volume integral 
with a surface integral over the boundary of the volume. 

1-31, Green's Second Theorem. In the last theorem make 


dV' dV' 

the substitutions u = V -^r-, ^ = F 

dx dy 


where F, F' 
oz 


are functions which with their first and second derivatives 
are finite and continuous through the region considered. 
Substituting in 1*3 (2), we obtain 




aFaF^ aF aF^ dVdV' 

8x dx dy dy ^ dz dz 


j dxdydz 


-m 


dW dW' d<‘V' 
02® 




- v r w r w r 

dx dy dz 


j dxdydz 
.( 1 ). 


But if djdrii denotes differentiation in the direction of the 

/1/Tr/v*^ ^ .O Vi o ttq 7 /vin nn ~ ^ 


ouiward-drawn normal to dS, we have l,m,n = 


so that 


duy^ dfii* dn^* 


dx dy dz dx drii dy drii dz dn^ 


Hence (1) is equivalent to 


rrr/dvdv' dvdv^ dVdV'\^ , , 

JJJ\Sa: dx'^ dy dy dz 02 j * ^ ^ 

“ J J J J J ^ 

and since V and V' are clearly interchangeable, this also 




dW dW 0®F\ 

-(2). 


* George Green (1793—1841), published his Essay on the AppU<xiHon 
of M.athevnaHcal Analysis to the Theory of Electricity and JMagn&tism 
in 1828. At the age of 40 he came to Cambridge and was foxirth 
wraiigler in 1837. 
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where the volume integrals are taken throughout the region 
considered and the surface integrals over its boundary. 

1*4. Solid angles. The solid angle of a cone is measured 
by the area intercepted by the cone on the surface of a sphere of 
unit radius having its centre at the vertex of the cone. 

The solid angle subtended at a point by a surface of any form 
is measured by the solid angle of the cone whose vertex is at 
the given point and whose base is the given surface. 



Let PP' be a small element of area dS which subtends a 
sohd angle do) at O. 

Let the normal to dS make an acute angle y with OP, and 

let OP = r. Then the cross-section at P of the cone which dS 

subtends at 0 is dS cos y, and this cross-section and the small 

area dw intercepted on the unit sphere are similar figures, 

so that , cr j 2 1 

aocosy: dm — r^ : 1. 


Whence 

or 


dw^{d8 ooBy)/r^ 
d8==r^BeGy dw 


{!)• 


It follows that the area of a finite surface can be represented 
as an integral over a spherical surface, thus 


8 



with suitable limits of integration. 


(2), 
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1*41. rfco in polar co-ordinates, du) is an element of the 
surface of a unit sphere. Let the 
element be PQRS bounded by 
meridians and small circles, where 
the angular co-ordinates of P are 
d, <f>. Then since the arc PS sub¬ 
tends an angle at the centre of 
a circle of radius sin^, therefore 
PS = sin 6dcl >; and PQ = dd, so that 

dcu = PQ . PS = sin ddddtf). 

1*42. Volumes in terms of solid angles. Consider a cone 
of small sohd angle da> and 
length B, Let PP\ QQ' be 
cross-sections at distances r, 
r + dr from the vertex 0. 

Then the area PP' = r^dco, 
and, neglecting an infini¬ 
tesimal of higher order, the 
volume of the frustum PQQ'P' 
is r^doidr, 

rR 

Hence the volume of the cone = dca r^dr 

Jo 

= ^B^do). 

Hence if 0 be an origin inside a region bounded by a con¬ 
tinuous surface such that lines from 0 each meet the surface 
in one point only and we divide up the region into narrow 
cones with vertices at 0 and B denotes the radius vector or 
length of the cone of solid angle dw, then 

will represent the volume of the region. 

1*5. Scalar functions of position and their gradients. 

Let F (Xf z) be a single-valued continuous function of the 
position of a point in some region of space. Suppose that the 
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function V is not constant throughout any region, so that 
the equation 

V {x, y, z) = const. 

represents a surface. We assume that through each point of 
the region in which V is 
defined, there passes a sur¬ 
face F = const. We also as¬ 
sume that at every point P 
on this surface there is a 
definite normal PN and that 
the tangent plane at P varies continuously with the position 
of P on the surface. 

It follows from the definition of V that two surfaces 



V {x,y,z)=^a and V {x,y,z) — b 

cannot intersect; for if they had a common point it would be 
a point at which V had more than one value, in contradiction 
to the hypothesis that F is a single-valued function. 

Consider two neighbouring surfaces 

F = a and F = a-hda. 

Let P, P' be points on each and let the normal at P to 
F = a meet V — a + da in N, For small values of da, PN will 
also be a normal to F = a -f da. 

Then using Vp to denote the value of F at P, we have 

Fp.-Fp da F;,~Fp F^~Fp PN 

ppr pp, pp, ppf pp, 


. Vs-Vp 
PN 


COS d. 


where 6 is the angle NPP'. 

Now if PP' = ds and PN = dn, and we make da and therefore 
also ds and dn tend to zero, the limit of —Ip)/PP' is 
the rate of increase of F in the direction d^ and is denoted by 
9F 

and similarly the limit of (^ — Pp)/P^ is the rate of 
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increase of F in the normal direction dn and is denoted by 

and we have 
dn 


dv dv ^ 

_=_COS0 

ds dn 


.( 1 ). 


Thus we have proved that the space rate of increase of V 
in any direction ds is the component in that direction of its 
space rate of increase in the direction normal to the surface 
V ~ const.; or that if we construct a vector of magnitude dVjdn 
in direction PN, then the component of this vector in any 
direction is the space rate of increase of V in that direction. 

This vector is called the gradient of V and written grad F. 

To recapitulate: F is a continuous scalar function of position 
having a definite single value at each point of a certain region 
of space, and the gradient of F is defined in this way: through 
any point P in the region there passes a surface F = const., 
then a vector, normal to this surface at P, whose magnitude 
is the space rate of increase of F in this normal direction, is 
defined to be the gradient of V at P, and it has the property 
that its component in any direction gives the space rate of 
increase of F in that direction. It is clear that the gradient 
measures the greatest rate of increase of F at a point. 

1*6. Vectorial methods. We shall denote vectors by letters 
in Clarendon type. The magnitude of P is denoted by | P |. 
We recall that vectors are compounded by the parallelogram 
law and that they obey the commutative and associative laws 
of addition, viz. that 

P + Q = Q + P, and that P + (Q-f R) = (P + Q) + R. 

The product of a vector P and a scalar m is a vector whose 
magnitude is the product of the magnitude of P and m, i.e. 
m I P |, whose direction is that of P or the opposite according 
as m is positive or negative. Multiplication by scalars follows 
the ordinary rules of algebra, viz. 

m (n’P)=^n(mP) = {mn) P; 
m(P+ Q) = mP-f mQ 
and (m + n)P = mP + nP. 




VECTORIAL METHODS 


9 


1 - 62 ] 

A vector s is called a unit vector when | s | = 1. Any vector 
whose direction is that of a unit vector s can be represented 
by ms, where m is a scalar. 

For convenience in representing vectors by rectangular 
components we introduce three fundamental unit vectors 
i, i, k whose directions are mutually perpendicular; these 
directions being those of the co-ordinate axes Ox, Oy, Oz, 
Then, Hl Py, be the projections of the vector P on the 
axes, we have a vector equation 

F=Pj+Pyi + P,}L .( 1 ). 

1*61. Scalar product. Let P, Q be two vectors and 
the angle between their positive directions, then the expression 

|P||Q|coS(^ 

is called the scalar product of the vectors and is denoted by 
(PQ) or (P, Q), or simply PQ. 

The scalar product of two vectors at right angles is zero, 
so that for the fundamental vectors i, j, k we have 
ij = jk = ki = 0 and ii = jj = kk=l. 

Hence, from 1*6 (1), 

FQ^{PJ + Pyi + P,k) (QJ + Qyi^ g,k); 

and on multiplying out this gives 

lPQ = Pa.Q^+PyQy^P^Q^ .( 1 ). 


1*62, Again if V denotes a function of x, y, z, it follows from 
1*5 that the projections on the axes of the gradient of V are 
dVjdx, dVjdy, dVjdz. Hence we have the relation 


.37 .37 ,37 


and the equivalence of the operators grad and i^^^ + j^+kx-, 

ox oy oz 

when the latter is applied to a scalar function. This latter 
operator is known as the HamUtoniau operator* and denoted 
by V. When applied to a scalar function it turns it into a vector. 


♦ After Sir William Rowan Hamilton (1806-1865), Astronomer 
Royal for Ireland, greatest of Irish Mathematicians, inventor of 
Quaternions* 
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We shall see in 1*65 that when V is applied to a vector it turns 
it into a scalar. 

The reader who wishes for further information on the 
subject of vector algebra may find it in the author’s Dynamics, 
Part II, Appendix. 

1*63. A vector field. If to every point of a given region 
there corresponds a definite vector A, which may vary in 
magnitude and direction from point to point, then the region 
is called a vector field, or the field of the vector A; e.g. 
gravitational field, electric field. 


1*64, Flux of a vector. If a surface 8 be drawn in the 
field of a vector A and denotes the component of A normal 

to an element d8 of the surface, then the integral ^A^dS is 

called the flux of A through 8 . Since a surface has two sides 
the sense of the normal must be taken into account, and the 
sign of the flux is changed when the sense of the normal is 
changed. The flux of a vector through a surface is clearly 
a scalar magnitude, in fact if we use the symbol dS to represent 
a vector of magnitude d8 directed along the normal to the 
surface, then the scalar product of the vectors A and dS 
denoted by (AdS) is the same thing as the subject of integra¬ 
tion A^dS, so that the flux of the vector may also be repre¬ 
sented by f(AdS). 


1*6S. Divergence of a vector field. Let A denote a 
vector field which has no discontinuities throughout a given 
region. Let 8v denote any small element of volume containing 

a point P in the region and let ^A^d8 denote the outward 

flux of A through the boundary of hv, then the limit as Sr 0 


of ^ jhv is defined to be the divergence of A at the point 


P and denoted by div A. 

Let P be the point {x,y,z). Let A^, Ay, denote rectan 
gular components of the vector A. 
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Let I, m, n be direction cosines of the normal to d8. 


Then 
so that 


A„ = lA^+mAy + nAg, 
jA^dS =J [lA^ + mAy + nA^ dS 


=J( 


dA, dA„ dAA , 


by Green’s theorem (1*3), where the volume integral is taken 
through 8v, The value of the integral is the product of 8v and 
the mean value of the integrand, say the value at a point Q- 
As 8v->0, this mean value tends to the value at P. 


Hence 


(A^dS 
divA= lim —r- 

dt? 


dA^ dA^j dA^ 

dx dy dz 


(1), 


and we note that this result is independent of the form of the 
element of volume 8v which surrounds the point P, and also 
independent of the way in which the element contracts to 
zero. 

Further, in the notation of 1*61, 1-62, we have 

^dA^ dAy SA^ 
dx dy dz 

= div A, 

so that when the operation V is performed on a vector A the 
result is the scalar function div A. 

Also if A = gradF (1*5), where F is a scalar function of 


X, y, z, so that 

A A A ^ ^ 

dV 

dV 


dy' 

dz’ 

we have div grad F = div A 



dA^ 

dA.. 

dA. 
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02 02 02 

The operator + ^ + ^ usually denoted by V®, for 
it is the same as 


a . 9 , 


5* 


l)( 


.a .a , a\ 


The reader will note that the operator V may be applied either to 
a scalar, turning it into a vector, or to a vector, turning it into a scalar; 
i.e. that V means ’ when applied to a scalar and it means ‘ div ’ 

when applied to a vector. But the operator ^grad ’ only has a meaning 
when applied to a scalar function, and ‘ div ’ only has a meaning when 
applied to a vector. Also the order of the operations is not inter¬ 
changeable, thus 


d^V dW dW 

div grad = + +-^,- = VW, 


but grad div V is meaningless. 

Again, div grad A is meaningless, but 


grad divA = 


.a .a , a 

■=-+>5,+'‘ 


’ dx 


i)( 


e^a , dA^ 
dx dy 


dz /* 


1*66. If we assxime that the divergence as defined in 1 *65 is indepen¬ 
dent of the shape of St?, it is easy to obtain the result (1) by taking for 
dv a rectangular parallelepiped with edges parallel to the axes and of 
lengths 5a;, dy^ dz with P at its centre. The method is given in full in 
the author’s Electricity and Magnetism, pp. 6, 7. We shall not reproduce 
it here, as in the next article we shall apply the same method to a 
more general case. 


1*67. Curvilinear co-ordinates. Gradient and Diver¬ 
gence. Let the equations 

f 2 (x,y,z) = ^, fs(x,y,z)==y .(1), 

in which a, y denote variable parameters, represent three 
families of orthogonal surfaces; i.e. surfaces which cut at 
right angles wherever they intersect one another. 

We assume that one and only one surface of each family 
passes through any point {x,y,z), so that the parameters 
a, jS, y may be regarded as curvilinear co-ordinates of the 
point whose cartesian co-ordinates are x, y, z. 

Let denote elements of arc of the curves of 

intersection of the pairs of surfaces jS, y; y, a and a, j8. Then 
since a is the only parameter that changes along the curve 
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jS, y, it follows that ds-^ must be a multiple of da. We may 
therefore suppose that 

ds^^^h^d^f ds^^^h^d^y 


where are in general functions of the co-ordinates.* 

Since the surfaces a, j8, y through any point are at right 
angles, a vector A at the point may be resolved into three 
components normal to these surfaces and denoted by 

-4y. 

Consider an element of volume dv bounded by surfaces 
a ± |da, j8 ± ^dj8, y ± Jdy, Let the surface a cut this element of 
volume in the curvilinear rectangle 
PQRS whose centre 0 is the point 
(*> y)- 

To find an expression for div A in 
curvilinear co-ordinates we have to 
find the flux of A out of dv. The 
component of A at 0 normal to 
PQRS is A^. Hence, if we assume 
what can easily be provedf, that the 
average value of a function over a 
small plane area is the value at the centroid, then the flux 
across PQRS is (area PQRS)^A^.h^d^,h^dy. 



* If the values of h^y are required, they may be found by 
differentiating each of equations (1) with regard to a, ^ and y, and 
solving the resulting equations for dxldtx., dxjdp, etc. Then 


so that 




\da) ’ 


t If/ {Xy y) be the function, S the area and (x, y) its centroid, we have 
jf(x,y)dS=jf{x+x'.y + y')dS 

= j |/ (®* ^+3^' ^ + higher powers of x', j/'j dS 


sinoe x'dS=Jy'dS= 0, 

where i is of order of the greatest linear dimension of S. 




14 


PRELIMINARY MATHEMATICS 


[1-67-- 

Now the parallel surfaces ABCD, EFOH are the surfaces 
a — |da and a + Jda, and the outward flux of A across ABCD is 

— ^Ji2\dpdy — ^^{AJi2h^d^dy) da-f , 

where is an infinitesimal which 0 as rfa 0 . 

Similarly the outward flux across the face EFGH is 

1 9 

AJi^h^d^dy + ^ {Afji^h^d^dy) den 4 - e^doedpdy, 

2 doc 

where €3 has a like meaning to 4 Ej, . 

Therefore the contribution of this pair of faces to the flux 
of A out of dv is 

(* 2 ^ 3 ^a) + dend^dy, 

where e 0 as a -> 0 . 

In like manner for the other pairs of faces of dv^ so that 
the total outward flux of A is 


a 

0 a' 


0 0 \ 

^ (^2^3-4^) + ^ (^3^1^^) + g- (^1^2-^)/) + dend^dy, 

where 0 as dt; -> 0. 


flux A 


Hence, since divA= Um —, and dv^h^h^hodxdpdy 

dv-^O 


therefore 


1 f 0 0 0 

.(3). 

Further, when A is the gradient of a scalar function V, 
we have 

. . . 0F 97 0F 10F 1 9F 197 

■“«> -“( 8 ) Ay = -, 


dsi’ 9 « 2 ’ Ssg kj^dx’ hgdp’ hgdy 


•(4), 


and by substituting (4) in (3), we have divgrad F or 
TOT/_ ^ f ^ 9F\ 9 / fe8fei9F \ 9 


.(6). 
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1*68. Polar and cylindrical co-ordinates. When using 
polar co-ordinates r, 6, the families of orthogonal surfaces 
are concentric spheres r = a, cones 0 == j 8 , and planes ^ = y, and 
it is easy to see that 

ds 2 , ds^^dr, rdd, r&in6d(f>, 


so that hi, ^ 2 , A 3 S 1 , r, rsin^, 

and 



When using cylindrical co-ordinates r, 9, z, the families 
of orthogonal surfaces are co-axial cylinders r = a, planes 0 = 
and planes 2 ; = y, and in this case we have 

dsi, ds^, dsQ = dr, rdd, dz, 
so that hi, ^ 2 , ^ 3 = 1 , r, 1 , 


and 



or 


1^/ m Id^ dW 
rdr\ dz^ 


( 2 ). 


Either of the results ( 1 ) and ( 2 ) can of course be obtained 
directly by applying the method of 1*67 to the appropriately 
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shaped element of volume, and we suggest as an exercise to 
the student the independent verification of results (1) and (2) 
in this way.* 

1*7. If a vector field A is continuous throughout a region R 
bounded by a surface S, then the outward flux of A through the 
surface S is equal to the volume integral of div A throughout JS. 

For we may subdivide R into a large number n of elements 
of volume hv^, ... The total fiux of A out of all these 
elements will be the outward flux across S\ because over all 
boundaries which are not parts of tbe flux will be taken 
twice in opposite directions. 

Hence if 8^. denotes the boundary surface of we have 

^A^d8 = 


n 


= S (div ASv^ + e^Sv,), 

r*=l 


where 0 when Sv,. 0. The latter step follows from the 
definition of div A as a limit. 

Therefore, by increasing n so that the Sv’s tend to zero, 

we get r r 

xA^dS or (AdS)= divAdv .(1). 


1-71. The theorem of 1-7 may be extended to apply to the case 
in which the region R contains a surface S' across which the vector 
field is discontinuous, i.e. there is a sudden change in A in crossing 
the surface S', 

Suppose that S' divides R into two parts R^ and i?, in each of which 
we assume that A is continuous and let ^ 2 n denote the normal 
components of A in R^, R 2 close to S', taken positively towards S' in 
both regions. Then we can apply the theorem of 1*7 to each of the 
regions R ^, R^ separately and add the results. This gives 

/^divA dv= I^A„dS + J^JAi„+At„)dS' .(1). 

Across any surface at which A is continuous, we have, of course, 
ahd the result reduces to 1*7 (1). 

♦ See the author’s Electricity and MagneHam, pp. 8, 9. 
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1*72. Green’s Second Theorem. The theorem of 1*31 may 
be obtained in vector form from 1*7, thus: 

Let A = F grad F', where F, F' are scalar functions which 
with their first and second derivatives are finite and continuous 
through the region bounded by the surface S, 

Then J(A dS) = Jf (grad F'dS), 

and div A = VA = V (FVF') 

= VFVF'+FV2F' 

= grad F grad F' + F div grad F'. 
Therefore 1*7 (1) is equivalent to 


Jgrad F grad V' dv=jv (grad F'dS) — JF div grad V'dv 

.( 1 ), 

or if we use the Hamiltonian operator 


JvFVF'dv = Jf (VF'dS) - jvVW'dv .(2); 


and the symmetry of the left-hand side shews that F and F' 
might be interchanged on the right-hand side. The reader 
will have no difficulty in identifying (1) and (2) with 1*31 (2). 


1*8. Convergence of volume integrals. In the theory of 
attraction and kindred subjects we frequently have to consider 

integrals of the form where dv denotes an element of 


volume and p a scalar function of position, and the region 
of integration R is bounded by a surface S which includes 
the origin 0 from which the radius vector r to a point of dp 
is measured. It follows that the integrand becomes infinite 
within the range of integration, and we have to investigate 
the circumstances under which the integral can be said to 
have a meaning. 

Draw a closed surface 8^ inside S and surrounding O, then 
so far as the subject of integration depends on r it is finite 
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C ^ P dv 

between 8^ and 8, and -— has a definite value, where the 

limits indicate that the integral is 
taken through the volume bounded 
by 8 and 8^. If as 8^ contracts upon 
the point 0 the last integral tends to 
a definite limit which is independent 
of the shape of 8i, then 


Um 



is defined to be the value of the given integral through JB 
and the integral is said to converge. If the integral does not 
converge, it is said to diverge. 

A condition for convergence is that corresponding to any 
small positive number € there must be a surface 8^, defined as 
above, such that for all surfaces 8^ inside 8^ and surrounding 0 



( 1 ). 


For if a limit A exists we can choose 8^ such that 


Js. IJs. 




md therefore 

Js, Js. Js. 


pdv 


<€. 


Conversely, if condition (1) holds, then a limit exists, for 
by taking 8-^ small enough we can make the fluctuation of 
the value of the integral for different cavities inside as 
small as we please and therefore small without limit. 

Now let M denote the greatest value of p within a finite 
region round 0. Let rj denote the greatest value of r on 8^ 
and r 2 the least value of r on S be the sphere of radius 

7j and a the sphere of radius rg. Then 
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but, as in 1*42, w© may put dv = r^dwdr, so that the contribu¬ 
tion to the last integral of the con© of solid angle dco, when 
n < 3, is 


doj I 

J rt 


dr or 


doj 
3 — 71 


(n 


,3—w 


Consequently 


or 


Z-n ■ 

47771 ^^”' 

^ 3-7^ ’ 


and 


C^^pdv 4:7 tM7I^^ 
^ 3-71 • 


Hence, provided ti < 3, for any arbitrarily chosen positive 

number c however small we have only to take 17 less than 
1 

{(3 — 7 i)€/ 47 rJf} 2 ~" in order to get a surface such that for 
all surfaces S 2 inside 8^ 


C^^pdv 

}s,~^ 


<€. 


Consequently the given integral is convergent if < 3. 

W© note that the form of the outer boimdary 8 does not 
come into consideration, because it is only in the vicinity of 
the origin 0 that the integrand tends to become infinite. 

It can be shewn that for ^ 3 the integral is divergent. ♦ 


♦ For further information on the subject see J. G. Leathern’s Tract 
on Volume and Surface Integrals used in Physics, Camb. Univ. Press, 
1906. 
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Chapter II 

GRAVITATIONAL ATTRACTION AND 
POTENTIAL. SIMPLE APPLICATIONS 

2*1, The Law of Gravitation. The law discovered by 
Newton* is as follows: Every particle in the universe attracts 
every other particle with a force which is directly proportional 
to the product of the masses of the particles and inversely propor- 
tional to the square of the distance between them,^ 

Thus if m, m' denote the masses of two particles and r their 

distance apart, the force of attraction between them is y ; 

where y is known as the gravitation constant, and measures 
the attraction of two particles of unit mass at unit distance 
apart. 

It is evident that some limitation must be imposed upon the 
dimensions of the particles in order to avoid a dilB&culty in 
defining the distance between them. For this purpose we 
may define a material particle as a body so small that, for 
the purposes of our investigation, the distance between its 
different parts may be neglected. 

We shall prove later that a uniform solid sphere attracts 
external particles as though its mass were concentrated at 
its centre; and this may be regarded as some justification for 

* Sir Isaac Newton (1642-1727), Lucasian Professor of Mathematics 
at Cambridge and Fellow of Trinity College. His great work Philo- 
Sophias Naturalis Principia Mathematica was published in 1687 by the 
Royal Society, then under the Presidency of Samuel Pepys. He made 
discoveries in Optics and invented his Methodus Fluxionum which 
constituted the foundation of the Calculus. In later life he left 
Cambridge for London, becoming President of the Royal Society and 
Master of the Mint, 

t An outline of the process of reasoning by deduction, from Kepler’s 
laws of planetary motion through the propositions of Newton’s 
Principia to the law of gravitation, is given in the author’s Dynarmcs, 
Part I, pp. 167-8, 
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the hypothesis we propose to make, namely, that for the 
purpose of calculating attractions we may regard particles 
as quantities of matter concentrated at points. 


2*11. Numerical value of y. An approximate value for y 
may be determined thus: Regarding the earth as a sphere 
of radius r and assuming it to attract a small body on its 
surface as though its mass were collected at its centre, then if 
E denotes the mass of the earth in grammes and M that of 
the body, 

y —^ = weight of the body = Mg .(1). 


So that, if p denotes the mean density of the earth, E — 
and from (1), 




4:TTpr 


.( 2 ). 


But, in c.G.s. units, 

9 = 981 cm. sec.~^ p = 5'67 approx, 
and ^Trr = a quadrant of the earth’s circumference 


= 10^ cm. 


3 X 981 1 

’'“ 8xl0»x5.67 - 15.500,0W) ^PP^^. 

when c.G.s. units are used. 

The mean density of the earth has been determined by a 
series of experiments in which the attraction of the earth on 
a given body has been compared with the attraction on the 
same body of another body of known mass.* 


2*12. Astronomical units. We can avoid the continuous 
repetition of the constant y throughout our work by choosing 
units such that y = 1. Such units are called astronomical or 
theoretical units. 

We observe that the acceleration/produced by the attrac* 
tion of a particle of mass m on a particle at a distance r is 

given by/=y^; so that y = 1 when/, m and r are all unity. 


♦ See Encyd. Brit, article ‘Gravitation*. 






22 ATTRACTION ANB POTENTIAL [2*12- 

Henoe the astronomical unit of mass is the mass of a particle 
which by its attraction produces unit acceleration at unit 
distance. 

We can find the astronomical unit of mass in grammes by 
taking the above formula for acceleration, which holds good 
in all systems of units, and putting r = 1 cm. and/— 1 cm. sec."“^ 
so that 1, or m= l/y= 15,500,000 grammes. 

In what follows we shall omit the constant y, save that in 
the collections of examples from examination papers we shall 
print the questions in the form in which they were set, in 
order that readers may become familiar with the use of y. 

2*2. Field of force. When we speak of the attraction of 
a body or system of particles at a point external to itself, 
we mean the force of attraction which the body or system 
would exert on a particle of unit mass placed at the point. 
There must be a definite value (including zero at a point of 
equilibrium) for this force at every point at which a particle 
can be placed. Thus we arrive at the conception of a field of 
force, or region of space with every point of which there is 
associated a force which is definite in magnitude and direction. 

Inasmuch as a force belongs to the class of physical magni¬ 
tudes known as vectors, the field of attraction of any system 
of bodies is an example of a vector field. 

2*21. Attraction of a system of particles. It is part of 
the Newtonian theory of attrac¬ 
tion that the force exerted at a 
point by a system of particles 
is the vector sum of the forces 
exerted by each of the particles 
separately. This is called the 
principle of superposition of 
fields of force. 

Let particles of masses mj, 

^2 > ^3 > • • • be situated at points 
A^, A^, ... whose co-or¬ 
dinates referred to rectangular axes are c^; a 2 , Cg; 




SYSTEM OF PABTIOLBS 


23 


2*22] 

®8> ^8» •••• I* required to find expressions for the com¬ 

ponents X, Y, Z of the attraction of the system of particles 
at any point P whose co-ordinates are x, y, z. 

Let rj, r^, r^, ... denote the distance PAi, PA^, PA^, .... 
The attraction at P of nii &t A^ is mjr-^ directed along 


PA^-, and PA^ has direction cosines 



h-y 
»•! ’ 


where 


^1 = («i - + (^1 - y)^ + (Ci - 3)*- 


• Therefore the particle exerts at P a force whose com 
ponents parallel to the axes are 


^ ^ _ 


The other j^articles m^, make like contributions, so 

that, by the principle of superposition of fields of force, we 
have 



where 5=1, 2, 3, ... and the summation extends to all the 
attracting particles. 


2'22, Potential. Let us now define a function V by the 
formula 

F = .(1), 

where + .....‘....( 2 ), 

in the notation of 2-21. 

Thus defined, F is a function related to a system of attracting 
particles having a definite value at every point P of space 
external to the particles. It is a function of the co-ordinates 
(a?, j/, z) of P and is clearly a single-valued function, in the 
sense that it cannot have more than one value at each point 
P; for it represents simply the sum of the masses of the separate 
particles divided by their respective distances from P. Further, 
though we have expressed F as a function of a?, y, z, it clearly 
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represents a sum which does not depend on the particular 
system of axes of reference. 

Now, by differentiation, we have 

dx r/ dx ‘ 
dv 

But, from (2), »-g ^ = - (a, - a;), 


SO that 



and by 2*21 (1) this is equal to X. Similarly 


dV_ 

dy 


Y 


and 


dz 


Z. 


(3). 


The function V defined by (1) is called the potential of 
the attracting particles, or the potential of the field of 
force. We have proved that the derivatives of V with 
regard to a;, y, z give the components of attraction at P in 
the directions of the axes. And since the directions of the 
axes can be chosen arbitrarily, it follows that the space- 
derivative of F in any direction gives the component of 
attraction in that direction. This may be verified thus: let 
djds denote a differentiation in a direction ds whose direction 

^z 

cosines are l,m,nor7r,^,^- Then 
ds ds ds 


ds dx ds ^ dy ds^ dz ds 
,dV dV dV 

= 4-5—|-m-5- + »3- 

dx dy dz 
— IX+mY + nZ 

=force component in direction ds .(4). 

In the language of vectors the force of attraction is the 
gradient of the potential (l-S), or, if R denotes the resultant 
attraction at a point at which the potential is V, 

R = gradF 


( 6 ). 
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The potential function V plays an important part in the 
theory of attractions, since, as we have just seen, if the 
potential of any given distribution of matter can be determined 
the attraction at any point can be found immediately by 
differentiating the potential. 


2*23. Physical interpretation of potential. With the 
notation of 2*22 we have, for the total differential of the 
potential, jj, 

= Xdx+ Ydj/+Zdz, 

from 2*22 (3). 

Hence, by integrating along any path from P to Q, we get 


Fn-F, 




dy ^ dz\ 


ds 


.( 1 ). 


But the integral represents the work which the forces of 
attraction would perform upon a particle of unit mass as it 
moved along this path from P to This gives us a measure 
of potential in terms of work per unit mass, in that the potential 
at any point Q exceeds the potential at any other point P by the 
work which the forces of attraction would perform upon a particle 
of unit mass as it moves along any path from P to Q, 

It is clear that the addition of a constant to the potential 
will not affect the values of the force components since these 
are obtained by differentiating the potential. Also, when the 
potential is determined by integration, from known force 
components, as in (1) above, the constant of integration may 
be so chosen as to make the potential vanish at an infinite 
distance from the attracting matter, since this will accord with 
the definition of 2*22 (1). On this hypothesis we see that the 
potential at a given point due to a given attracting system is 
the work that would be done by the oMractions of the system on 
a particle of unit mass as it moves along any path from an infinite 
distance up to the point considered. 

We have seen that the definition of potential as S(m/r) 
leads to this expression for potential in terms of work done 
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per unit mass, and it is easy to demonstrate the converse, thus: 
Let m be the mass of a typical 
particle of the system situated 
at the point A. Let PP' — ds be 
an element of any path from an 
infinite distance to the point Q, 
and let = r, and AP' — r-\-dr. 

Then so far as the field of force 
depends upon the particle m at 
A its value at P is rrijr^ directed along PA\ and the work 
done by this force on the unit particle as it moves from P to 



P' is 


m 




ds, or 


— mdr 


Hence the total work done by 


the attraction of the particle m on a unit particle moving 
from an infinite distance to the point Q is 




-xdr 

y.2 


or 


m 

AQ' 


By the principle of superposition of fields, the total work 
done by the attractions of all the particles of the system is 
obtained by adding for their separate effects, so that 


F = S 


m 

AQ 


gives the potential at Q\ where the formula represents the 
sum of the masses of the separate particles each divided by 
its distance from Q, The interchangeabihty of the two defini¬ 
tions of potential is thus completely established. 


2*24. Dimensions. OraviMi(malpotentialmwtnot\)eoori- 
fused with potential energy. They have di^erent physical 
dimensions; because the dimensions of potential energy are 
those of work, i.e. in terms of the fundamental units 

of mass, space and time; but the dimensions of gravitatioruil 
potential as explained in 2*23 are those of work per unit mass, 
i.e. L2T-2, or the square of a velocity. 

In this connection it is important to remember that we are 
using astronomical units (2*12) and omitting the gravitational 
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constant y, and though this does not aflFect the argument 
when potential is defined as work per unit mass, yet when we 
use the formula for potential of 2*22 (1) if we want to find the 
dimensions of V we must restore the constant y and write 

F = yS (m/r), because y has dimensions. By definition 

represents a force and is therefore of dimensions MLT“^, so 
that the dimensions of y are and the dimensions 

of y— are as above. 

^ r 

It is important also to note that though potential energy 
decreases when work is done, gravitational potential increases. 

2-25, The definition of potential in terms of work (2*23) leads at 
once to the relation 

R = grad F 

of 2*22. Thus if P, P' are any two neighbouring points at a distance da 
apart, F, F + dF denote the potentials at P, P' and P the component 
of attraction at P in direction PP', we have 

(F + dF) — F = increment in potential from P to P' 

= work done hy attractions on xmit mass passing from 
P to P' 

^Fda. 

dY 

Therefore F^ — , and this being true for components of attraction 
in every direction round P it follows that 

R = grad F. 

2-26. Equipotential surfaces and Lines of force. Re¬ 
garding the potential F of a given attracting system as a 
function of co-ordinates x, y, z, the equation 

F {Xy y, z) = const, 

represents a surface over which the potential is constant. 
Such surfaces are called eqpiipotential surfaces. It is 
evident from the definition of potential that only one such 
surface passes through any point of space, so that no two 
equipotenitial surfaces can intersect. Also since F has a con¬ 
stant value over such a surface no work would be done by 
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the attractions on a particle moving on such a surface, there¬ 
fore at every point the resultant attraction is normal to the 
equipotential surface through the point. Thiels also obvious 
from the relation R = grad V (1*5). 

A hne so drawn that its tangent at each point is in the 
direction of the resultant attraction at the point is called 
a line of force. It is evident that a hne of force is at right 
angles to the equipotential surfaces at all their points of 
intersection; and, conversely, a surface which cuts all lines 
of force at right angles must be an equipotential surface, 
because at no point on the surface is there a component of 
force tangential to the surface, so that no work would be done 
on a particle moving on the surface and there could therefore 
be no variation in the potential. 


2*3, Continuous bodies. We have now to pass from the 
attraction components and potential of a system of separate 
particles to the attraction components and potential of 
distributions of matter regarded as continuous bodies. We 
remark in the first place that such bodies are not reaUy 
continuous but consist of separate molecules and ultimately 
of indivisible electrons. By the principle of superposition 
therefore there should be no difficulty in obtaining the attrac¬ 
tion components and potential of such a body provided that 
we have a means of summing the contributions of the separate 
particles. It is natural to look to integration to efiFect the 
summation, but when we consider what the process of inte¬ 
gration involves, we find that it does not fit the physical 
conditions of the problem precisely, and it is only by giving 
a special interpretation to our conception of ‘body’ that we 
can jjpstify the use of integration. Thus it is usual to represent 
the potential of a continuous body by a volume integral 



where dv denotes an element of volume of the body at 


a distance r from an external point P and p denotes the density 
of matter in dv. The process of integration implies that the 
density of the body is continuous; but this is not the case 
for a real body with a molecular structure. To justify the 
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use of such an integral it is necessary therefore to suppose 
that it is applied to a hypothetical continuoua distribution of 
matter occupying the same region as the body and having 
at each point a suitably chosen density; this density being 
found by considering a small but finite volume surrounding 
the point and taking the average through this small volume 
of the masses of the particles of the real body contained therein. 


2*4. Before proceeding further with the general theory of 
the subject we shall calculate the attraction components and 
the potential of some simple bodies. 

Attraction of a uniform straight rod. Let m denote the 
mass per unit length of a uniform rod AB, It is required to 
find the components of attraction of ^ 5 at an external point P. 



Let the perpendicular from P to AB meet AB in M, which 
for simplicity we take on produced though this is not 
necessary for the argument. Let MP=p. Consider an element 
QQ' of the rod, where MQ==^x and QQ' = dxy and let the angle 
MPQ be 0. Then x^p tan 0, and the mass of the element QQ' 
of the rod is mda; = mp sec^ 0d0, 

The attraction at P of this element is therefore 


or 


mp sec^ ddd 

PQ^ 


along PQ, 


mp sec^ Odd ^mdO 
p^ ^ec® 0 ^ p 


along PQ 


( 1 ). 


Denoting the components of attraction of the rod parallel 
and perpendicular to BA by X and 7, and the angles MPAy 




30 ATTEAOTION AND POTENTIAL [ 2 ' 4 - 

MPB by a, we have 


X=f —sin0(i!0=—(cos8—cosa) 

J/3i> ' 

= —8ini(a + ^)sini(a-^) .(2) 

Ir 

and F=f — costfd0 = —(sina—sinfl) 

hP V 

= ^8m|(a-^)coa|(a + ^) .(3). 

The resultant attraction is therefore ^sin|(a-~j8), or 

sin \APB, and it makes an angle tan"^ “ or J (a + jS) with 

PM, i.e. it acts along the bisector of the angle APB, 

Another convenient form for the component parallel to the 
rod is 

X = -(co8^-co8a) = ^-^.(4) 


in the sense parallel to BA, 

We note that if a circle of centre P and radius PM cuts 
PA, PQ', PQ, PB in D, R\ R, E, then the attraction at P of 
the element RR' of a rod in the form of a circular arc DE of the 
same line density (mass per unit length) as JIJ? is 


mpdd ^mdd 
~~¥ 


= attraction of QQ!. 


Hence this circular arc DE exerts the same attraction as 
the rod AB, 

Cob. If the rod is infinitely long, the angle APB is two 
right angles and the resultant attraction is Imjp perpendicular 
to the rod. It would appear from this result that if thef 
attracted particle were close to the rod the attraction would] 
be infinite; but this conclusion is not justified because in the 
foregoing argument we assumed that every point of an element 
QQ' of the rod was at the same distance from the point P, 
and for this to be true when P is close to the rod it would be 
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necessary for the rod to have no thickness. To find the attrac¬ 
tion at a point close to a rod of finite thickness it will be 
necessary to take account of the cross-section of the rod. 
The attraction of a uniform long rod whose cross-section is 
a circle is found in 2*42 below. 


2*41. Potential of a uniform straight rod. With the 
notation of 2-4 the potential at P is given by 

F— r 

]yQ 

= m j^log tan \{6 + ^tt) J 

— mlogQot \PAB cot \P BA .(1). 

If 21 denotes the length of the rod AB and PA = r, PB = r', 
we easily find from (1) that* 

jT . I* + 2? 

V = mlog- 




.(2).y 


It follows that if .4, J5 are foci of an ellipse passing through 
P and 2a is its major axis, 

.(3), v/ 


T7 ,l-f“C . 

or F=mlogj-^ .(4), 

where e denotes the eccentricity of the ellipse. 

Hence the potential is constant over any prolate spherokj . 
of which J., B are the foci; i.e. this family of confocal prolate 
spheroids are the equipotential surfaces. Since the normal 
to an ellipse at any point bisects the angle between the focal 

* Iff + r' + 2Z=2d, 

cot iPAB cot iPBA = ^ {(« _ t) (» -21)' (a - r') (a - Si)} 

_ 8 _ r-^r' + 2l 

r' —2r 
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distances and the resultant attraction at a point is normal 
to the equipotential surface, it follows that the resultant 
attraction at P bisects the angle APB, as we found in 2*4. 

If the rod be of great length and P in the neighbourhood of 
its centre, then in (2) we may put r + r' = 2 where p 

is small compared with 1. This gives 


V 


= mlog 




= 2m log 


V(P±^±l 

P 


= 2mlog{(2i+g)/^>j, 

or neglecting 

V = 2m log 21 — 2m logp 


(5). 


By differentiating the last result we get for the attraction 
of the rod in the direction p increasing, agreeing 

with 2*4 Cor.; and conversely, if we put ^ and 

integrate, we find for the potential V «= (7 — 2mlogp, in agree¬ 
ment with (5). 


2-42. Attraction and potential of a long thin uniform cir¬ 
cular cylinder. Take a cross-section of the cylinder about the middle 
of its length. Let 0 be the 
centre of the cross-section and 
P any point inside it. Through 
P draw chords QPR, Q'PR' 
making a small angle d6 with 
one another, and intercepting 
small arcs QQ\ RR' on the 
circle. We can conceive the 
cylinder to be composed of 
long parallel rods, and take 
QQ', RR' as the cross-sections 
of two of them. Then if m 
denotes the surface density (mass per tmit area) of the cylinder, mQQ' 
and mRR' denote the mass per unit length of the two rods. Let the 
equal angles OQP, ORP be denoted by Then using the result of 
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2*4 Cob. for the attraction of a long rod, we have 
attraction at P due to rod through QQ' 

= = 2m sec <j> dO along PQ, 

attraction at P due to rod through BR' 

~ = 2 m sec d6 along PR, 

Hence these two rods exert equal and opposite attractions at P; 
and by dividing up the whole cylinder into similar pairs of rods we 
see that its resultant attraction at any internal point is zero. 

It follows that the potential must be constant at all points inside 
the cylinder sufficiently far from its ends, and is therefore equal to 
the potential at O, and by 2*41 (5) the potential due to the rod QQ' is 
2mQQ' log 21 ~ 2mQQ' log OQ, 

where 21 is the length of the cylinder. Therefore, for the whole potential, 
F = 2 Mlog 2 Z- 2 Mloga .( 1 ), 

where M = 27rma, a being the radius. 

Again, for an external point, take P' the inverse of P, so that 
OP. OP' = and, by similar triangles, OP'Q = OQP = = OP'B, Then 
the attraction at P' of the rod through QQ' 

2mQQ' 2mPQdd8ec<l> « j ® i r>//^ 

“”P^^^ ^■ OP' ^ ®• 

But the resultant attraction of the cylinder at P' is clearly along 
P'0, and by resolving the attraction of the rod through QQ' ia this 

direction we get dd. Similarly, other rods give like results and the 

whole attraction at P' is therefore ^^7 = ; %nd we observe that 

this is the same as if the whole mass of the cylinder were condensed 
into a rod of equal mass along the axis of the cylinder. 

To find the potential, we may put r for OP' and write 

dr”' r * 

so that V = const. — 2M log r, 

where, in order that V may take the form ( 1 ) when r = o, (3-7), the 
constant is 2 ikf log 21, so that 

F = 2ilflog2Z-2Mlogr. 

It is clear that the same results will be true for the attraction and 
potential of a tmiform long solid cylinder at external points, for the 
cylinder can be divided up into coaxial thin cylinders for each of which 
the results are true. 


RN A 


3 
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2*5. Attraction of a uniform circular disc at a point on 
its axis. Let 0 be the centre of the disc, P a point on its 
axis Oz at a distance z from O, and 
let m denote the surface density or 
mass per unit area. Divide the disc 
into concentric rings. Let OQ=^x be 
the radius and QQ' = dx the breadth 
of one of these rings. The mass of the 
ring is 2TTmx dx, and the attraction at 
P of each element is got by dividing 
the mass of the element by PQ^. But 
the resultant attraction of the ring is clearly in the direction 
PO, so that its magnitude is 



27Tmxdx 

p¥ 


cos 6, 


where d is the angle OPQ, But x^ztaxid, so if a is the angle 
which a radius of the disc subtends at P, we have for the whole 
attraction of the disc 


J, 


“ tan 6 sec^ 0 cos 0 

z^ sec^O 


d0 = 27rm(l —cosa) ...(1). 


For an infinite plate we may put a == so that the attraction 
of an infinite plate is 27rm at right angles to itself. But we 
majr not infer that the attraction of a circular disc at a point 
on its axis close to itself is 27rm unless the thickness of the 
plate is negligible, so that the particles of a ring in the 
immediate vicinity of the attracted particle may be regarded 
as equidistant from it. 

A frustum of a right circular cone bounded by planes at 
right angles to its axis may be regarded as a pile of discs all 
subtending the same angle at the vertex, and the attraction 
at the vertex of the cone due to such a frustum is easily seen 
to be 27Tm^(l —cosa), where t denotes the thickness of the 
frustum, m the mass of unit volume, and a the semi-angle 
of the cone. 
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2*51. Potential at any point on the axis of a uniform 
circular disc. With the notation of 2*5, we have 


-j: 


2TTmxdx 

PQ 


, where a is the radius of the disc, 


= 27nn 


i: 


xdx 


^y{z^+x^) 


= 2TTm{^/{z^ + a^) — z} .( 1 ). 


We note that this gives for the attraction in direction PO 


dV 

dz 


= 27rm {1 — 




in agreement with 2*5 (1). But if we find the potential from 
the attraction by integrating the last formula, we get 

V = 27rm — + 

where (7 is a constant of integration. In the physical definition 
of potential in terms of work, it is assumed that the potential 
is zero at an infinite distance from the attracting matter. 
And if in the case under consideration V is to vanish when 
5 ; 00 , we must have (7 = 0. 

If zja is small, then by (1) 

F = (7 — 2TTmz, .. (2), 


where C — 2'nma, a constant. 

It should be noted that the formula (1) is equivalent to 
F = 27rm {BP — OP) [fig, 2-5], and that this will give the value 
of F on either side of the disc if 8P, OP denote numerical 
lengths. And for negative values of z, instead of (1) we 
should hare F = 2^{V(^*+a<*) + 3 } .(3). 

The expressions (1) and (3) are both comprised in the 
formula V = 27 rm{V{z^ + a^)-\z\]. 


2*52» Attraction at a point on the axis of a uniform circular 
cylinder. Let I be the length, a the radius and p the density of the 
cylinder. Let the point 0 at which the attraction is to be foimd be at 
a distance c from the nearer end. The cylinder may be regarded as 


3-2 
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a pile of discs, and a disc of thickness dx at a distance x from O has 
a mass pdx per unit area, so tliat its attraction is 

~ -/(»* +a»)} 

and the whole attraction 

rz+c [ £r 1 , 

= 2wpj^ + 

= 2 tp {«--/[(« + c) 2 + a2] + -^/(c* + a“)}. 

If the cylinder is of infinite length in one direction 
we have to find the limit of the last expression as 
i -> 00 , i.e. the limit of 

2,rp {«- r (1 + J - y + + a“)} 

= 27rp{-v/(c® + a‘*)-c}. 

If the particle be at the centre of the end of such a cylinder the attrac¬ 
tion is 27rpa, 

2*6. Use of solid angles. There are many applications of 
the theory of attraction in which calculations are simplified 
by the use of the sohd angle. For example, to find the com^ 
ponent of attraction perpendicular to itself produced by a plane 
plate of any form. 

Let dS be an element of area olthe plate at a distance r 
from the point O and subtending a solid angle c?a> at 0. Then 
if m denotes the mass per unit area of the plate, the attraction 
of this element at O is mdSjr^, and 
resolving this at right angles to the 
plate we get mdS cos 9jr^, where 6 
is the inclination of r to the normal. 

But dS cos d/r^ = do) (1’4). Therefore 
mdo} is the contribution of tl^is 
element dS of the plate to its whole 
attraction at 0 at right angles to itself; and the attraction of 
the whole plate in the same direction is mu>, where m is the 
sohd angle which the plate subtends at 0, 

2-61. Bacamples. (i) Shew that all parallel frusta of a uniform cone 
of the same thickness exert equal attractions at the vertex. 

Since a cone of any form can be subdivided into a number of cones 
of small solid angle having the same vertex, and since the attraction of 



3 


as ^ 00 
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a thick frustum is got by adding the attraction of a number of thin 
ones, it is sufficient to prove the theorem for parallel thin slices of 
a cone of small solid angle. 

Let PQ be a slice of the cone of small thickness t, and let p be the 
density, O the vertex, OQ = r and 
6 the angle the normal to PQ 
makes with OQ, Then the area 
PQ = r^8eo6dw (1-4), and the 
mass of the frustum = prH sec Odco, 

The attraction at O of this 
frustum is therefore ptsecddcuy 
and this has the same value for all parallel frusta of the same thickness, 
.. . 

(li) Prove that the potential of a solid of uniform density p at an external 

point P can be represented by a surface integral ipj cos 0dS over the surface 

of the solids where Q is the angle between the inward normal to dS and the 
line joining dS to P. 



Lot a cone of small solid angle do} 
and vertex P cut the surface of the solid 
in elements of area dS ^, dS^ at .4, B, 
where the inward normals make angles 
, $2 with the line BAP, Let 

AP, BP = ri, rj. 

The mass of an element of volume of 
the cone at a distance r from P is, as in 
1*42, pr^dwdr. Hence the mass of the 
cone between A and B produces at P a 

= cos ^2 + ^i) 



If we take the sum for all cones which intersect the solid we shall get 
Jp fcos BdS integrated over the surface. 


EXAMPLES 

1 . A sphere of 150 kilogrammes, placed with its centre 30 cm. 

vertically below that of another sphere, is found to cause an apparent 
increase of the latter’s weight by the 1*14/10® part. What value does 
this imply for the constant of gravitation? [M. T. 1918] 

2. Three equal uniform rods, each of length 2a and of mass m per 
unit length, form an equilateral triangle. Shew that the attraction 
on a particle of unit mass at a point on the straight line through the 
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centre of the triangle at right angles to its plane, and at a distance z 
from it, is 

[C. 1906] 






3. Prove that any plane closed rectilinear polygon of thin uniform 
attracting wire, all of whose sides touch one and the same circle, 
exerts no attraction at the centre of the circle. 


[London Univ. 1927] 


4. A imiform wire of mass M has the form of a semicircular arc 
of radius a together with the diameter of the semicircle. Find the 
attraction of the wire upon a particle of unit mass placed at the other 
end of the diameter through the middle point of the arc of the semi¬ 
circle. [London Univ. 1933] 


5. Prove that the attraction of a thin imiform cylindrical shell of 
radius a and lengt-h 2 at a point on its axis at a distance 6 from one 
end and Z — 6 from the other {h<\l) is 

MV I 11 

^ ^ A/K+(i-6m’ 


where M is the mass of the shell. 

[The ends of the shell are open and circular.] 

[London Univ. 1926] 


6. Shew that the attraction of a imiform solid hemisphere at the 

centre of its plane base is where M is its mass and a is its radius. 

7. Prove that the potential of a circle of uniform surface density a 
and radius a at a point on its circumference is 4yo-a. » 

[London Univ. 1932] 


8. Prove that the potential of a plane annulus bounded by concentric 
circles at a point on its axis is 27ryp (V — l)i where p is the mass per unit 
area and Z, V are the distances of the point from the inner and outer 
edges. 

9. A uniform circular lamina is divided into two segments by 
a chord which subtends an angle 2a (< tt) at the centre. Prove that 
the attraction at the centre of the circle due to the minor segment is 

2ycr [log tan ( Jtt + ioc) — sin a], 

a being the surface density of the lamina. [London Univ. 1931] 

10. Prove that, if the surface density of a circular disc at any point 
distant r from the centre is Ar, then the attraction at a point on the 
axis at a distance d is 

'v/(a®4-d2)-f-a 2a 1 

[London Univ. 1932] 
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11. Four uniform rods, each of the same mass per \mit length, 
form a square A BCD, Any point P is taken inside the square and 
E, F, G, H are the feet of the perpendiculars from P upon AB^ BC, 
CDf DA respectively. Prove that the potential at P due to the square 
A BOD is the same as that due to four rods EF, FG, GH, HE each of 
the same mass per unit length as that of the rods forming the square. 

[London Univ. 1935] 

12. Two imiform rods AB and (7Z), of length a and mass m per unit 
length, are placed so that CD lies along a right bisector of AP, and G 
and D are distant a and 2a respectively from the rod AB. Shew that 
the attraction of either rod on the other is of magnitude 

2ym’^{log,(l+ V5) + log,2-logJl+ 

[London Univ. 1938] 

13. Find the attraction of a imiform circular ring, of mass M 
and radius a, upon a thin uniform rod, of mass Af' and length 2^, 
lying along the axis of the ring with its middle point at a distance 
x(<l) from the centre of the ring. 

If the ring is fixed and the rod released from rest when its centre is 
at a small distance x from the centre of the ring and moves under the 
mutual attraction only, shew that the period of small oscillation is 

27t{a^-\-l^)^l(yM)^. [London Univ. 1935] 

14. Find the attraction of a uniform thin shell, in the form of the 

portion of a circular cylinder between two planes perpendicular to 
the axis, on a particle at the axis; and prove that a particle slightly 
displaced on the axis from the position of equilibrium will perform small 
oscillations of period 27rAf(a^ + c®)*, where 2a is the length eind c 
the radius of the cylinder, and M its mass. [C. 1892 

16. Prove tiiat the potential due to a thin rod AB, of uniform line 
density m, at any point P in the line A B produced is ymlog (PA/PB). 
Two such rods, each of length 2a, are plaoed so that they lie in the 
same straight line with their nearer ends at a distance 2a apart. 
Prove that a particle, released from rest at a point in the line of the 
rods at a small distance x from the point midway between them, begins 
to move with acceleration 16yma?/9a*, approximately, f 

[London Univ. 1936 

16. Prove that the resultant attraction at P of a straight rod AB 
of uniform line density m is 2ymelPG diropted along PGy where e is 
the eccentricity of the ellipse of foci A, B passing through P, and G 
is the point in which the normal to the ellipse at P cuts AB. 

17. The ends of a thin uniform rod are at the foci of a hyperbola, 

and a similar rod is placed parallel to it with its ends on the curve. 
Shew that the attraction between the two rods is inversely proportional 
to their distance apart. [C. 1908] 



40 


ATTEACTION AND POTENTIAL 


18. Shew that the potential of a closed hemispherical shell, made 
of uniform thin material of surface density cr, at a point on the axis 
of symmetry at a distance a from its plane face remote from its ctirved 
surface is 27 Tyaa, where a is the radius of hemisphere. ^ 

[London Univ. 1937] 

19. ABy A'B' are two parallel rods of line densities m, rti' and are 
placed with the line joining their middle points perpendicular to the 
rods and of length c. Shew that the attraction between the rods is 

^^(AB'-AA'). [C. 1014] 

c 

20. Two uniform rods AB, AC, each of mass m and length a, are 
freely jointed at A, Shew that the couple required to keep them at 
right angles against their mutual attraction is (2 — y" 2 ) ym®/a. 

[London Univ.] 

21. Two uniform thin rods AB and CD are of linear density A 
and A' and length 2a and 2c respectively. The midpoint E of AB and 
the midpoint F of CD are at a distance b apart, and AB, CD and EF 
are mutually perpendicular. Shew that the attraction between the 
rods is 

4yAA' sm-1 {((, 84 . 53 ) ( 6 a + c!>)}i ’ 

where y is the gravitational constant. [P. 1934] 

22. Shew that the attraction of a uniform solid circular cylinder, 
of radius a and density p, at an external point P on the axis of the 

cylinder is 27 rypa (tan Ja — tan J/S), 

where a, jS are the angles subtended at P by the radii of the nearer and 
farther ends of the cylinder. [London Univ. 1926] 

23. A finite solid, of uniform density p, is bounded by the paraboloid 
a:* -f 2/2 = 4 a 2 : and the plane z = a. Prove that the intensity of attraction 
at the origin is 

2nypa [1 - -v /6 + 4 log, J {1 + Vs)]. 

[London Univ. 1938] 

24. A quantity of homogeneous matter of density p is in the form 
of a portion of a paraboloid of revolution boimded by a plane perpen¬ 
dicular to the axis. Prove that its attraction at its vertex tends, as c 
the length of the axis increases, to the limit 47 rpalogc/ae, where 4a 
is the length of the latus rectum of the generating parabola. 

26. A uniform solid, of density p, is in the form of the finite portion 
of a paraboloid of revolution of latus rectum 4 a cut off by a plane 
perpendicular to the axis at a distance 2 a from the vertex. Prove that 
the attraction at the centre of the circular end is 

47 r)/a{-\/ 2 ~log(l + \/2)}. 


[London Univ. 1933] 
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26. Shew that for a long rod of length 21 and line density m, the 
component attractions parallel and perpendicular to the rod are 
approximately 

“0 ■T('-fi)- 

where the origin is at the middle point of the rod and the axis of x 
along it, and terms of order equal to or higher than (djl)^ are neglected, 
where d is the distance from the middle point. [M. T. 1919] 

27. A homogeneous prism, infinite in length, whose cross-section 
is an equilateral triangle ABC attracts a particle at A; prove that the 
resultant attraction is 4,TrMISa, where M is the mass of a unit length 
of the prism and a is the length of a side of the triangle. [C. 1908*] 

28. ABC m a, normal cross-section of a uniform triangular prism 

of infinite length and of mass m per imit length. Shew that the resultant 
attraction at the point A is compounded of ^tymAja perpendicular 
to BG, and 4ym log {blc)la parallel to BG, the angle A being measured 
in radians. [London Univ. 1925] 

29. Two infinitely long imiform circular cylinders, each of mass m 
per unit length, are placed with their axes parallel. Prove that the 
potential, F, at a point, outside both cylinders, whose distances from 
their axes are r, r\ is given by 

V = const. — 2ym log (rr'). 

Further, prove that the lines of force outside the cylinders are arcs of 
rectangular hyperbolas. , [London Univ. 1936] 

30. A tetrahedron is made from a uniform sheet of thin metal. 
Shew that, if the law of attraction were the inverse cube, a particle 
would be in equilibrium if placed at the centre of the inscribed sphere. 

31. Two thin straight uniform rods AB, GD are pivoted together 
at their middle points. Prove that the attraction between them 
reduces to a couple of moment 2mm' (AGBG) cosec a, where m and 
m' are their line densities and a is the angle between them. [C. 1903] 

32. The cross-section of a long uniform solid cylinder is a semicircle. 

Shew that the direction of the resultant attraction at a point on an 
edge far from the ends of the cylinder makes an angle tan“^( 2 / 7 r) 
with the plane face. [C. 1928] 

33. Shew that a uniform thin square lamina of mass M, whose 
side is of length Z, exerts an attrewtion at a point P on the perj>en- 
dicular to the lamina through one corner, whose component parallel 
to the lamina is 

A{i + (fe“ + 2Z*)‘} J’ 

where h is the distance from P to the lamina. [London Univ. 1931] 
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34. A plane sheet of matter of uniform suiface density is bounded 
by two infinite parallel straight edges and extends to infinity in both 
directions. Find the components of the attraction parallel to and at 
right angles to its plane; and shew that the loci of the points for which 
these are constant are two sets of orthogonal cylinders. [C. 1898] 

36. Shew that, if a uniform straight rod CD were free to slide in 
a smooth groove parallel to a fixed similar rod A B, its time of oscillation 
about its equilibrium position would be 

where m is the mass of imit length of the rod AB, and AD, BD are 
distances in the equilibrium position. [C. 1926] 

" 36. Determine the potential of a solid uniform right circular cone 
at the centre of its base. [London Univ. 1931] 

37. A vertical solid right circular cylinder of height h, radius a 
and density p, bounded by plane ends perpendicular to the axis, is 
divided by a plane through the axis into two parts. Shew that the 
horizontal attraction of either part at the centre of the base is 

2phlog[{a-{-V(a^ + h^)}lh]. 

38. Shew that two finite surfaces having the same surface density, 

wliioh axe geometrically inverse to one another, exert equal attractions 
at the centre of inversion. [C. 1906] 

39. If the law of attraction of a thin uniform straight rod is that of 
the inverse cube, shew that the equipotehtial surfaces are generated 
by the revolution of a family of curves, whose polar equation is 

(r* -f- cos 26)^ = 2ar sin B . cosec (cr sin B) 

about the initial line, 2a being the length of the rod and c a parameter. 

[C. 1907] 

40. The mass per unit length at a point P of a rod AC B, of which A 

and B are the ends and <7 a particular point, varies as CPjViAP . PB), 
Prove that the resultantjtoaction of the rod at any point at which 
AC, CB subtend equal acts through C. [M. T. 1902] 

41. An infinitely long thin, rod, of uniform line density w, and a thin 
uniform rod of mass M and finite length 2a are so placed that the 
shortest distance between them is of length h and bisects the finite 
rod. Prove that the resultant attraction of either rod on the other 
is equivalent to a wrench of intensity 

2yMm<l>l(a sin B) 

and pitch b cot B (tan ^ l<l>, 

where B is the angle between the rods and ^ = tan“^ [(asin B)/b]» 

[London Univ. 1938] 
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42. A solid uniform hemisphere is of mass M and radius a. Prove 
that its attraction at any point P of its axis on the convex side of the 
curved surface is equal to 

where r and R are the distances of P from the centre and rim of the 
base. [C. 1902] 


Mf^ . i2« + 2r«~ 


43. Shew that, if a denotes the radius and p the density of a hemi¬ 
spherical hill and R the radius and a the mean density of the earth, 
the difference of the latitudes as measured by a plumb line of stations 
on the north and south sides of the base of the hill is increased by the 
attraction of the hill by approximately ap/Ra, 

44. Consider the effect on the plumb line of a deep crevasse running 
east and west, and shew that if the density of the earth’s crust is 
greater than two-thirds of the earth’s mean density the north edge 
of the crevasse will have a lower apparent latitude than the south edge. > 


46. Prove that at the edge of a very long straight depression of 
rectangular section, whose depth h is small compared with its breadth 
c, the plumb line will be deflected through an angle approximately 
equal to 2p,h{l + \ogclh)lg from its mean position, where p is the 
attraction at unit distance of unit volume of the superficial strata 
supposed condensed to a point. [M. T. 1887] 

. 46. The surface density of a circular disc of radius a is cto (1 — r^/a*)* 

at a distance r from the centre of the disc. Prove that the potential 
at a point on the axis of the disc at a distance a cot from its centre is 

irya^a (<)i cosec* ~ cot <j>), 

[London Univ. 1927] 


47. Prove that a lamina, of uniform density a in the form of the 
segment of a circle, attracts a point P on the remaining arc of the circle 
with a force tending towards D whose magnitude is 


|{DPl<,g 


DP — 2Z>aI 

DP-DA 


where AB is the base of the segment, D the middle point of its arc and 
a the radius of the circle. [C. 1906] 


48. Prove that a solid regular tetrahedron of uniform density p 
exerts at a vertex an attraction of intensity 

yph (3 cos“^ i — 

where h is the perpendicular distance of a vertex from the opposite 
face. [London Univ, 1936] 
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49. Prove that a thin triangular lamina ABC oi uniform surface 
density a exerts, at a point on the line through A perpendicular to 
the plane of the triangle and at a distance k from A, an attraction the 
intensity of whose component perpendicular to the lamina is 


ya 


A — sin“i 


h cos B 


— sin“^ 


A; cos C 

VW^)s 


where p is the perpendicular distance of A from BG, 

[London Univ. 1939] 


50. If a \miform lamina is the portion of the plane xy containing the 
origin and bounded by the hyperbola —1, shew that 
the z component of its attraction at any point on the ellipse 

+ b^) -f z^/b^ = 1, 2 / = 0 is , where h is the z co-ordinate 

of the point. [C. 1900] 



Chapter III 

ATTRACTION AND POTENTIAL AT 
INTERNAL POINTS. SPHERES 


3*1. Attraction and Potential at internal points. In 

the previous chapter we confined our attention to the attrac¬ 
tion and potential at points external to the attracting matter. 
We have now to consider what meaning to attach to the same 
terms at points inside the attracting matter. We have no 
reason to suppose that the law of inverse squares holds good 
at small distances comparable with the size of the molecules, 
and we have explained in 2*3 what, for mathematical con¬ 
venience, we understand by the term ‘continuous body’. 
Our definitions of attraction and potential at an external 
point imply the existence of a separate attracted particle at 
the point under consideration. Such a particle cannot exist 
inside a continuous body because two particles cannot 
occupy the same space simultaneously. We therefore imagine 
that there is a small cavity in the body surrounding the 
attracted particle placed at the chosen point. We assume that 
we can calculate the attraction and potential at this point by 
our former rules, since the attracted particle is not in contact 
with the matter. Then we define the attraction and potential 
at the same point in the continuous body to be the limits to 
which the attraction and potential at the point in the cavity 
tend as the cavity decreases in size and ultimately vanishes. 


3*11. Convergent integrals . By reference to the test given 
in 1*8 for the convergence of volume integrals of the form 

J , namely that the integral converges at points within 

the region of integration when rK 3, we see that the integral 

which represents the potential, converges at points 


I 


within the attracting matter; and the same is true of the 
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integrals of the form 

tion components, si; 

shewn, however, that when the attracting matter consists of 
a surface distribution the integrals which represent tangential 
components of attraction at a point in the matter are not 
convergent in the sense of 1*8 but depend on the shape of 
the vanishing cavity. 


/ 


P which represent the attrac- 


ace 


\x — x' 


*.3 


is of order It may be 


3-2. Uniform thin spherical shell. Let m be the mass 
per unit area of a thin spherical shell of radius a and centre O, 
Attraction at an internal point. With any internal point P 
as vertex construct a cone of small solid angle dto intersecting 
the surface in elements QQ\ 

RR of areas dS, dS\ The 
attractions at P of these ele¬ 
ments are mdSjQP^ and 
mdS'/RP^ in opposite direc¬ 
tions. But 

dS^QP^sec OQPdco, 

and 

dS' = sec ORPdo), 

and the angles OQP, ORP are equal, so that the elements 
QQ\ RR' of the shell exert equal and opposite attractions at 
P. Since the whole shell can be divided into similar pairs of 
elements by taking cones in all directions round P, it follows 
that the resultant attraction of the shell at P is zero.* 
Attraction at an external point. Let P' be any external 
point and P its inverse, so that OP .OP' — a^. The resultant 
attraction at P' is by symmetry along P'O; and with the same 
construction, the element QQ' exerts an attraction mdSjQP'^ 
along P'Q, and resolving this in direction P'O we get 

mQP^ sec OQPdco 
--COS OP'Q. 



* This proof was given by Newton, Principia, Proposition LXX. 
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But by similar triangles QP : QP' = OQ : OP', and the angles 
OQP, OP'Q are equal. Therefore the element QQ' contributes 
an amount ma^dcDjOP'^ to the resultant attraction. By taking 
cones in all directions round P, we get for the attraction at P' 
of the whole shell 4:7rma^lOP'^, or MjOP'^, where M denotes 
the mass of the shell. It follows that the attraction of the 
shell at external points is the same as if its mass were collected 
into a particle at its centre. This result is also due to Newton. 


3*21. Potential at an internal point. Since the attraction 
is zero throughout the interior of the shell, there can be no 
variation in the potential, or the potential is constant. The 
potential at every point in the interior is therefore the same 
as the potential at the centre, i.e. Mja, where M denotes the 
whole mass, since every element of M is at the same distance 
a from 0. 

This result may also be obtained directly without assuming 
a knowledge of the attraction, thus: The contribution to the 
potential at P of the elements QQ\ RR' is 

md8IQP-\-mdS'jRP 

= mQP sec OQPdw + mRP sec ORPdu>, 


and since the angles OQP, ORP are equal, this 
= mQR sec OQPdco = 2mada). 

If we take all such double cones round P we take the shell 
twice over, we therefore put 27t for do), and get for the whole 
potential 4:7rma, or Mja, 

Potential at an external point. Let OP' = r, then since the 
force at distance r is Mjr^ in the direction in which r decreases, 
therefore 


dy 

dr 


M 
■ ^2 


and 



+ 0 . 


But the potential vanishes at an infinite distance, therefore 
0^0 and V^^Mjr, 


♦ Frincipia, Proposition LXXI. 
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This result may also be obtained independently, thus: The 
contribution of the elements QQ\ RR to the potential at P' 


= md8IQP'^mdS'IRP^ 

_ mQP^ sec OQPdoj mRP^ sec ORPdo) 

Wp'~~ 


— mdo) 


( 


madu) 


^ sec OgP + 
.QR sec OQP = 


aRP 

'OR 


sec ORP 


j, by similar triangles, 


2ma^dco 

~~OP^^ 


and, as before, for the contributions of the whole shell we 
put 271 for do), and obtain 

_ ^irma^ _ M 


3*3. Attraction of a uniform solid sphere. Such a sphere 
may be regarded as composed of a series of concentric thin 
spherical shells, and the required results may be obtained by 
summation. Let a be the radius and p the density of the 
sphere. 

Attraction at an internal point. Take ^ point P at a distance r 
from the centre. Imagine a thin shell of matter of external 
and internal radii r -f € and r — e to 
be removed and consider the at¬ 
traction at a point P in this cavity. 

By 3*2 the concentric shells external 
to the cavity exert no attraction at 
P, and those internal to the cavity 
attract as though their masses were 
concentrated at the centre 0. Hence 
the attraction at P is the limit as 
€ ~> 0 of ^rrp {r —€)^lr^ = ^7Tpr; i.e. 
the attraction of a uniform solid 
sphere at an internal point is directly proportional to the 
distance from the centre. 

Attraction at an external point. Since each of the concentric 
shells attracts at an external point as though its mass were 
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collected at its centre, the same is true of the solid sphere, 
and its attraction is represented by Mjr^, where M is its mass 
and r the distance of an external point from the centre. 


3*31. Potential of a uniform solid sphere. 

At an internal point. Adopting the method of 3*3, let R 
denote the radius of a shell external to the cavity. Its mass 
is 4:7rpR^dR, so that by 3*21 the potential it produces at a 
point P inside itself is ^TrpRdR, and the potential at P due 
to all such external shells is given by 


r 4:7TpRdR = 2'jTp{a^ —+ 
J r+c 


Also, by 3*21, the shells of radius less than that of the cavity 
produce the same potential as if the mass were collected at 0, 
i.e. f7T/) (r —€)^/(r —€). Hence the whole potential at P is the 
limit as € -> 0 of 


f77/>{(r - €)2 +1 [o* - (r + e)®]} = fwp (3a* - r*) .. .(1). 

At an external point. Since each of the concentric shells 
produces at an external point a potential equal to its mass 
divided by the distance of the point from the centre, the same 
is true for the solid sphere. 

To deduce the potential from the attraction. At an external 
point, as in 3*21, we have leading as before to 


V^Mjr .(2). 

At an internal point, since the attraction is %7rpr, we have 

dV ^ 

SO that F = (7 — ITTpr^ .(3), 

where the constant may be determined from the fact that the 
potential at the surface r~a must be the same whether 
obtained by putting r a in the formula that gives the internal 


RN A 


4 
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or the external potential. But from (2), when r = a, F = 
and substituting in (3) gives C^^irpd^, and 

F = f7rp(3a2~r2), 

as before. 

It is evident that the results obtained for the external field of force 
of a uniform shell bounded by concentric spheres, or for a imiform 
solid sphere, are also true for heterogeneous shells or solids provided 
they can be divided up into thin shells or strata of imiform density. 


3’4. When the potential of any distribution of matter on the 
surface of a sphere is known at 
internal points it can be found at 
external points. Let P, P* be in¬ 
verse points with regard to the 
sphere, so that, if a is the radius, 

OP, OP' — a^. Then, if Q be any 
point on the surface, by similar 
triangles, 

QP\QP'^OQ\OP'^a,OP', 



But, if m denotes the mass of an element at Q, the potentials 
it produces at P, P' are mjPQ and mjP'Q and these are in 
the constant ratio OP': a for all positions of Q, Hence, if 
F, F' denote the potentials at P, P' of the distribution of 

matter on the sphere, we have ^ 


For example, the potential at an internal point of a uniform 
thin shell of mass M and radius a is given by V — Mja (3*21), 
therefore the potential at an external point at a distance r 
from the centre is given by F' = iff /r. 


3*5. Attraction at an internal point of a heterogeneous sphere 
in which the density is a function of the distance from the centre. 
Let a be the radius and let/ (P) denote the density at a distance 
R from the centre. Consider the attraction at a point P whose 
distance from the centre is r. As in 3*2, shells of radius greater 
than r exert no attraction at P and shells of radius less than r 
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attract with a force Mass/r^. Hence it is easy to see that the 
required attraction is 

3«61. Bxeunple. Find the law of density when the attraction at an 
internal point of a sphere of given mass is constant and eqtial to its value 
at the surface. 

Let a be the radius and p the mean density so that the mass is ^npa^. 
Then with the notation of 3*5 

1 jUnf(R)It’‘dR = i^pa, 

80 that j'f(l{)Ji^dR=iparK 

Differentiating with regard to r, we get 
/(r)r3=:|/)ar, 

and f(r) = ^^ gives the law of density, 
o r 

3*52. Uniform shell bounded by non-intersecting spheres. 

When the boundary spheres are not concentric, let ^4, B be 
their centres and a, 6 their radii. Let p 
be the density of the matter between 
the spheres, and let r, r' be the distances 
of any point P from A and B. The field 
of force due to the given shell is the 
difierence of the fields due to two uni¬ 
form spheres of radii a, b and centres 
A, B. Hence, using the results of 3*31, 
the potential has the following values: 

outside the shell li = |7rp 

in the substance ^ 2 = 

and in the hollow — fTrp {Sa^— + 

Also, using the results of 3*3, the resultant attraction has 

(p 

components as follows: when P is outside the shell 




4-2 
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along PA and ^7rp-j^ along BP; in the substance jirpr along 

PA and i^rp^ along BP; and in the hollow ^Tvpr along PA 

and |7Tpr' along BP. And in this last case the forces compound 
into a resultant force ^rrpBA parallel to BA; so that the force 
inside such a hollow shell is constant in magnitude and 
direction, i.e. a uniform field of force. 

3*6. Long uniform solid cylinder. We may regard a solid 
cylinder as composed of thin cylindrical shells having the same 
axis. We have seen in 2*42 that the attraction of a uniform 
long thin circular cylindrical shell is zero at internal points, 
and at external points is the same as if the mass were con¬ 
centrated uniformly along the axis of the cylinder. For the 
attraction at an internal point P, all the shells external to P 
exert no attraction at P, and the resultant attraction at P 
is therefore due to the shells nearer to the axis than P, so that 
if M be the mass of unit length of the cylinder and a its radius, 
and r the distance of P from the axis, the cylinder attracts 
at P like a rod of mass Mr^ja^ per unit length along the axis; 
i.e. the attraction at P is 

2 {Mr^ld^)lr = 2Mrla^. 

For the potential we have 

dV_ 2Mr 
dr ^ ^ 

so that V = const. — Mr^jd^, 

where r<a. But from 2*42 the potential at an external point 

“ V=2Mlog2l-2Mlogr, 

where 21 is the length of the cylinder. And in order that these 
expressions may give the same value for V on the surface 
r = a (3*7), the internal potential must be 

V = 2Mlog2l-^2Mloga-\- M (1— 

3*7. Variation in the attraction in crossing a surface on which 
there exists a thin layer of attracting matter. Let Pj, Pg be two 
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points on opposite sides of the surface close to one another 
on the same normal to the surface. Let 
m be the mass per unit area of a small 
circular disc of the surface between 
and Pg. 

The points may be regarded as so 
close to the disc that, in accordance 
with the result of 2*5, it exerts equal attractions 27rm at P^ 
and Pg both directed towards the surface. 

The attraction of the rest of the matter on the surface (or 
elsewhere) may be represented in the neighbourhood of P^ 
and Pg by a normal component X taken in the sense P 1 P 2 
and a component Y in some direction tangential to the surface. 
And the values of these components at P^, Pg will differ by 
infinitesimals of the order PiP^- For if X be the value at 

“dX 

Pi and Pi Pg = Sn, the value at Pg is X -f ^ 8n + higher powers, 

and the infinitesimal may be neglected if we neglect the 
thickness of the layer of matter. 

We now compound the attraction of the small disc and the 
attraction of the rest of the matter, and let X^, and 
Xg, Fg denote the components of the whole attraction at 
P^, Pg, taking X^, Xg in the sense P 1 P 2 , and we have 

Xi = X-f27rm, Fi=F 

and Xg = X — 27rw, Fg = F. 

Therefore Xg — X^ = — 47rm, Fg = F^ .(1), 

or the tangential component of attraction is unaltered and 
the normal component is diminished by 47rm in crossing the 
surface. 

As regards the potential, it is clear from 2*51 that a thin 
disc of matter produces the same potential at points on 
opposite sides close to itself, and the potentials of the rest 
of the matter at P^, Pg only differ by an infinitesimal, which 
may be neglected if we neglect the thickness of the layer of 
matter. We conclude that the potential is continuous across 
such a surface as we are considering, in the sense that there 
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is no jump in its numerical value. This is illustrated by the 
case of the thin spherical shell of 3*21. The internal and 
external potentials are represented by different functions, but 
they have the same value on the sphere r = a. 

Reverting now to the result (1), if we use to denote 
the potential functions on opposite sides of the surface, and 
dn to denote an element of the normal taken positively from 
Pi to Pg, we have, at the surface. 


and 




dn dn 


= — 47rm 


( 2 ). 


We remark that m denotes the mass per unit area or the 
‘ surface-density ’ of the layer, and the thickness of the layer 
is neglected; this implies that a finite mass is condensed into 
a volume of no thickness, so that the volume density of the 
matter is infinite. 


3*71 • Attraction o! a thin uniform spherical shell at a 
point of itself. The attraction at a point of itself of a thin 
layer of matter depends on the shape of the gap in the surface 
in which the attracted particle is placed. We may define the 
'principal value of the attraction as the limiting value of the 
attraction at the centre of a circular 
hole when the radius tends to zero. 

We proceed to calculate this for a 
spherical shell. Let m be the mass 
per unit area of the shell, and omitting 
a small circular element of the shell 
surrounding a point P, consider the 
attraction of the rest of the shell at P. 

Let an element QQ^ of area dS subtend 
a solid angle dca at P. The attraction of this element at P is 
mdSjPQ^ along PQ, and resolving this along PO, which is 
clearly the direction of the resultant attraction, we get 
mdScosOPQfPQ^, or mdScos OQPjPQ^, which is equal to 
mdo). If now we allow the gap in the shell round P to shrink 
to vanishing point, we see that we have to take the sum 
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'Lmdw for all cones on one side of the tangent plane at P, so 
that the resultant attraction is 27rm. 


3*8. Mutual work of two attracting systems. Let there 
be a system of particles of masses , mg, mg, ... at the points 

A^, ... and a second system of masses m^', mg', mg',... 
at the points j?g, .... 

If we suppose that initially the particles of the second 
system are scattered at an infinite distance and are brought 
from thence up to their final positions , B^, JBg, ... by the 
attractions of the first system, the work done by these attrac¬ 
tions is represented by 

T^mi' + T 2 ^ 2 ' + ^^ 3 " +••• SFm', 

where ... are the potentials of the first system at 

B^ , B^ , B^ , .... 

Similarly the work done by the attractions of the second 
system in bringing the particles of the first system from a state 
of infinite diffusion up to their final positions A-^, A^, A^, ...is 
represented by 

P^'m^-f l^'mg + Ta'mg-h... or SF'm, 

where Pi', P^', P 3 ', ... are the potentials of the second system 
at A-^^ A^^ A ^, .... 

If we denote the distances A^B^y A^B^, ... A^B^, ... by 
^ 11 . ••• so that 


V^ 


m, m, 


-■f- 


2 


mq 




'11 


'21 


'31 


m, mo m« 
F2=:-^ + -^ + —^4-..., 


^12 


'22 


it is evident that each of the expressions SFm', SF'm is equal 
to the sum Smm'/r extended to the product of masses of 
every pair of particles, one of each system, divided by their 
distance, so that 


SFt^'==SF'm 


(1). 
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For continuous masses we have in like manner 

[Vdm’={v'dm .( 2 ). 


where dm, dm' denote elements of mass of the systems and 
F, F' the potentials of the first system at dm' and of the 
second system at dm. 

These expressions may be called the mutual work of the 
two systems, or the exhaustion of their potential energy. 
For potential energy can be considered to exist in matter in 
a state of infinite diffusion, and work is done and potential 
energy lost when the matter assumes a more condensed form, 
whether under the attractions of another system or under the 
attractions of its component parts. When we speak of the 
mutital potential energy of two attracting systems we should 
therefore prefix a negative sign to either of the expressions 
( 1 ) or ( 2 ). 


3*81. Lost potential energy of a gravitating system. 

Let the system consist of particles of masses m^,m^,m^, ... at 
points .dj,. 42 ,^ 3 ,.... Let us find the work done by the mutual 
attractions of the particles in assuming these positions starting 
from a state of infinite diffusion. Let be in its final position 

A^; the work done by the attraction of the particle m^ on the 
particle m^ as it comes from an infinite distance to .4 2 is 
^ 2 /^ 12 * Then these two particles attract the third and do 
work mimJr^^-{-m 2 mJr 22 in bringing m^ to ^43 and so on; so 
that the total work done may be represented by 
where the summation extends to every pair of particles. 


Let 


Fi = —+ — 


'12 


'13 


V,-. 


m. 


'21 


m, 

'•23 


Then Vi is the potential at .4^ of mg, m 3 ...; is the potential 
at A 2 of m^, m 3 ...; and so on. 
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It follows that the total work done is equal to iSmF, 
the factor | being introduced because the sum SmF contains 
every term such as twice. 

This expression ^SmF represents the work done by the 
mutual attractions of the system of particles or the exhaustion 
of their potential energy as they come together from a state of 
infinite diffusion. 

If the system forms a continuous body, its lost potential 
energy can be found by dividing it into elements dm and 
taking the limit of dm. Here F is the potential at the 
position of dm of the rest of the body, and its Hmit is, by 
definition, the potential of the body at the point to which dm 
contracts. Thus the lost potential energy is 

ijvdm, 

where F is the potential of the body at a point of itself and 
dm is an element of the body at the point. The potential 
energy of the body may be denoted by the same expression 
with a negative sign. 

It used to be supposed that the duration of the sun’s heat might be 
calculated by assuming that, as the matter forming the sun condensed 
from a state of infinite diffusion, the mechanical potential energy lost 
was converted into heat, and then estimating how long it would take 
for this heat to be lost by radiation. 

3*9. Examples, (i) The density of a sphere varies as the depth below 
the surface; shew that the resultant attraction is greatest at a depth equal 
to i of the radiusj and that the valus there is J of the value at the surface. 

Let a be the radius, and let the density p at a distance r from the 
centre be given by p = A (a — r). 

The attraction at a distance E from the centre, when R<a,is due 
to the sphere of radius E and its value is 

= 4^A - (fa - «)»}. 

The greatest value of F occurs where i? = fa, i.e. at a depth equal 
to fa, and this greatest value is |7rAa*; whereas the value at the surface 

= a is seen to be f TrAa*. 
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V^(ii) Prove that, if the law of potential he ~ the potential ofaunifonn 

thin spherical shell of radius a at an external point is the same as that of 

a particle of mass sinli^y ^ that of the shell placed at the centre. 


Let P be an external point at 
a distance c from the centre 0, 
Consider a narrow zone of the 
shell of breadth QQ' at a distance 
r from P. If Q denotes the angle 
QOP, the area of the zone is 
27ra^sin^d^; so that if m is the 
mass of unit area, the potential 
at P due to this zone is 


27rma® - e~^^^ sin Odd, 
r 


and the whole potential is 


[M. T. 1918] 



: 2TTma'^A 


sin 6 do. 


But = o® + — 2ac cos 6, so that rdr = ac sin OdO, and 

V = 27Tm-A — 

c Jc^a c ^ 


= 47rm - A Ae*"®/'' sinh ^ • 
c A 


Also the potential at P of a mass 4Trma® at 0 is, with the same law of 
potential, 4:7Tma^-~e~^i^; hence the result follows. 

The reader will notice that this method of direct integration can be 
applied for finding the potential inside or outside the shell with the 
ordinary law of potential. 

(iii) Prove that the attraction of a uniform elliptic disc at a focus is 
27tw{ 1 — \/(l if the vanishing cavity round the focus is circular 

with the focus as centre, and zero if the cavity is hounded hy a similar 
and similarly situated ellipse with the focus as centre of similarity. 

Let Z/r = 1 — e cos 6 be the equation of the boundary ellipse and r = a 
that of the small circular cavity round the focus. 

The attraction is easily seen to be given by 

f 27 T fl}il—ecoBd) cosO 
I I -—^-mrdOdr, 

Jo Ja r> 

and integrating with regard to r, this 

f27r I 

= m / cos 0 log — 7 ^;- 

Jo ®a(l-ecos^) 
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Then, integrating by parts, this 

r . .. i 


r . l f^^esin^Bde 

m sinolog—7-=-+m / -r-2, 

® a (1 — e cos Jo Jo 1 — e cos d 


— g 1 ’ 


cos ^ H- 


e e(l~eco8 0). 


dd 


But if the inner boundary, which gives the lower limit of integration 
with regard to r, be I'jr = 1 — e cos 6, the integral of ^ with respect 

to r is cos 6 log and integrating with respect to B the result is zero. 

(iv) Prove that the gravitational potential energy of a thin uniform 
circular disc of mass m and radius a is 
C — Sym^j^TTa, where C is independent of o, 
and y is the constant of gravitation. 

[M. T. 1913] 

If a denotes the surface density, the 
potential at any point P on the edge of such 
a disc is easily seen to be 
f\7r f2acoa6 rdSdr 

r 

fin 

= ya 2a 00s BdB = 4yaa. 

-ilT 

Now suppose that the particles which make up the disc are initially 
in a state of diffusion at an infinite distance and that C denotes their 
potential energy in this state. Let the particles be assembled in con¬ 
centric rings. When the radius of the disc is r, the potential at its edge 
is, from above, 4y(Tr, so that the work done in assembling a narrow 
ring of breadth dr is 4yor. 2TTardr. It follows that the work done by 
its own attractions in bringing together the whole disc is 

STryo^J^ r^dr = |7ryo-*a*. 

But w = 7Tao2, SO that the work done =8ym2/37ra, and the potential 
energy remaining = (7 — Sym^lZira. 



EXAMPLES 

1. Taking the earth as a imiform solid sphere of 4000 miles radius, 

shew that a pendulum clock would lose 10*8 sec. per day if taken to 
the bottom of a mine 1 mile in depth. [M. T. 1921] 

2. A small smooth cylindrical hole is drilled through a solid uniform 
sphere. Shew that a particle set free at the surface will oscillate in 
a period independent of the direction of the hole. 
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3* The density of an infinite circular cylinder of radius a at a distance 
r from its axis is proportional to (a^—r*)*. Shew that the attraction 
at an internal point is <> 


where M is the mass of unit length of the cylinder, and y is the gravita¬ 
tional constant. [London Univ. 1926] 

4. A sphere is described on a radius of another sphere as diameter. 

If the latter be a homogeneous sphere of mass M and radius a, shew 
that the resultant attraction on the portion included in the smaller 
sphere is [C. 1914] 

5. A solid sphere of radius a is such that its density at any point is 
proportional to the nth power of the distance of that point from the 
centre of the sphere. Find the potential at any external point, and 
shew that the potential at an internal point at a distance r from the 

vAf /n-}-3 

centre is - -—), where M is the total mass of the sphere. 

n-f2\ a a"+V 

[London Univ. 1931] 

6 . A imiform solid sphere of mass M is cut in two by a diametral 
plane; shew that the resultant attraction between the halves is 

where a is the radius of the sphere, and y the constant of 
gravitation. [C. 1906] 

7. Prove that the force required to separate the two parts of a solid 

uniform sphere divided in any manner is yMM'.OO'ja^t where ikf, M' 
are the masses, and G, O' the centres of gravity of the two parts 
and a is the radius of the sphere. [C. 1906] 

8 . P is a point on a uniform sphere of mass M and radius a. The 
sphere is divided into two parts by a plane perpendicular to the 
diameter through P. Shew that, if the radius of the section by the 
plane subtends an angle a at P, the attraction at P of the portion 
of the sphere which contains P is 3yM cos^ a (1 — f cos 0 L)/a^, 

[London Univ. 1937] 

9. A solid right circular cylinder is of density p and length I, and 
the radius of each end is a. Prove that the potential, due to it, at the 
centre of each end is 

7Typ {I V(a^ -f P) sinh~^ (^/a)}» 

where y is the constant of gravitation. 

If in the above cylinder l^SajS, shew that the potential at the 
centre of the cylinder is 

27rypa^ -f log <5 3). [London Univ. 1935] 

10. Prove that the potential of a uniform circular cylinder, of 
length 2a, radius of base a, and mass M at the middle point of its axis, is 

yM { VI -1 + log, (V 2, +1 )}/o. 

[London Univ. 1938] 
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11. P is a point outside a thin uniform spherical shell at a distance 
c from the centre. Prove that the sphere of centre P and radius c 
cuts the shell into two portions whose potentials at P are equal. 

[London Univ. 1933] 

12. A thin uniform spherical shell of density p has an internal radius 

a and external radius a + a, where (oLja)^ is negligible. Find the attrac¬ 
tion of the shell at a distance a-\-kaL from the centre, where 0<A;< 1. 
Hence shew that, if the outer radius contracts so that the shell becomes 
a thin layer of surface density a and radius a, the attraction of this 
layer at a point in itself is ^Trah, How does this illustrate the discon¬ 
tinuity of normal force in crossing a surface? [C. 1928] 

13. Prove that if a solid sphere of uniform density p is divided into 
two parts by a plane which cuts the surface in a circle of radius 6, the 
attraction between the two parts is iTr^yp^b*. [London Univ. 1936] 

14. A solid homogeneous sphere, of radius R and mass M, is divided 
into two portions by a plane at a perpendicular distance iE from the 
centre. Shew that the mutual attraction between the two portions 

27 

of the sphere is y . [London Univ. 1933] 

15. From an infinite plate of thickness 2a of uniform homogeneous 
matter a spherical portion touching the two plane faces is cut away. 
Prove that the attraction of the plate at a point where the sphere 
touches either of its faces is f of its former value. 

[London Univ. 1927] 

16. A uniform solid sphere of mass M is placed near an infinite 
plate whose surface density is uniform and equal to a. Prove that the 
sphere attracts the plate with a force 2TryMcr. [London Univ. 1932] 

17. Shew that, if M be the mass of a sphere of radius a, the loss of 

gravitational potential energy in assembling the particles from a state 
of diffusion at an infinite distance is fyiUT^/a. [C. 1926] 

18. If the radius of a solid gravitating sphere be a, the density p, 
and an internal point P be at a distance b from the centre O, prove that 
the difference of the potentials at P due to the two portions into which 
the sphere is divided by a plane through P perpendicular to OP is 

(aa _ (oa _ b*)t}. [C. 1908] 

19. The density at distance r from the centre, in a sphere of attract¬ 

ing matter, is po — Ar, where po is the density at the centre, and A is 
constant. Find the potential at any internal point, and prove that the 
potential at the centre is ^ \ 

where a is the radius of the sphere, pi the density at the surface of the 
sphere, and y the constant of gravitation. [London Univ.] 
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20. If AL B is a diameter of a spherical shell of uniform surface density, 
find the position of a plane perpendicular to AB which divides the 
shell into two portions such that their potentials at A are equal. 

If the surface density of the shell is a, find the attraction at its 
centre due to each portion of the shell. [London Univ. 1938] 

21. Shew that the initial rate of decrease of in descending a mine 
shaft would be equal to gja if the density (p) of the earth (radius a) 
were uniform. But if the earth had a spherical nucleus of different 
density and of radius 6, the density of this nucleus must be 

r l~Aa«) 

where Xgja is the initial rate of decrease of g in descending the shaft. 

[C. 1930] 

22. A sphere of radius a, mass Af, and density varying directly as 
the distance from the centre, is built up of matter brought in small 
amounts from a state of complete dispersion at infinity. Shew that the 
work W made available during the process is given by W 

where y is the constant of gravitation. 

Shew that, if the matter is now redistributed so as to form a sphere 
of the same radius but of uniform density, a further amount of work 
is made available. [P. 1936] 

23. Find the total work done by the gravitational forces in collecting 
from infinite dispersion six equal particles of mass m and placing them 
at the corners of a regular hexagon of side 2a. [London Univ, 1932] 

24. Six particles, each of mass m, are placed symmetrically on the 

circumference of a circle of radius r and are allowed to move from rest 
under their mutual attractions. Find the velocity of each particle 
when they are on a circle of radius Jr. [London Univ. 1933] 


26. Assuming the earth to be a homogeneous sphere of density a 
and radius a, surrounded by a spherical shell of density p and thickness 
6— a, shew that the weight of a body given by a spring balance is 
unaltered at first on descending vertically from the surface if 
p = 2o*or/(2a® + 6®). The rotation of the earth may be neglected. 

[M. T. 1914] 

26. A solid sphere of radius a and uniform density is surroimded 
by a concentric shell of internal radius a and external radius 6 in which 
the density is a function of the distance from the centre. Find the law 
of density if the force of attraction is constant in the substance of the 
shell, and prove that the mass of the shell is to the mass of the sphere 
as 6®—a® is to a®. [C. 1932] 
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27. The density of a sphere of radius a and mean density /> is a 
function of the distance from the centre such that the attraction at 
any point inside the sphere is proportional to the square of the distance 
from the centre. Express the density at a distance r from the centre 
in terms of p, a and r, and prove that the potential at an interna 

point at a distance r from the centre is ^ ^ (4a® — r®). [P. 1933] 

y a 


28. If 0 be the centre of a solid uniform sphere and P an external 
point, shew that the sphere on OP as diameter cuts off a portion whose 


attraction at P is equal to 


1 M 

2 OP®’ 


where M is the whole mass. 

[C. 1905] 


< 29. Find the potential at any point in the interior of a sphere of 
r£tdius a, the density of the sphere at a distance r from the centre being 
p„(l-fcr"). 

Shew that, if the sphere be divided into two parts by a plane through 
the centre, the mutual pressure between the parts is 




[London Univ. 1934] 


30. A portion of a uniform solid paraboloid of revolution is cut 

off by a plane at right angles to the axis OX. Prove that the attraction 
at the focus S will be zero when sin(JPiS'X)= 1/V^» ^ being a point 
on the rim, [M. T. 1905] 

31. If p be the density of a homogeneous solid sphere and a the 

radius, prove that the force in astronomical units exerted by the rest 
of the sphere on a conical portion of itself bounded by a right circular 
cone of angle 2 a with its vertex at the centre is sin® a. 

[C. 1915] 


32. A solid of imiform density p is in the form of the spheroid 
obtained by rotating an ellipse, of latus rectum 21 and eccentricity e, 
about the major axis. Prove that the potential and the intensity of 
the attraction at a focus of the generating ellipse are 

27rypZ®/( 1 - e®) and (27rypZ/e®) [log {(1 + e)/( 1 - e)} -- 2 e] 
respectively. [London Univ. 1936] 

33. Find the potential and attraction, at any point, of a homo¬ 
geneous sphere with a spherical excentric cavity, and deduce the 
corresponding values for a thin spherical shell whose small thickness 
at any point is proportional to its distance from a plane. [C. 1899] 

34. Two uniform solid spheres are of density p and radius a, and 
the distance between their centres is c (> 2a). Prove that they attract 
each other with a force 167 r®yp%®/ 9 c®, 

where y is the constant of gravitation. 


[London Univ. 1926] 
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36. A solid sphere of radius a and mass M is of uniform density. 
Prove that the attraction on an octant of the sphere due to the rest 
of the sphere has components parallel to the edges of the octant each 
equal to where y is the gravitational constant. [C. 1914] 


36. Prove that the work done by the gravitational forces when 
a uniform spherical shell, of mass M and of external and internal 
radii a and b respectively, is changed to a solid sphere of the same 


density is 


where c^ = a^-‘ 6®. 


3yM^ 

10 c« 


[2c5~2a‘'-365-f 5a263], 


[London Univ. 1939] 


37. The intensity of gravitational attraction at any point on the 
surface of a certain uniform solid sphere, of radius a and centre 0, is gr. 
A spherical cavity of radius 6 is made with its centre at a point C, 
where OG — c, (6 + c<a). Prove that the magnitude's of the intensities 
at the points where the line 00 (produced) meets the surface of the 
sphere are , ^8 ] 

and that at the ends of any diameter through O perpendicular to OC 

they are ^ 26^ 6® 1 i 

(a2 + c2)8’^a*(a2 + c^)2J 

[London Univ. 1935] 


38. A uniform spherical shell of gravitating matter, in which the 

volume of the matter is equal to that of the hollow it encloses, contracts 
under its own gravitation into a uniform solid sphere of the same 
density. Prove that the work done by the gravitating forces is 
3 (2i— l)yM2/10a, where M is the mass and a the radius of the solid 
sphere. [London Univ.] 

39. The density of a sphere, of mass M and radius a, at a distance r 

from the centre is proportional to a — r. Prove that the potential at 
this distance is yM(2a* — 2ar2-|-r®)/a^, and that the potential energy 
exhausted in collecting the elements of the sphere is 26yM®/36a, where 
y is the gravitational constant. [London Univ. 1926] 

40. Find the negative potential energy W for a spherical distribution 

of mass M and radius a in which the density varies uniformly from pi 
at the surface to po centre, and deduce that, if pi = 0, 

W=:26GM^I36a; 

while if po = 0, pi=#0 this value of W is reduced in the ratio 10/13. 
G denotes the gravitation constant. [M. T. 1928] 

41. Prove that, if a segment of a uniform thin spherical shell has 

an angular radius a, its attractions on the two sides, at points close to 
the pole of the segment, are 2^^y<7(l ±sin^a), where a is the surface 
density. [M. T. 1922] 
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42. Defining the principal value of the attrs/ction at a point of an 
attracting surface as the limit of the attraction at the centre of a 
circular hole of which the radius tends to zero, shew that its tangential 
components a<re continuous with those of the attraction near the 
surface, while its normal component is the mean of the limiting normal 
components on either side. 

Shew that the attraction at the focus of a small elliptical hole of 
eccentricity e in a thin uniform spherical shell makes an angle 

with the radius of the shell. [C. 1929] 

43. Prove that the work done by a gravitating system in attracting 

a new body from infinity is equal to iJp(V + V') dr, where the integral 

is taken through the space occupied by the new body, and F, F' are 
the gravitational potentials before and after the new body arrives. 

Calculate the work done in filling up a spherical hollow of radius b 
in a uniform sphere of radius a, the distance between the centres being 
c and the density p. [C. 1914] 

44. Find an expression for the mutual gravitational potential energy 
of a uniform sphere of mass M and a uniform rod of length 2a and mass 
m when they are placed in any relative position. 

Shew that this is constant when the centre of the sphere moves on 
the spheroid formed by rotating an ellipse whose foci are the ends of 
the rod about the rod. [M. T. 1923] 

46. A uniform rod BO of mass M is free to turn about its middle 
point O which is fixed, and a uniform sphere of mass m has its centre 
at A, Shew that the couple on the rod tending to turn it into the 

line OA is \yMm > where 6, c denote the lengths of 

AC, AB and a the angle BAG, [M. T. 1930] 

46. The profile of a long straight moimtain range of uniform 
density p is an arc of a circle less than a semicircle. Prove that if h 
is the height of the range and h its breewlth at its base, a plumb line at 
its foot will be deflected towards the range through an angle approxi- 

mately equal to 3o6 , 2A 

.n. tan~^ , 

27rpQB b 

where B and p^ are the radius and mean density of the earth. 

[London Univ. 1927] 

47. Shew that, if the potential of a particle of mass m at distance r 
is rriift"(r)lr, then the potential of a uniform solid sphere of mass M 
and radius a at an external point at distance r from the centre is 

g.{a[f (r+a) + ^'(r-a)]-[^(r+o)-^(r-a)]}. 


RNA 


5 
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where the functions (x) are successive derivatives 

of <l> (x). 

Consider in particular (a:) = where A^O. Shew that the 
external field of the sphere is the same as that of a single particle at 
its centre; and that the mass of this particle is 

3M (Aa cosh Aa — sinh Xa)IX^a^ 

if A>0, orAfif A = 0. [M. T. 1924] 

48. Shew that the potential of a uniform circular disc at a point of 
itself is 

r27r aa(a — r cos 6) dd 
Jo V(^^ — cos^ + 7-2)* 

where r is the distance of the point from the centre and cr is the mass 
per unit area. 

Deduce that the work done in removing all the particles to infinity is 
iTTffV. [C. 1924] 

49. If the attraction between two elements of matter of mass m 
and m' is mm' f(R)y where R is the distance between them, shew that 
the inward acceleration due to a uniform thin spherical shell of mass M 
and internal radius a, at an internal point distant x{<a) from the 
centre, is F, where 

If F is zero for all values of a and x, subject to the condition xKa, 
shew that the only possible form of f{R) is AjR^^ where A is some 
constant. [M. T. 1937] 



Chapter IV 


THEOREMS OF LAPLACE, POISSON 
AND GAUSS. GENERAL THEORY 


4*1. Laplace’s equation for the potential, (i) Let V be 

the potential of a system of attracting particles at a point 
{x,y,z) not in contact with the particles, so that V^Itinlr, 
where m is the mass of the particle at (a, 6, c) and 

= (x — (y — b)^ (z — c)^. 


Then 


dV __ mdr ^ m{x — a) 

dx dx 


and 


dx^ r® 


Similarly 

and 


dy^ r3 ’ 

dW m m(z-c)^ 
dz^ + 


Whence, by addition, 

dW dW dW _ 

(ii) Let V be the potential of a continuous body or bodies 

at a point (x, y, z) outside the body or bodies, so that V = 

where p is the density of the 
element of volume dv at 
{x\ y\ z'), and 

+ (2/-2/')^+ (2^-2:')^. 

Since r is measured from an 
external point P there is noth¬ 
ing to make the integral other 
than convergent, and the re¬ 
sults of dififerentiation with regard to x, y and z under the 




5-2 
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integral sign are continuous and yield convergent integrals, 
therefore we may dififerentiate under the integral sign and 
write 


and 

Similarly 

and 


dx J 

=-Jt 


p{x — x') 


dv 


dx^ 

022 = 


(p_ 

^p{x — x'Y\ 


J 


^pi.y-y'f\ 


r® j 


1 

1 


r® 1 


'^dv. 


dv 


dv. 


Whence, as before, by addition, 

dw dw 

0X2^91,2 

or V2F = 0) 


.( 1 ). 


This result is known as Laplace’s equation.* It is satisfied 
by the potential of an attracting system at every point at 
which there is no matter. 

Poisson’s equation for the potential. Now let the 

point P of co-ordinates x, y, z be inside the attracting matter. 
Describe a sphere of small radius e and centre (a, 6, c) con¬ 
taining the point P, taking € so small that we may regard 
the density p of the matter in this sphere as uniformly distri¬ 
buted. The matter which produces the potential F at P may 
now be divided into two parts, viz. the matter outside and 
the matter inside the small sphere. Let denote their 

contributions to the whole potential F at P. Since P is not 
in contact with the matter which produces the potential Tj, 
* therefore, by 4*1, = And being the potential at 


* Pierre Simon, Marquis de Laplace (1749-1827), distinguished 
French mathematician. His great work Micanique Celeste in five 
volumes was published between 1799 and 1826. The theorem above is 
to be found in Micanique Cihste, t. i, p, 137, Paris, An VII. 
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a point (x^y,z) inside a small sphere of radius e, we have, 
from 3-31, 

where r is the distance of (x, y, z) from the centre (a, 6, c), 

= §7rp -(x- a)2 - (y ~ 6)2 - (2 - c)^}, 

dV 02F 

Hence = - ^np (x - a), and -|7rp. 


Similarly 


dy^ dz^ 




a^F, a^F, a*F„ 

ax^'^ay*"''a22“ 


But F = 1^ + 14; therefore at every point at which there is 
attracting matter of volume density p 

S2y S2y 327 ^ 

dx^dy^dz^ .( 1 ), 


This result is known as Poisson’s equation.* 

^^•12. Attraction = grad F, inside the body. We may 

adopt the line of argument of 4*11 to shew that the relations 

X, r, Z = ^^, which we have seen to be true at points 

outside a system of attracting particles, are also true at all 
points outside or inside a continuous distribution of matter. 
The definitions of F and X for a continuous distribution are 


and X 


J 


* Sim6on Denis Poisson (1781-1840). French mathematician and 
physicist, author of more than 300 papers. The proof given in 4*11 is 
very nearly in the form given by Poisson in Nouveau Bulletin des 
Sciences par la Societi Philomaihique de Paris, t. iii, p. 388; 5® Ann6e, 
1812. He observes that, in the case of a body of continuously varying 
density, since V®Fi = 0, it follows that V^F and also V^Fj must be 
independent of the form and dimensions of the element of which F* 
is the potential, so that this element may be taken to be so small that 
its variation in density may be neglected. An alternative proof is given 
in 4*21 and 4*2. 
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where the symbols are as defined in 4-1, and where for an 
internal point (*, y, z) the definitions need the amplification 
of 3-1. 

When the point (*, y, z) is outside the matter, differentiation 
imder the integral sign is permissible, and since 


r^=(x—x')^ +{y—y'Y + {z-z'f. 


dx 




But when the point {x, y, z) is inside the matter we proceed as 
in 4*11 and use suffixes 1 and 2 to denote contributions arising 
from the matter outside and the matter inside the small 
sphere. Since the point {x^y^z) is not in contact with the 
matter outside the small sphere, we have as above Xj = dVJdx, 
and for the matter inside the sphere 


and 


V, = §7rp (362 ^i^x-af^{y-bf-{z^ c)^} 

xr A X — a A / X 

^ 2 = -^pr.-y-= -^p(x-a) = ^. 


Hence, for the whole attraction component, 




dx dx dx ’ 


at points inside the attracting matter. 


4*13. To summarize our results up to this point: 

(i) The attraction is the gradient of a potential function V 
both outside and inside the attracting matter (2*22, 4*12). 

(ii) In empty space the potential satisfies Laplace’s equation 
V2F = 0 (4-1). 

(iii) At any point at which there is matter of volume 
density p the potential satisfies Poisson’s equation VW = — 47rp 
(4-11). 

(iv) When there is a surface distribution of matter the 

potential function assumes different forms Pi, on opposite 
sides of the surface, but at the surface they satisfy the con¬ 
ditions T7_T7 

ri—vz 
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dV, 0Fi 


= — 47r?n, 


where m is the surface density of the matter, and dn is an 
element of the normal directed from region 1 to region 2 (3’7). 

V^4-2. Gauss's theorem.* The outward flux of the force of 
attraction over any closed surface in a gravitational field offeree 
is equal to —Afn times the mass enclosed by the surface. 

This theorem can be deduced in two steps, from Green's 
theorem 1'3 and Poisson's equation 4*11. 

Thus, in the notation of 1*3, the outward normal com¬ 
ponent of force across an element dS of the closed surface is 
IX + mY + nZ and, using dv for an element of the volume 
enclosed, 


j(ix + mY+nZ)d8=: + 


dY 

dy^dz) 


^ 

[dx^ ^ 03 *) 


= —4TTjpdv (Poisson’s equation) 

= — 47r (mass enclosed). 

Conversely, Poisson's equation may be deduced from 
Gauss's theorem. For if we write jpdv for the mass enclosed 
by Sy we have 

— 4t7rjpdv=j (IX + w F+ nZ) dS (Gauss) 

^ r/dx dY dz\^ 


for the mass enclosed 


(0a: 0y 02 j 


(0a:*^0y*'^02*, 


Therefore 


f0*F 0*F 0*F , 




* Carl Friedrich Gauss (1777-1855). German mathematician and 
physicist. 
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And this is true for integration through the volume enclosed 
by every closed surface, i.e. for all ranges of integration, so 
that the integrand must vanish, or V^F = — 47r/>.* 

4-21. Independent proof of Gauss's theorem. Regarding the 
field as due to a system of particles, let m be the mass of a particle situated 
at 0. Let a cone of small solid angle dm and vertex O cut a closed 



surface 8 in elements of area dS^^ dS 2 y dS^f ••• at P^, P^, Ps, 
and let the outward drawn normals at Pj, Pg, Ps, ••• make angles 
^ 2 » ^8> ••• with the line ... PgPgPiO. Then so far eis the outward 
flux of force from 8 depends on the particle m at 0, the contributions 

of the elements d 8 i , d/S' 2 , d 82 , ... are each of the form cos 6 . But 


♦ In vector symbolism, we have the resultant attraction R = grad V , 
and the flux of R out of a region bounded by a closed surface 8 

nd8 


=jdivHdv (1-7) 
=JdivgradVdv 




V>Fdv=-47r 


J pdv 


and conversely, if 

— 47r Jpdv=j Ef^dS 


= —477 (mass enclosed); 


=/■ 

=/ 


divRdt;= j div grad 
V*Fdt;, 


Vdv 


then J( V* F + ^vp) dv = 0 , for all ranges of integration, and. V*F = — irrp. 
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this expression is equal to ± mcfoj; + or — according as 0 is acute or 
obtuse, i.e. + or — according as the cone proceeding fix)m O is entering 
or leaving the region bounded by S. It follows that when 0 is outside 
the region (as in the figure), since the cone leaves the region as often 
as it enters it, therefore the total contribution to the outward normal 
flux arising from the intersections of this cone and the surface is zero. 
But when O lies inside S the cone leaves the region and if it re-enters, 
it finally emerges, so that it leaves the region once more often 
than it enters it, and the total contribution in this case is — 
Hence by taking cones in all directions round O, we see that a particle 
of mass m contributes zero to the outward normal flux when it lies 
outside the surface, and — 47rm when it lies within the surface. And 
by summing for all particles of the system we obtain the required 
result, viz. total outward normal flux= ~ 47r (mass enclosed). 

4*22. When the field is due to the attraction of a continuous body 
and the surface S intersects the body, 
we may regard the body as divided into 
distinct parts of masses Mi and M ^,, inside 
and outside S, and suppose that these are 
separated by an infinitesimal distance so 
that Mi lies wholly within S and Mq 
wholly without it. The outward normal 
flux of force from S is then —^nMi as 
before, and this will continue to bo true 
as the distance between Mi and M^ 
decreases without limit. 

The result is not affected by the ‘ matter 
on the surface because we are con¬ 
sidering a body of finite volmne density 
and, as explained in 3*7, it requires a finite surface density of matter 
(i.e. an infinite volume density) to affect the normal component of 
force across the surface. 

There is, however, an exceptional case in which the theorem needs 
modification, namely when the matter is a distribution of finite surface 
density on the surface S, Reverting to 4*21 we observe that if the 
point O lies on the surface S, then only cones on one side of the tangent 
plane at S intersect the surface, so that the contribution of the mass m 
to the total outward flux is — 27rm, and the same being true for all 
matter lying on the surface, the total outward flux due to a total mass 
ilf is — 27rM, 

It is easy to verify, for a uniform distribution of mass M on the 
surface of a sphere /S', that the total outward flux of force 

across a concentric sphere inside S is zero (3*2), 

across /S is — 2itM (3*71), 

and across a concentric sphere outside /S is — ^ttM (3*2). 
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^4*23. Applications of Gauss’s Theorem to spheres and 
cylinders. 

(a) Spherical shell. Let M be the mass of a uniform spherical shell 
of radius a. Apply Gatiss’s theorem to a concentric sphere S of radius r. 
By symmetry the outward normal force R has the same value at all 
points of S, 

When r<o, fig. (i), jS encloses no matter, so that 

4tTrr'^R = 0, and i? = 0 at all points inside the shell. 

When r>a, fig. (ii), S encloses a mass M, so that 
AiTtr^R = — 47TAf, and i? = — Af/r* at all points outside the shell (3‘2). 




00 


(6) Solid sphere of uniform density p. To find the attraction at an 
internal point, apply Gauss’s theorem to a concentric sphere S of 
radius r < a, fig. (i). S now encloses a mass so that 

4TTr^R==: — An.^vpr^ and —f^v-pr (3*3). 

The attraction at an external point may be found as for a uniform 
spherical shell. 

(c) Long uniform cylinder. Let the figures now represent cross- 
sections of a iiniform long cylinder of radius a and mass M per imit 
length. At points sufficiently distant from the ends of the cylinder 
we may neglect force parallel to the axis, and as above let R denote 
radial force. Then we may take for S a imit length of a coaxial cylinder 
of radius r and ignore any force across its plane ends. 

When r<a, fig. (i), /S' encloses no matter, so that 

2irrR = 0, and i? = 0 at all points inside the cylinder. 

When r > a, fig. (ii), S encloses a mass M, so that 
2rrTR = — 47rAf, and jR = ~ %Mjr at all points outside the cylinder (2«42). 
/ 

^4*24. Tubes of force. If the lines of force (2*26) be drawn 
through every point of a small closed curve drawn on an 
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equipotential surface, these lines form a tube of which the 
equipotential surface is a cross-section, called a tube of 
force. 

Apply Gauss’s theorem to the region bounded by a portion 


a tube of force cut off be¬ 
tween two equipotential sur¬ 
faces and not enclosing any 
matter. If denote the 

force of the field at the ends of 
the tube and co ^, ojg the areas 
of the ends, since there is no 
flux of force through the sides 
of the tube and no matter 
inside it, therefore 



2a>2 = 0, 


or the force varies inversely as the area of the cross-section 
of the tube. 


Example. If the field of force be due to the attraction of a uniform 
infinite plane sheet of matter of mass m 
per unit area, it is evident that the equi¬ 
potential surfaces are planes parallel to 
the given plane and the tubes of force are 
straight, so that the force is constant. If 
B denotes the force directed away from ^ A 7\ H—> 

the plane and we apply Gauss’s theorem R ^^^ ^ 

to a portion of a tube of force of cross- 
section w extending to both sides of the 
plane, since by symmetry the force has 
the same value on both sides of the 
plane 

2Rcjd = — 4^770)171, or i? = ~ 27rm (2*6). 


V^3 


-^•3. Theorems on the potential, (i) The potential cannot 
have a maximum or a minimum value at any point of space 
unoccupied by matter. 

For if the potential had a maximum value at a point 0 
we could surround O by a small sphere 8 such that at every 
point of 8 the potential would be less than at 0, The outward 
normal force would therefore be negative at aU points of 
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S, and js^dS taken over S could not vanish, as it must if S 


contains no attracting matter. Similarly for a minimum. 

As a consequence of this theorem we note that at a point 
in free space, if the potential be not constant in the neighbour¬ 
hood of the point it increases in some directions and decreases 
in others, so that, although a particle placed at the point 
free from constraint might be in equilibrium its equilibrium 
could neither be stable, nor unstable for all displacements. 

This latter theorem is due to Earnshaw* 

(ii) If the potential be constant over a closed surface contain¬ 
ing none of the attracting matter, it must be constant throughout 
the interior. For otherwise it would have a maximum or 
minimum value in the interior, which is impossible by (i). 

(iii) Gauss's Mean Value Theorem. The mean value over 
a spherical surface of the potential of any attracting system is 
equal to the potential at the centre of the sphere due to the mass 
outside the sphere, plus the mass inside the sphere divided by 
the radius. 

Let a be the radius of the sphere, m the mass of a particle 
at an external point A and m' 
the mass of a particle at an \9 

internal point A'. 

The potential at any point Q A 

on the sphere is 

^ m m' 

The mean value of F over the sphere is -^^jvdS, where 
dS denotes an element of area at Q. 



This 


m r dS m' C dS 


r dS 

But J is the potential at A due to a spherical layer of 
unit surface density and therefore = 4^a^lOA . 
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r ds 

Similarly, J is the potential at due to the same layer 
and therefore = ^rra^ja = 47ra. 

Ttt 7Yl/ 

Hence the mean value of F = 7^7 ^- 1 

OA a 

In like manner the mean value of the potential due to any 
attracting system = S -f , where the first sum repre- 

iyJL CL 


sents the potential at 0 due to the matter outside the sphere, 
and thus the theorem is proved. 


(iv) If the potential be constant throughout any region T, 
unoccupied by matter, it has the same constant value through 
all space that can be reached from T without passing through 
matter (Gauss). 

Let V be the value of the potential throughout T, If the 
potential differs from V in neighbouring space there must be 
some regions contiguous to T where the potential is greater 
than F and others where it is less. In the neighbourhood 
of a place where the potential is greater than F draw a sphere 
with its centre 0 inside T, Then on the part of the sphere 
inside T the potential is F and on the part outside it is greater 
than F, so that the mean value of the potential over the 
sphere exceeds F, the value at the centre O of the sphere, 
which is impossible by the last theorem since the sphere con¬ 
tains no matter. Hence the region T of constant potential 
must be extended to include the whole of the sphere, and 
so on step by step through all space that does not include 
attracting matter. 

(v) In any case of symmetry about an axis, if the potential 
be constant through any finite distance along the axis, however 
short, not within the attracting matter, it has the same constant 
value through all space that can be reached from this portion 
of the axis without passing through matter* Without loss of 
generality we may suppose that the potential is zero at all 
points of a segment of the axis. Let P be a point of the 


♦ Thomson and Tait, Treatise on Natural Philosophy, Art 4d8, Oxford 
1867. The proof given above is due to Dr S. Verblunsky. 
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segment. On account of the axial symmetry, there is no force 
in any direction through P perpendicular to the axis. Hence 
the axis is normal at P to a surface on which V is constant; 
i.e, F = Tp = 0. Since the argument applies to any point P of 
the segment, it follows that F = 0 throughout a region. The 
theorem therefore follows from (iv). 

V 4'31. Comparison theorems, (i) If two different distribu- 
tioTis of matter have equal potentials over any closed surface not 
including any attracting matter, they have equal potentials 
throughout the space enclosed and throrigh all external space 
which can be reached without traversing any of the attracting 
matter. 

Let the attraction of one of the distributions be changed 
to repulsion. This changes the sign of its potential and makes 
the potential due to the two distributions taken together zero 
over the given surface. Therefore by 4*3 (ii) the potential is 
zero throughout the space enclosed, and by 4*3 (iv) zero 
throughout all external space which can be reached without 
passing through matter. Then changing back from repulsion 
to attraction, it follows that the potentials of the two distribu¬ 
tions are equal. 

(ii) If two different distributions of matter produce equal 
potentials over any surface enclosing both, they produce equal 
potentials throughout all space external to the surface. 

Let the attraction of one of the distributions be changed 
to repulsion. The two distributions combined then produce 
zero potential over the given surface, and they also produce 
zero potential over the infinite sphere. Therefore their com¬ 
bined potential is zero throughout the space between the given 
surface and the infinite sphere, for this space contains none 
of the matter (4*3 (ii)). Reverting to the original distributions, 
it follows that they have equal potentials throughout the 
same region. 

(iii) If two finite distributions of matter have the same equi~ 
potential surfaces outside both distributions, their attractions at 
every external point have the same direction and are proportional 
to the masses of the distributions. 
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The attractions at any external point are in the same 
direction because they are along the inward normal to the 
equipotential surface through the point. 

Let V, V' denote the potentials of the two distributions. 
Since V' is constant whenever V is constant, there must be 
a functional relation between them; say F' ==/(F). Therefore 


dx 


=f'(V) 


dx 


and 


dx^ 




dx ^' 


Adding the similar results obtained by differentiating with 
regard to y and z, we get 

But, outside the matter, 

V2F = V2F' = 0. 




Therefore / 


Hence, unless V is constant everywhere to make the second 

factor vanish, we have ^ 

J ( K j — U, 

giving, on integration, 

F'=/(F) = AF + JS, 

where A, B are constants. But V and F' vanish at infinity 
in the ratio of the masses M, M' of the distributions, so that 

M' 

J5 = 0 and A = M' jM, Hence F' = F, and the result follows. 


4*32. Condition that a family of surfaces is a possible 
family of equipotential surfaces in free space. A family 
of surfaces in general is not a family of equipotential surfaces 
unless its equation satisfies a certain condition. To find the 
condition that the equation 

fix, y,z) = const. 

may represent a family of equipotential surfaces. 
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If the potential V is constant whenever f(z, y, z) is constant, 
there must be a functional relation between F and/ (a:, y, 2 ); say 

V=<f>{S{x,y,z)). 




By adding similar results in y and z, we get 


Then 


and 



But in free space V^F = 0; therefore 

Vif) 


dH dH dH 


\aa;/ '^[dy) ■^\02/ 




= a function of /, say x if) 


.( 1 ). 


This is the necessary condition, and when it is satisfied the 
potential F can be expressed in terms of f{x;y,z) thus: 

We have F = ^ (/), where 

^k’U) 


Therefore 


or 


^'(f)+X(f) 0 - 

log<f>'{f)^A-jx{S)df 

4 ,' 

It follows that V=<f>(f) = AL-ix^^fdf+B .(2). 

/ 

n/ 4*33. Bxample. Shew that a family of right circular cones with a 
common axis and vertex is a possible family of equipotential surfaces^ 
and find the potential function. 

Taking the axis of z for the common axis, the equation of the family 

of cones is x* + v* 

/(!r,y, 2 ) = --^ = const.(1), 
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The condition 4'32 (1) in this case becomes 

2,2, 6(x» + y*) 

_ 2+3/ 

■^'(/) 4x\4y\ 4 (!»»+ !,»)« 2/(1+/)’ 

-4 ,4 + -4 


and the condition is satisfied since the expression on the right has 
reduced to a function of /. Hence 

^'(/) V^2(/+l) 

and, by integration, 

log <t>' (/) + log/+ i log (/+!) = const. 

Put/= tan^ B, and we get dcf) = 


Hence F = <^ (/) = 20log tan \d-{- O'. 

It is clear that B is the half-angle of a cone, and F is constant when B 
is constant. 

The truth of this proposition is evident at once if we consider 
Lapletce^s equation in polar co-ordinates, namely 


i^/^aFx 1 a/. aF\ i dw 

dr \ h‘)'^ r^ sin B dB V dB sin* B d<l>^ 


= 0 , 


and whether it can have a solution which makes F constant over the 
surface of a right circular cone. There is clearly such a solution if F 
is independent of r and <j6 and satisfies 


• A 

sme^=^. 


where ^ is a constant, and this leads to 

F =:^logtan j^B+B, 

as before. 

^4-4. When the potential is given, to find the distribution 
of matter. We have seen that potential functions may become 
infinite at points or lines where matter is concentrated, and 
that otherwise potential functions are finite and continuous, 
even when crossing surfaces on which matter is concentrated 
(3*7); but in crossing such a surface there is a discontinuity 
in the gradient of the potential. Thus a closed surface on 
which matter is concentrated divides space into regions in 


RN A 


6 
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which the potential function takes different forms, though 
these forms take equal values on the surface itself. 

When the potential is given at all points of space, Poisson’s 

equation, in the form p— serves to determine the 

4[7T 

volume density of matter wherever a finite volume density 
exists. If the potential is given by different functions , Pg 
on opposite sides of a surface S, this implies that there is 
on 8 matter of surface density a determined by the relation 


477 dn ) 


(3*7), 


where the direction of the normal dn is from 1 to 2. 

If the potential becomes infinite at a point, the amount of 
matter concentrated there can be found by applying Gauss’s 
theorem to a small sphere having its centre at the point. 

If the potential becomes infinite at all points of a line, 
the line density of the matter can be found by applying 
Gauss’s theorem to a short cylinder whose axis is the line and 
whose radius tends to zero. 

y 

4-41 • Examples, (i) The potential outside a certain cylindrical 
boundary is zero; inside it is V^x^-^xy^ — ax^-^Say^. Find the 
distribution of matter. 

We have first to find the boundary. Since the potential is continuous 
across the boundary and zero outside, the boundary must be given 
by 



and outside /> = 0 since Fj = 0. 
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At P on AB 


or 


47r \dx dx J a.„o 
= j (o» - 33/») = ~ {MA<‘ - MP*), 

4:7r 47r 

a=%AP.PB. 


At Q on AO (a;= V3^) 

^Ya= 1_ f _ 1 

47r ^ 47r ( 2 2 dy Ja,„viy 

= + f 2/^ + aa; - 3 V 3x?/ + 3 V 3a2/}«.v'iy 


or 



(a-x), 



OQ,QA 


and similarly on OB. This shews that a solid prism of uniform density 
ajiT would produce the same external field as a distribution of matter 
of ^rface density ZAP.PBj^Tr on each of the faces of the prism. 


Vvii) Use the theorems of Laplace and Gauss to verify that the function. 


V = m log 


r-Hr'-f2a . 


is the potential due to a uniform line density m on' 


^ r -h r' ~ 2a 

the straight line joining the points A, B from which r, r' are measured. 
Take the origin at the middle point of A-B and the axis Ox along BA, 
and let ^ denote distance from Ox. 



We notice that F is a finite and continuous function except on the 
line between B and A where it becomes infinite. 


We have r®=(ir —= {a?+«)* +S*, 


and 



r + + 2 a — r + + 2a \ 

r-hr' —2a* —r-fr'-f 2a/ 


=:mlog 


~ + r'^ + 4a^ H- 4 gr^ _ . r' -f a: 4- a 
— r* 4* r'* — 4a*-f 4^ ~ ^ r -k-x — a ’ 


6-2 
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Hence 


dV_ m + ® . |\ ^ 1 

dx^r' + x + a\ r' ) r-\-x — a\ r 


and 

dW m(xAa) m(x — a) 

dx» r'» ' r* . 

.(1). 

Similarly, 

dV _ m { m f 

-hx-ha^r' r-j-x — a‘r 



_m /r' — x — a r~-x + a\ 

“ { \ r' r ) 


or 

ydV f x+a x—a\ 

) . 

.(2); 

and 

1 d m(x-ha) m(x — a) 

r'o r» . 

.(3). 


By adding (1) and (3) we get 


1^/ aF\ 


0 , 


shewing that, since F is symmetrical about Ox, it satisfies Laplace’s 
equation in cylindrical co-ordinates x. This holds good at every 
point save at points between A and B on the line AB where F is 
infinite, so there is no volume density of matter; and since there is no 
discontinuity in the form of F there is no surface density. Hence the 
matter is confined to the line AB, It remains to shew that the line 
density is uniform. 


B 



A 


Apply Gauss’s theorem to a small cylinder of length c and radius { 
with its centre at any point P onAB and its axis along AB, The flux 
across the curved surface is 


27r€{^=27r€m(^-^^4-^y^^, from (2), 

-►--47r€m as {->0; 

02 F 

and the flux across the plane ends is , which tends to zero 

with I, It follows from Gauss’s theorem that the mass on any element 
of AB of length c is so that the line density is riniform and equal 
to i 


to m. 


^4*5. Simple applications of Laplace’s equation. 

(i) Field due to a uniform plane sheet. Taking the axis of z at right 
angles to the sheet, F is clearly a function of z only, so that Laplace’s 
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equation reduces to dWjdz^^O, which gives on integration dVjdz^A, 
where A is the constant value of the force in the direction in which z 
increases; and, by 2*5 or 4*24, this is — 27rm, where m is the mass per 
unit area of the sheet, so that 


and F = — 27Tmz -h const. (2*51). 

(ii) Symmetry about a point. When the field is symmetrical in 
such a way that F is a function of r only, in three dimensions, we can use 
the transformation of Laplace’s equation into 
polar co-ordinates given in 1 *68, putting dV jdB ~ 0 
and dVld<f> = 0, or we can obtain the required /y^ 
form directly by applying Gauss’s theorem to f/ x/ \\ 

the space between concentric spheres of radii r j j \ \ 

and r 4- dr. Assuming that this space contains 11 II 

no matter and that the force at distance r is \\ // 

dV/dr, the flux of force out of the region across 

dV 

the surface of the inner sphere is — 477r® ^ -. 

The outward flux across the sphere of radius r dr is therefore 

to the first power of dr, and the vanishing of the total outward flux gives 




Hence, by integration, r® 


. dF A 
= A, or -,- = -g, 
dr r^ 


where A and B are arbitrary constants to be determined to suit each 
particular problem. For example, consider the field due to a uniform 
thin spherical shell of radius a and mass m per unit area. Denote the 
potentials inside and outside the shell by F^ and Fj. The conditions 
to be satisfied by Fj and F^ are given in 4«13, namely 

A 

where r < a, V^Fj = 0, so that from (1), Fj = — + 5^; 

where r>a, V®F2 = 0,. F8=— 

where r=a, Fi = Fj and —47rm. 

* “ dr dr 
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The le«t two conditions give 

+ + .(2) 


and 


-4® 

a2 




(3). 


We also know that, for large values of r, V tends to zero like ^tifnia^/r, 
so that jBa = 0 and 4:TTina^, Hence from (3) = 0, and from (2) 

Bi = dTTtna; so that = 47rma and Fg = 47rma2/r. 


(iii) Symmetry about an axis. When the field is such that F is 
a function of r only, where r denotes distance from an axis; we may 
either quote the transformation of Laplace’s equation into cylindrical 
co-ordinates given in 1‘68, putting dVldd=0 and dV/dz = 0, or we may 
apply Gauss’s theorem to imit length of the space bounded by cylinders 
of radii r and r -f dr. In this case, the flux out of this region across the 

cylinder of radius r is — 27rr so that the corresponding flux across 

the cylinder of radius r + dris, ^nr ^ -f ^ ^27rr dr to the first power 

of dvy and since there is no matter in this region 


d_ 

dr 



or 





so that F = -41ogr+j5, 

where B are arbitrary constants to be determined to suit each 
particular problem. 

If the matter be a long uniform thin cylinder of radius a and mass m 
per imit area, we may take for the internal and external potentials 

Fi = ^ilogr-f Hi, Fgrz^glogr-f-Hg. 

When r = a, we have Fi = Fg, so that 


Ai log a + Hi = ^2 log a + Hg 


( 1 ); 


and when r = a, 


dr 


dFi 

dr 


— 47rw, so that 


a 


~ = — 47rm 


( 2 ). 




Also applying Gauss’s theorem, as in 4*23 (c), to a cylinder of radius 
>a), we have 

2rfr = —4n. 2nmaf 
dr 


so that Ag = — 47rma. 

Then from (2) Aj = 0, and from (1) Hj = — 47rma log a + Hg, shewing 
that the forms for Fi and Fg are in agreement with the results of 2*42, 
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4*6. Potential of a body at a distant point. MacGuUagh's 
Formula.* Let 0 be the centre of gravity of the body and P 
a point whose distance R from G is large compared with the 
dimensions of the body. Let m be the mass of an element 
of the body, situated at Q, where OQ=r and the angle 
QQP^e, 



The potential at P is given by 

PQ “(i?2-2i2rcos0 + r2)i 
^m/_ 2r . 


=s 


r cos 9 cos^ 9 

1 + —^ ^2+ 2JR2 


...) 


2m 2mr cos 9 — SLmr^ sin^ 9 


But if M denotes the mass of the body, A, B, C its principal 
moments of inertia at 0 and I its moment of inertia about 
OP, we have 

2m = M, liinr cos 0 = 0 (since O is the centre of gravity), 
22mr2 = ^ + j 5 + C and 2mr^ sin^ 0 = /. 


Therefore 


_M ^ + B + 
~ i? 2i2» 


is an approximation to the value of the potential. 

4*61. Attraction components at a distant point. Taking 
axes through the centre of gravity, let x, y, z be the co-ordinates 
of the distant point and r its distance. In 4*6 this means that 

♦ James MacCullagh (1809-1847). Irish mathematician and 
physicist. 
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OP is now of length r and its direction cosines are xjr, yjr^ zjr, 
80 that j ^ ^ ^^2 4 . Cz<^)jr2, 

and 

V=y + ^{iB + C-2A)x^ + {C + A-2B)y^ 

+ {A + B-2C)z^} .( 1 ), 

to this order of approximation. 

The attraction components at P are therefore given by 

dx ^ 


-^^{{B-{-C-2A)x^ + {C + A-2B)y^ + (A + B-2C)z^} 

.( 2 ). 

and similar expressions for Y and Z, 

Since action and reaction are equal and opposite it follows 
that a unit particle at P exerts on the body M forces equal and 
opposite to X, r, Z acting at P. To measure their effect on 
the body it is convenient to transfer them parallel to them¬ 
selves to act at O by introducing suitable couples. Thus the 
reaction on the body is represented by a force at 0 with 
components — X, — T, — Z, and a couple with components 


zY--yZ^ 


xZ — zX = 


Z(C^B)yz\ 
3{A — C) ZX 

I 


yX-xY^ 


3 {B—A)xy 


(3). 


4*62. As an example of 4*61 (3), if the body be the earth regarded 
as an oblate spheroid, and a distant body be regarded as a particle of 
mass Jkf' at distance r and angular distance S above the equator, we 
may without loss of generality suppose the distant body to be in the 
plane zx^ then the couple exerted by the distant body on the earth 
tending to cause rotation about a diameter (the y axis) is 


ZM'(A-C)xz 


or 


ZM'(A-G) 


sin 8 cos 8, 


and this has its greatest value when the declination 8 is 
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4*63, Mutual potential energy of two bodies at a great 
distance apart. Let Jf, M* 
be the masses and G, O' the 
centres of gravity of the 
bodies, and let 00' — r. Then 
if m'is the mass of an element 
of the second body at P and 
V the potential of the first 
body at P, the mutual potential energy is, by 3*8, equal to 





M ^ + P + 
GP'^ 20 


from 4*6, 


where Ip denotes the moment of inertia of the body M about 
OP; and this 


7i>fV A. B C — 3/ ^ f 

■^OP -ss- 


correct to the order , where / is the moment of inertia of the 
body M about 00'. 

fti' 

But is the potential at O of the body M', and by 4*6 


M' A'-^B' + C'^sr , , 

IS equal to-1---, where A , B , C , I bear 

r 2r^ 

the same relation to the body M' that A, B, C, I do to the 
body M. 

Hence we get for the mutual potential energy 

Li/* j 2r^ J 


^M' 


+ 


M(A'^B' + C'-^r) M'{A + B+C--3I)' 


2r® 


+ - 


2r3 


']• 


4*7. Gentrobaric bodies. When the action of terrestrial 
or other gravity on a rigid body is reducible to a single force 
always acting through a point fixed relatively to the body, 
that point is called the centre of gravity of the body and the 
body is said to be centrobaric. 
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It follows that uniform spheres and spherical shells are 
centrobaric. In elementary statics it is usual to assume that 
the attractions of the earth on the elements of a body are 
a set of parallel forces with a resultant acting through a point 
fixed relatively to the body, so that on this hypothesis every 
body has a centre of gravity and its position is defined by 


relations 


Mx — Ymx^ My —limy, Mz — lmz 


( 1 ). 


We shall now prove the following theorem: If the resultant 
attraction of a body at all external points passes through the 
same fixed point, that point must be the centre of gravity, as 
defined by {!), and the body must be such that every axis through 
the centre of gravity is a principal axis and all principal moments 
of inertia are equal. 

Take the fixed point as origin 
0. Let fjL be an element of mass 
of the body at Q whose co¬ 
ordinates are x, y, z. Let P be 
an external point at distance r 
from 0 and let I, m, n be the 
direction cosines of OP. Then 



for the potential at P we have 


^ _ fX 

PQ {{It — xf^ -f {mr — yY + {nr — zYy 

2{lx-\-my-\-nz) + 

-r - 

_M 'LiJ.{lx+my + nz) 1’Efi {x^ + + z^) 


3 2/x (lx+my + nz)^ 
+ 2 ^ 

M , M(lx+my + nz) , 1 ^ , 

= 7 -+ . - +:^{P^ii(2x^-y^-z^)+ 

+ 6mrtS)a3/2J+ ... 4- 


-f %Dmn 4*... 4“. *.} 4-..., 
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where A, B, C, D, E, F are the principal moments and the 
products of inertia. 

The resultant force at P is to pass through 0 for all positions 
of P, so that, if we write Z = cos m = sin 0 cos </>, n — &mdsin<t>, 
we must have dVjdO and 3F/3<^ both zero for all values of 6, (/>. 
This can only be so if the coefficients of the different powers 
and products of I, m, n vanish separately, so that 

= ^ = 2 = A^B^C and P = P = P = 0, 
which proves the proposition. 


4*8. Points of equilibrium. If in the presence of attract- 
ing bodies the resultant force on a free particle is zero, the 
position of the particle is a point of equilibrium. If V denotes 
the potential at the point P of co-ordinates x, y, z, then P 
is a point of equihbrium if 


dx dy dz 


( 1 ), 


for these are the conditions that the resultant force at P is 
zero. 

Regarding F as a function of x, y, z, the equation 

F {x, y, z) == const.(2) 

represents a family of equipotential or level surfaces. Using 
suffixes to denote partial differentiation, the direction cosines 
of the normal at (x, y, z) on such a surface are proportional to 
Vx,Vy,V^, which shews that the resultant force at the point 
is normal to the surface (2*26). It follows that if the point is 
a point of equilibrium, it is a singular point on the equipotential 
surface; i.e. there is not a unique tangent plane at the point, 
but every straight line through the point cuts the surface (2) 
in two coincident points. 

Let F' be the potential at a neighbouring point P' of co¬ 
ordinates x + i, y-\-7], z + C, then 


+ +.(3), 

where the derivatives denote values at (x, y, z). 
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The result (3) incidentally affords another proof of the 
theorem 4*3 (i) that the potential cannot have a maximum 
or a minimum value at any point unoccupied by matter. For 
a first condition that V' may be greater (say) than V for all 
small values of rj, ^ (positive or negative) is the vanishing 
of 1^, 1^, and the further conditions require that Vyy, 
shall all be positive; which is not possible since VW = 0. 
At a point of equilibrium there will therefore be some directions 
in which V'>V and others in which F' < V, so that equilibrium 
would be stable for some displacements and unstable for 
others [4*3 (i)], and these will be separated by the surface for 
which F' = F; and from (3), since = = = this is, 

for small values of a quadric cone 

+ Vyyri^ + F,,^2 ^ 2VyM + 2F,,C^-f 2V,yir^ = 0 .(4). 

Further, since V^^-i-Vyy-i-T^^=0, it is a cone with three 
perpendicular generators. 

If two sheets of the same level surface intersect one another 
along a line, every point of this line is a point of equilibrium, 
because there are two normals to the surface at such a point 
each of which would be the direction of the resultant force 
if any. In such a case the tangent cone to the surface, repre¬ 
sented by (4), breaks into two planes; so that (4) is equi- 

valent to + my + n^) (I'i + m'ff + n'Q = 0, 

where ll' + mm' + nn' = f^x + Vj,^ + 1^g = 0, 

so that, if two sheets of a level surface intersect, they intersect 
at right angles. 

4-9. Examples, (i) A particle is placed on one of the plane faces of 
a uniform circular cylinder at a small distance from the centre of that 

fa>ce; prove that it will make small oscillations of period 2 (^^ + ^*)*» 

where a is the radius of the cylinder^ h its height and p its density, and p. 
is the attraction between two unit masses at unit distance apart, 

[M. T. 1889] 

The attraction of the cylinder at an external point on its axis is 
foiind from 2*52 to be 

2'rrpp{h- V[(^ + «)® + o®]+ V(25* + a®)} .(1), 

where z is the distance of the point from the centre of the nearer end. 
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Taking the centre of the plane face as origin and the axis of x 
through the particle and the mass of the particle as unity, its equation 
of motion is 

•• 

.(2). 


where dVjdx is the force of attraction at (a;, 0,0). By MacLaurin’s 
theorem 

dV /dV\ ^ (dW\ . 

= Ui jo)/h.gher powers. 

{dV\ fdW\ fdW\ 

And at the origin ^ — j = 0, by symmetry, and ‘ 


/dw\ , 

fd^V\ , 






therefore 




dx 


'\dx^/o 






•(3), 


where the expression (1) represents --dVjdz at (0,0, s). Hence, on 
the axis 


dW 

dz^ 


= 27r^p |- 


h-\-z 


o»)}’ 


(8W\ _ _ 2„^ph 

\dz^J,~V(h‘+a^) . 


and at the origin 
Hence, from (2), (3) and (4), the equation of motion is 


x + 


‘jryiph 


X=0y 


giving oscillations of the period stated. 

^ (ii) A solid bounded by the plane xy and the half of the ellipsoid 
x^ld^ + y^/b^ + z^/c^ = 1 o/i the positive side of the plane has a volume 
density pz(lla^+ 1/62 + 3/0^). Prove that it produces the same field of 
attraction outside itself as would a layer of surface density 

ip.(l-x^la^-y^/b^) 


over the plane boundary and pzjp over the curved surface, when p is the 
central perpendicular on the tangent plane. [C. 1928] 

To prove the equivalence of the two fields of attraction in external 
space, it is only necessary to shew that if we reverse the sign of one 
distribution then the two distributions together produce a zero field 
at external points. 

The potential in this combined field may therefore be assumed to 
be zero outside the solid and some function inside the solid which 
vanishes over the surface. We therefore write 


V — Az{l-~x^/a^’-y^jb^—z^lc^), inside, 
and F' = 0, outside. 
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Hence we get 


?Z- 

dx " 


2Azx 2-42^ 


6* 


and 


dy 

2Az W 
a® * dy^ 

ifp is the volume density, we have 


^2 


a® 


dx^ ~ 


2^2 W 
'6* ’ ^2^" 


6^2 


or, putting ^ = 2^7r, P =(^2 + ^2 +inside,! 


and 


p = — -p- V^F' = 0 outside 

47r 


.(1). 


Again differentiating along the outward noimal, on the plane 
boundary we have 

47ro’ = 


and, putting A = 2/x7r, 


/ dV 

dV') 


'^\ dz 

dz j 


= -A fl 

x^ 

y\ 

V 


W 

- 

a* 



Also, on the curved boundary 

^ dV dV' 
^ dn 


__ /xdV ydV zdV\ 

\a^ dx~^b^ dz) * 

which reduces to — 2Azlp9 so that, putting A = 2p7r, we have 

o^-yzjp .(3). 

Since the distributions (1), (2), (3) give a zero field outside the body, 
it follows that when the sign of the surface distribution is changed it 
produces the same external field as the volume distribution. 


EXAMPLES 

1. If the potential in the interior of an attracting mass bounded by 
two concentric spherical surfaces be proportional to l/Vr, where r 
is the distance from the centre, find the law of density, and the potential 
exterior to the attracting matter. 

Shew that at any point of space where the density is p and potential F, 
the magnitude of the force is given by the expression 

:i{8^pF+V«F>}i. 


[C. 1892] 
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2. Find the distribution of matter which produces a potential 

/Lt (x® + 2/® + — 2ax) at all points inside the sphere x® 4- 2/® -f* 2^® — 2ax = 0 

and zero potential at all points outside it. [C. 1908] 

3. Find the distribution of matter, none of which is external to the 

sphere r =;: a, which produces a potential — 27rpx^ at points within or 
on the sphere, the origin being the centre, x the abscissa, and p a 
constant. [C. 1914] 

4. What distribution of matter will produce potential 


ox\ 
3rV ’ 


according as r is less or greater than a? [C. 1902] 

5. The potential outside an infinite elliptic cylinder 
x*/a* 4- 2/^/6® —1 = 0 

is zero; inside it is F = x^/a® 4- — 1. Find the distribution of matter. 

[C. 1915] 

If the potential within an attracting medium having spherical 
^''^junmetry is given by V = (c^ where r is the distance from the 

centre and c a constant, prove that the density is given by 
p=3cV47r(c2 4*r2)l, 

and that the potential energy of the whole medium is — 37r/32c. 

[M. T. 1915] 

7. Find the distribution of matter that gives rise to potentials 

Fi = JgTra® 4-1^5/X7ra( 2 x^ — 2 /^ —r<a 

arid Vi = ^tiiT^-+^-^IJi,iT^{2x^-y^-z^), r>a. 

8. Prove (by taking x^ + 2 /^ + 2 :^ — 1 as a potential function) that, if the 

space enclosed by the surface x* 4* 4- 2 :^ = 1 is filled with attracting 

matter of density proportional to the square of the distance from the 
origin, its attraction at all outside points will be the same as if its mass 
were distributed over the boundary ernface with surface density 
inversely proportional to the distance of the tangent plane from the 
origin. [M. T. 1890] 

^ 9. Shew that the family of surfaces defined by 

a;* 4 - 2 /* = const. 

can be a family of equipotential surfaces in free space, and find the 
law of potential. [C. 1906] 

x/io. Shew that the system of co-axial cylinders 
X* 4* 2 ^* 4-2Ax 4-c® = 0 

can form a system of equipotential surfaces. [London Univ. 1934] 
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^1. Shew that the cylinders a:*-f y*=:2Aar are a possible set of eqni- 
potential surfaces in empty space, but that the spheres a;® + = 2Aa; 

are not. 

In the former case find the potential as a function of A. 

[London Univ.] 

1 2. Two equal circular cylinders, with their axes in the same straight 
line and their centres of gravity at a great distance r apart, attract 
one another according to the law of the inverse square of the distance. 
Shew that the resultant attraction is approximately 

M^r-^ + 6M(A-C)r~\ 

where M is the mass of one of the cylinders, C is its moment of inertia 
about its axis and A its moment of inertia about a perpendicular line 
through the centre of gravity. 

13. If two distributions of matter, A and B, have throughout 
a region T, to which both are external, the same equipotential surfaces, 
the attractions due to A and B throughout the region T are in the 
same direction and in a constant ratio. 

14. If two distributions of matter, A and B, have the same closed 
surface B, enclosing both, as an equipotential surface, every surface 
outside S which is an equipotential surface for A is also an equipotential 
surface for B. 


15. A finite solid of revolution is cut at right angles by a plane, 
and a particle is placed on the plane, a smooth surface, near to its 
centre O. Shew that the time of a small oscillation is 



where F»is the potential of the solid, the axis of the solid is the axis 
of z and dWjdzQ^ denotes the value at O. 

16. A gravitating solid of revolution is cut by a plane perpendioulax 
to the axis. A particle is fastened by a fine string of length Z to a point 
in the prolongation of the axis, so that when the string is vertical the 
particle just does not touch the plane face at its centre 0. Assuming 
the conditions such that when the particle is slightly disturbed the 
motion is that of a simple pendulum, shew that the time of a small 
oscillation is 

^"(sTpR') ’ 

where R is the force exerted by the solid on unit mass at O and R' 
is the derivative of that force at O, taken outside the solid, along the 
axis. [C. 1892] 
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17. A uniform circular disc is of mass M and radius a and has its 
centre of gravity at O. A particle of mass m is situated at P, at a great 
distance r from 0, and OP makes an angle B with the axis of the disc. 
If the attractive force of the particle on the disc is reduced to a force 
through 0 and a couple, shew that the moment of the couple is 
sin 26)/r^y and find the magnitude of the force. 

[London Univ. 1938] 


18. Prove that the a;-component of attraction of a gravitating 
medium on a small sphere of itself with its centre at (x, y^z) is approxi¬ 
mately 


1 


8V^ 




30y 




S'-)- 


where a is the radius of the sphere, p the density and V the potential 
at (XyyyZ), [C. 1932] 


19. Shew that, if the curves /(ir,^/, A) = 0 form a system of equi- 
potential lines in free space for a two-dimensional system, the surfaces 
formed by their revolution round the axis of x cannot be a system of 
equipotential surfaces in free space unless 


IdX 

y^y 




is a constant or a function of A. 

^ 20. Shew that the family of cylinders 


(x^ -f 2/2)3 _ 2^3 __ ^xy^) -f a® = const. 


[C. 1905] 


is a possible form of equipotential surfaces, and find the corresponding 

[London Univ. 1938] 

that (a? —c)2 + 2 /* = A[(a;-f c)2-f 2 /^] represents a family of 
equipotential surfaces, and that 



V = AlogA + P, 


where V is the potential function and A and B are arbitrary constants. 

[London Univ. 1939] 


© 


22. Determine a distribution of matter whose potential shall be 
con^ant over confocal elliptic cylinders. [C. 1915] 

3. Shew that the family of ellipsoids 
oH-A 6-f-A c + A 


where A is a variable parameter and a, 6, c are constants, is a possible 
form of equipotential surfaces, and express the potential in terms of 
a, by c and A. 


R N A 


7 
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24, The potential functions in the three regions bounded by the 
spheres r = 6 and r = a (> 6) are 

= + r>a; 

F, = .^jiyp[^|o*-Jr*-* + a>r>h-, 

and F 3 = jj 7 Typ[^a* —|6®-f i(22j2 —a;* — i/2)], 6>r>0. 

Find the volume and surface distributions of matter and the total 
mass of the distribution. [London Univ. 1939] 


26. In the interior of a spherical distribution of matter of radius a 
^e gravitational potential at a distance r from the centre is 

A sin (irrla ) 

Trrla 


- 4* const. 


Prove that the density at any point is 


ttA sin(7Tr/a) 


, and that the mass 


nr /a 

of the whole distribution is AajG^ O being the constant of gravitation. 

Prove also that the negative potential energy (compared with a state 
of infinite diffusion) is ^A^a/O. [M. T. 1926] 

26. Find the attracting systems whose potential is 

■^i^7rxyz{9a^-‘lr^) + 4nrx/{x^ + y^) for r<o, 


and 


^TTxyz , 4- 4naxl{x^ 4- y^) for r > a, 


where r* = a;® 4- 2/^ 4- 


[C. 1927] 


27. The axes being rectangular, determine what volume and surface 
distributions of matter will give rise to potential Kxyz (a — x — y -- z), 
where K and a are constants, at all points within the tetrahedron 
bounded by the co-ordinate planes and the plane a? 4- 2 / 4-2 = a, and 
potential zero at all points outside the tetrahedron. 

Shew that the total mass of the matter contained within the tetra¬ 
hedron is Ka^lSOny. [London Univ. 1938] 

28. An infinite uniform straight line occupying the positive part 
of the axis of z attracts according to the law of nature. Shew from 
elementary considerations that its potential is a function of r — z, 
where r is distance from the origin, and hence prove that its potential is 

clog (r~z)4-c'. 

29. Prove that, if the matter be distributed symmetrically about 
jbhe axis Oz, dWldzdr=:0, at eve^y point of this axis, where r is the 
distance from this axis. Prove" also that if the origin be a point of 
equilibrium, the lines of force near it are given by the equation 

zr* = const. [C. 1901] 
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30. Compare roughly the couples exerted on the earth by the sun 

and moon respectively, supposing them to be at the same angular 
distance from the equator and the orbits of the moon and earth to be 
circles. Take the ratio of the mass of the moon to that of the earth 
to be 1/80. [M. T. 1922] 

31. Obtain^the mutual potential energy of two gravitating bodies 
placed at a considerable distance from one another, neglecting powers 
of the inverse distance higher than the third, and obtain the specifica¬ 
tion of the moment, about the centre of mass of either body, of the 
attraction exercised upon that body by the other. 

Shew that the couples exercised by the sun and the moon on the 
earth neutralize one another at an instant when they are on the same 
meridian if 

sin2ai'_ S(E-\^M) (T'y 
srn2a> ~ M(E-\-S)\t) * 

where w and ay are the angular distances of moon and sun from the 
equator, E, M and S are the masses of the earth, moon and sun, 
T' and T being the lengths of the sidereal month and year respectively. 
The orbits of the moon and earth may be considered circular. 

[M. T. 1937] 


32. O is the centre of a uniform cube of mass M and P is an external 
point. The distance OP is large compared with the edge a of the cube, 
and the direction cosines of OP referred to three concurrent edges of 
the cube are Z, m, n. Shew how to express the gravitational potential 
at P in a power series in OP~^, and evaluate the potential in terms of 
OP, ilf, I, m, n and o, rejecting all terms in OP~^ of higher order than 
the fifth. [P. 1936] 


33. Shew that the equipotential surfaces of a uniform equilateral 
triangular lamina (of side 2a), which are at a great distance from the 
lamina, approximate to the surfaces of revolution obtained by rotating 
the curves 


1 a2(l-3cos2^) 

12r» 


= const. 


about the initial line. [London Univ. 1938] 

34. Prove that, in the case of two equal cubes at a great distance 

apart, the third order term in the expression for their mutual gravita¬ 
tional potential energy disappears and obtain the next non-vanishing 
term. [M. T. 1929] 

35. Prove that for two gravitating solids, if one is a fixed sphere 
and the other has unequal principal moments of inertia and is free 
to turn about its centre of gravity, the only positions of equilibrium 
are those in which a principal axis at the centre of gravity of the 
second body is directed towards the centre of the sphere. [P. 1933] 

7-2 
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36. Prove that a distribution of matter of surface density mpl^Ttr^ 

over an ellipsoid will produce the same potential at every external 
point as a homogeneous concentric sphere of mass m within the 
ellipsoid, r and p being the central radius and the perpendicular from 
centre on the tangent plane. What is the potential of this distribution 
inside the ellipsoid? [C. 1892] 

37. In a region of space which does not intersect the plane t/ = 0 
a fimction ^ satisfies V2^ = 0, and the surfaces with the equations 
<!> = constant are portions of the elliptical cylinders 


where A is a variable parameter. By transforming to co-ordinates 
Tj, f defined by 

a:; = a cosh ^ cos y = asmh^ sin rjy 2 =^, 

^>0, 27r>^^0, 

prove that 0 is of the fonn where A and B are constants. 

[P. 1937] 



Chapter V 

GREEN’S THEOREM 


5*1. In 1*3 and 1*31 the following theorem was established: 
If Vy V* functions which unth their first derivatives are 
finite^ continuous^ and single-valued through a singly-connected 
region"^ bounded by a surface S, 


J \ 3a: dx dy dy dz dz j 


=^V^-^d8~ jvV^V'dv 


VWdv .(1), 


where the surface integrals are over the boundary S and the 
volume integrals throughout the bounded region and djdn denotes 
differentiation in the direction of the outward normal. 

The proof in 1*3 was for a region bounded by a single 
surface, but it requires Httle consideration to see that it is 
also true for a region with several 
boundaries, such as that bounded 
externally by S and internally by 
Si and S 2 as in the figure. But in 
this case the surface integrals must 
extend to the three surfaces S, Si 
and /Sg, and care must be taken 
that the normal difierentiation is 
outwards from the region of volume integration over each 
surface, as indicated by the arrows. The theorem needs 



* A singly-connected region is one in which every closed curve can 
be contracted to a point without passing out of the region. The space 
between two concentric spheres is singly-connected, but the space 
inside an anchor ring is doubly-connected. 
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modification for multiple-valued functions and multiply- 
connected space,* but we are not concerned here with such 
modifications. 

In the applications of this theorem which we shall consider, 
V and F' wiU be potential functions of attracting masses. 
We know that the potential function F of a system of attracting 
particles becomes infinite at each particle, if we regard the 
particle as a concentration of matter at a point; also, that 
there is a discontinuity in the normal derivative of F across 
any surface on which there is a surface density of matter 
(3*7), consequently such points and surfaces of discontinuity 
must be excluded from the region of integration in applications 
of the theorem. But with this proviso the region might be 
regarded as bounded externally by an infinite sphere and 
internally by one or more surfaces adapted to the particular 
problem under consideration. In such a case integration over 
the infinite sphere must be included among the surface 
integrals in (1). But if F, F' are potentials of finite masses 
within a finite region, and Jt denotes the radius of a large 
sphere, then on this sphere F, F' are each of order l/H; 
dVjdn, dV'Idn are of order l/R^ and dS" is of order R^dco, 

r 0F' 

where doj is a solid angle. Hence F is of order l/R 

and may be neglected when R tends to infinity. Similarly 

r dv 

for V V' ^ dS, so that when F, F' are potentials of finite masses 

the surface integrals over the infinite sphere vanish. 

It is evident that very varied applications of Green’s 
theorem can be made, but in every case the first essential is 
to specify to what region of space the theorem is applied and 
to ensure that F, F' and their derivatives are continuous in 
the way required within that region. 

5*2. Applications. Gauss’s Theorem. Put F' = 1, and let 
the region be bounded externally by a single surface 8, and 
suppose it to enclose matter of finite volume density p, where 


♦ See the author’s HydrodynamicB, 4*81. 
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gauss’s thboebm 


p may be zero through any part or parts of the region. The 
theorem gives 




VWdv .( 1 ). 


But, by Poisson’s equation, V*F = — ^irp, wherever the matter 
'dV 


exists; therefore 
or Ji?„d;S = 


/ 


■dS = 


8.- 

477 (mass enclosed hy 8) .(2); 


which is Gauss’s theorem. 

It is instructive to see how the apphcation of Green’s 
theorem needs modification when the surface 8 encloses not 
only matter of finite volume density p, but also material 
particles, implying points at which F is infinite, and surfaces 
on which there is a surface density, implying discontinuity 
in dVjdn. 

It will be sufficient for our purpose to suppose that the 
surface 8 encloses a body of mass M whose density at any 
point is p, a particle of mass m at a point P, and a surface 8' 
on which there is matter of 
surface density a. In order to 
apply Green’s theorem to the 
region enclosed by 8, we must 
exclude the point P and the 
surface 8' from the range of 
integration. This is done by 
drawing a sphere S of centre P 
and small radius €, and by 
drawing surfaces parallel to 8' 
and close together on opposite 
sides of 8' as shewn in the figure. If these surfaces are 
sufficiently close together, end effects may be ignored. We 
now apply Green’s theorem in the form (1) above to the region 
bounded externally by 8 and internally by S and the surface 
which surrounds 8\ We have 



f l^ds+r (~+^]dS'=hwdv 

Json Jxon Js'\oni dnj J 


•( 3 ). 
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where in the second integral 'bn is directed out of ‘ the region 
i.e. into the small sphere, and in the third integral bn-^ and bn^ 
are both directed towards S\ 

In the immediate neighbourhood of P the part of V which 
depends on the particle m at P is the only part which will 
contribute to the surface integral, if € is taken small enough; 
and the second integral is therefore 


r m 

J - ^ 


€^d(i) = 4:TTm. 


Then in the third integral bVjbn^ and bVjbn^ are the attrac¬ 
tion components normal to 8\ and by comparing with 3*7 (1) 

^ . ... . 
and paying due regard to sign we see that 


so 


that the third integral is 4??Hence (3) is equivalent to 

j ^ dS + ^irm -f 4:7r jadS' = — 4r7T jpdv, 

J* R^dS = — 47rilf — 4:7Tm — 4:'n^adS' 

= — 47r (mass enclosed by S), 


5*21. Green’s Reciprocal Theorem. Let F, F' be the poten¬ 
tials of two different distributions of matter of finite volume 
densities p, p situated anywhere in the finite region of space. 
Apply Green’s theorem 5*1 (1) in the form 

to the region bounded by the infinite sphere. As explained in 
5*1 the surface integral over the infinite sphere vanishes, and 
since by Poisson’s equation VW = — 47rp and V^F'= — 47rp', 
therefore (1) reduces to 

jVp’dv^jv'pdv .( 2 ). 

This is the reciprocal relation obtained in 3’8, either side of 
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the equality representing the mutual work of the two distribu¬ 
tions of matter. 

We might infer that if the distributions of matter include 
separate attracting particles and surfaces on which matter 
is condensed, then the complete expression of the theorem will 
be a relation 

jvp'dv-\-I.{Vm')+i:jva'dS' 

V'pdv + 'L{V'm) + 'Ljv'adS .(3), 

i 

where m, m' denote the masses of typical particles of the two 
distributions and adS, ad S' elements of the surface distribu¬ 
tions and the 2's imply summation for all such particles and 
surfaces. But though this is a true result, the foregoing 
argument does not constitute a proof of it, because the region 
of integration in (1) (all space) now includes the particles, 
i.e. points at which V or V' becomes infinite, and the material 
surfaces across which dVjdn or dV'jdn is discontinuous, and 
we may not apply (1) without excluding such points and sur¬ 
faces from the region of integration, so that the result (2) 
establishes nothing concerning the case of particles or material 
surfaces. 

Let be the masses of typical particles of the two 

systems A, A'. Surround them 

by spheres of small radii 

, Cl' and centres A, A', Let S, S' 
be typical surfaces with surface 
densities <t, a' in the two systems 
and let them each be enclosed by 
parallel surfaces a small distance 
apart, as in 5*2. We can now apply 
the theorem (1) to the region bounded externally by the 
infinite sphere and internally by the small spheres and other 
surfaces so as to exclude the particles and material surfaces 
from the region. The surface integral on the left in (1) must 
now be understood to refer to the infinite sphere and to all 
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the inner boundaries, and, as before, over the infinite sphere 
it vanishes. Consider 





over the small sphere Sj. 

Since is small, on the principal part of F is 


SO that we may put - -- + 17 for F, where 77 is finite; and, taking 

m dV 

account of the sign of we may put --| +171 for ^, where 771 

Ctl 

is finito; and dSi = so that the integral is 


and if we make -> 0 , this reduces to — 47 r^'mi, where 
is the value of F' at -4. In like manner the same integral 
over the small sphere Sg reduces to 

Taking the same integral over the surfaces which surround 
S (neglecting edge eiffects), we have 





dS, 


and taking account of the signs of dn ^, dn ^, we see that by 3*7 

jy , ,0F' ar ^ , V, . 

;r—— = ‘iTTa; and - 1 - 5 ;— = 0 since by hypothesis none of 

v^2 VTh-^ ^^2 

the second distribution of matter resides on S and F' is 
continuous across 8. The integral therefore reduces to 


— 47 rjF'adAS. In like manner the similar integral over the 

surfaces surrounding S' reduces to 47 rjFcr'd!iS'. 

Then substituting from Poisson’s equation on the right-hand 
side of ( 1 ), dropping the factor 4?!, collecting like terms, and 
summing for aU particles and material surfaces, we obtain 
the result (3). 

We conclude that whenever we use Green’s theorem to 
prove something about attracting particles or surfaces it is 
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not through the volume integrals but through the surface 
integrals that we must expect their masses to enter into the 
result. 


5*3. Potential energy. Let V be the potential due to 
a system of particles, surfaces on which matter is condensed 
and finite bodies. It will suffice for our purpose to consider 
a single particle of mass m at a point P, a single surface 8' 
on which there is matter of surface density a, and let the finite 
bodies be represented by a distribution of volume density p. 
Take the infinite sphere as the external boundary of the 
region and the same internal boundaries as in 5'2, excluding 
from the region the point P and the surface /S'. 

Apply Green’s theorem in the form 5*1 (1), putting F' = F. 
For the reason given in 5-1 the integral over the infinite sphere 
vanishes, so that we have in this case 




VVWdv 


( 1 ), 


where the volume integrals are through all space, save the 
vanishing regions which we have excluded, and the surface 
integrals are over the small sphere round P and over the 
surface /S' as in 5*2. 

The integral on the left is where R is the resultant 

force at a point of the element dv, 

r 171 

The first integral on the right, as in 5*2, is F-g €^dco, where 

if € is small enough we may take for F its value at P, say 
Vp , so that the integral is 47rmlp. 

The second integral on the right, as in 5’2, is 47r Ford/S', 

J S' 

and the third integral on the right is 4:7rjVpdv, integrated 
wherever there is a volume density of matter. 
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Therefore 


jR^dv~47rmVp + 47Tj VadS' + Airjvpdv .( 2 ). 


It is clear that the presence of other attracting particles and 
surfaces would contribute other terms hke the first and second 
on the right. 

But from 3*81 the potential energy of this attracting system, 


there denoted 


is in this case 


by - */ Vdm, is 
-\{mVp + j VadS' +j Vpdv^, 

and, from ( 2 ), this is equal to xR^dv integrated through 

877 j 


all space. It follows that the potential energy of any attracting 
masses may be represented by an energy density of amount 
— R^IStt distributed through all space. 


5*31. If V is a solution of Laplace's equation which with its 
first and second derivatives is single valued and finite within 
a region bounded by a closed surface S, and dVjdn^-O over the 
surface, then V is constant throughout the region. 

In 5-1 ( 1 ) put F' = F, so that 

Since dVjdn — O over S, and 'SIW — 0 inside S, therefore the 
integral on the left of ( 1 ) is equal to zero. This requires that, 

dV dV dV 

at every point of the region, — = — = — = 0 , so that F must 

be constant throughout the region. 

It follows that if are two solutions of Laplace’s 

equation satisfying the above conditions within the bounded 
region and making dVJdn^dVJdn over S, then T 1 — I 2 also 
satisfies the same conditions and makes 3 (Ij — ^ over 

S, so that Ii —T 2 iRRsf constant throughout the region. 
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5*4. Green's equivalent layer. Let be a closed equi- 
potential surface of a given distribution of attracting matter 
of which some or all lies inside 8. For the sake of brevity 
we assume that the matter inside 8 can be represented by 
a volume density p. Let r denote distance from a point Q 
outside 8. Apply Green’s theorem to the region inside 8 in 
the form 


jVj^dS-jvVWdv=jv'^dS-jv'VWdv ...( 1 ), 


where V denotes the potential of the whole distribution of 
matter and F' = 1 /r. 

This gives 


-jv^^^dv= j^^dS-j^VWdv ...( 2 ). 

Let Vg be the value of F on the equipotential surface 8, 
then the first integral in (2) is Vg is 

potential due to a particle of unit mass at Q, and Q is outside 8, 
so that by Gauss’s theorem J d8 — 0; therefore the first 

integral vanishes. The second integral also vanishes since 
^ = 0 throughout the range of integration. Also V^F = — 47rp 
throughout the same range, so that (2) reduces to 


But 


j 




r dn 


d8 + 47r 


r pdv 


pdv 


is the potential at Q due to the matter inside 8, so 


that denoting this by (F)^ we have 




.(3). 


Therefore the part of the potential at any point outside 8 
which is due to the matter inside 8 is the same as would be 
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produced by a layer of matter of surface density 

1 dV 


4-77 dn 


.(4) 


spread over 8, where V is the potential due to the whole of 
the attracting matter. 

Further, the amount of matter in the layer is given by 


J' 


GdS = 


rw 


dS 


47 r J dn 
= 5*2(1), 

=jpdv 

= mass contained within S, 


We have therefore proved that the matter contained within 
any closed eqxiipotential surface in a given field can be 

spread over the surface 8, with a surface density — ^(9F/9n) 

47r 

at any point, without altering the potential at any point of 
the field outside 8, This distribution is known as Green's 
equivalent layer. 

If we suppose that this distribution over the surface 8 of 
the matter inside 8 has been effected, 8 is now an equipotential 
surface containing no matter, so that the potential at any 
point inside 8 is now Vg the value on the surface [4-3 (ii)]. 

This can also be established analytically thus: let r' be the distance 
from a point Q' inside S and let p' denote the density of the matter 
outside S, We apply Green’s theorem in the form (1) to the region 
between S and the infinite sphere, taking V as before to be the potential 
of the whole distribution of matter and 1/r'. The normal deriva¬ 
tives in (1) are now directed along inward normals to dSf so if dn is 
to retain the same meaning as in (3) and (4) we must change the signs 
of the surface integrals and write 

-/f !,*= -/i VF*...(6). 

The first term is — ®^d since r' is measured from a point 

Q' inside 8, --4it, (Gauss), and the first term in (6) is 
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47rVt, The second term vanishes as before, and the fourth term is 
so that 




Or, if we put a for the surface density of Green’s layer, as in (4), 


[adS fp'dv 
} r' r' 


( 6 ). 


We thus verify that the potential inside S due to the joint effects 
of Green’s layer and the original matter outside S is the constant 
potential of the equipotential surface. 


5*5. Dirichlet’s problem.* There exists a solution of 
Laplace's equation which with its first and second derivatives 
is single-valued, finite and continuous throughout a region 
bounded by a surface S and takes arbitrarily assigned values over 
the boundary 8. 

Dirichlet’s method of investigation is as follows: let V be 
any function which satisfies all the conditions save that it is 
not necessarily a solution of Laplace’s equation, then we can 
prove that, of all such functions, that one which gives a 
minimum value to the integral 

through the bounded region is a solution of Laplace’s equation. 

Let Vi be the function which gives a minimum value to W, 
If V be any function which satisfies the same conditions of 
continuity and takes the same values over 8, then F = 
where 1^ = 0 on 8, Conversely, let satisfy the same con¬ 
ditions of continuity, and let P 2 = 0 on 8, If e is any constant, 
then + cPJ takes the assigned values on 8 and satisfies the 
conditions of continuity. Thus the set of functions to be con¬ 
sidered are comprised in the formula PJ + €p^, where P^ = 0 on /S. 



. 

* Peter Gustav Lejeune Dirichlet (1806-1869). German mathe¬ 
matician. 
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^ince Wi is the minimum value of TF, W— Wi must be 
positive whatever be the sign of €, so that the second term 
on the right of (1) must vanish. But by Green’s theorem the 
coefficient of 2€ 

^jv,^^ds-jv,vw,dv, 

and T^ = 0 on /S, so that the first integral vanishes and 

therefore /• 

lV^V%dv = 0 

for all functions I 4 . We must therefore have = 0 through¬ 
out the bounded region; for if had a sign throughout 
any portion of the region we could take to have the same 
sign and the integral would not vanish. 

The remainder of the argument is an assertion that among 
the different functions which satisfy the conditions of con¬ 
tinuity and take the assigned values over S there must be 
one which gives W a smaller value than the rest, and, from 
what precedes, this is a solution of Laplace’s equation. But 
the argument is invalid, because Weierstrass has shewn that 
we are not entitled to assume the existence of a minimum, 
but only the existence of a lower limit to the possible values 
of W and we cannot assert that the limit is attained. 

Many mathematicians have produced solutions of Dirichlet’s 
problem but they are too lengthy for reproduction here.* 

The theorem is also applicable to the region between the 
surface S and the infinite sphere if we impose the condition 
that the required function shall vanish like 1 jr at infinity. 

We are entitled to conclude from the argument given above 
that if there is a solution it is the only one. For the above 
argument shews that, if = 0 and 1^ = 0 on 8, then the 
coefficient of e in (1) is zero, so that W = €^^ 2 * Similarly 

if V2F = 0, then since Pi==F —epg, we have, by the same 
argument, = TF-f Hence e^W^ = 0, and so €1^ = 0, i.e. 
F = and the solution is unique. 

♦ See H. Poincar6, Theorie du Potmtid Newtonien; and, for a 
bibliography of the subject, P. Appell, Traite de Micaniqm Eationelle, 
U ni, p. 98. 
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5*6] 

5'6. To find the solution o/ V^F = 0 inside a region bouhded 
by a closed surface S when the values of V and dVjdn over S 
are known. 

Take Green’s theorem in the form 


j^V^-V'^^dS=j{VVW'-V'VW)dv ...( 1 ). 


Put F' = 1 jr, where r is the radius 
vector drawn from a point P inside S. 
Take a sphere S of small radius e 
with P as centre and apply the 
theorem (1) to the region bounded by 
8 and S. Then, since (1/r) = 0 and 
V2F = 0 throughout the region, the 
right-hand side of (1) is zero, and 
therefore 





dn 


IdV 
r dn 




Taking account of the inward direction of dn over S we see 
that d ( 1 /r)ldn = 1 /c^, and since we may put dS = e^dw, we have 

=jvdo), and this being independent of e we may 

take € as small as we please and write for the integral 4cttVp , 
where Vp denotes the value of F at P. 

r 1 dv r 1^F 

Again, dS= -T^^e^dw -> 0 as €->0; so that (2) 

J^rdn J^edn 


reduces to 




IdV 

dn 


^ d(l/r) 

dn 


dS 


.( 3 ), 


giving the value of F at any point inside 8 in terms of the 
values of F and dVjdn on the surface. 

When the region is boimded by the surface 8 and the infinite 
sphere, it is easy to shew that the same theorem holds good 
provided we impose the condition that RV and dVjdn 
remain finite as R tends to infinity, where R denotes distance 
from an origin in the neighbourhood of 8. 


RN A 


8 
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[ 6 - 61 ~ 


5*61. Green's function. We leam from Diiichlet’s problem 
that there should be a unique solution V of Laplace’s equation 
which with its first and second derivatives is finite through 
a bounded region and takes a prescribed form over the 
boundary, but the procedure of 5*6 seems to imply that we 
need to know not only the value of the function V but also 
the value of dVIdn over the boundary in order to determine V 
at points inside the region. The problem can however be 
reduced if we can find a function of a more restricted form 
known as Green's function and defined as follows: 

Let if be a solution of Laplace’s equation inside the boun¬ 
dary 8, which takes the value — 1/r on 8, where r denotes the 
distance from an assigned point P inside 8, 

Let G = H + l/r; then G is called Green's function for the 
point P and surface 8 ; and by definition we see that G vanishes 
on 8 and satisfies Laplace’s equation at every point inside 8 
except at the point P where it becomes infinite. 

In 5*6 (1) put V' = Hy then since V^F = 0 by hypothesis and 
== 0, we have 




adding this, with a factor ~ 1/477, to the right-hand side of 
5*6 (3), we get 



But G vanishes on 8 so that 


Fp = 


47r J dn 


dS 


( 1 ). 


By this means Dirichlet’s problem is reduced to the deter¬ 
mination of Green’s function. 

The case where the region is bounded by a closed surface S 
and the infinite sphere does not need special consideration 
if we impose like conditions as before to the way in which the 
functions vanish at infinity. 
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5*63] 

6*62. Green’s function for a sphere. Let O be the centre of a 
sphere S of radius a and P an internal point at a distance / from O, 
Let P' be the inverse of P, so that OP .OP' Let r, r' denote the 
distances of any point M from P and P'; and for a point M' on the 
sphere let PM' = and P'M' = r^'. 


Then by similar triangles 


so that 


r^' OM' a 
r,- OP -/’ 

a _ 1 


( 1 ). 



Now consider the function //= —alfr'. It is a potential function 
due to a mass ( — aIf) placed at P', and therefore it satisfies Laplace’s 
equation inside the sphere; and because of (1) it takes the form — 1/r 
on S. 

Hence .<2) 

is the required Green’s function for the sphere S and point P. 

The same formula can be shewn to hold good when the point P 
lies outside the sphere at a distance / from the centre, P' then being 
inside and the symbols having the same meanings. 

5*63. Reciprocal property of Green’s function. If P, Q 

are two points inside a region bounded 
by a surface S and 0 (P, Q) denotes 
the value at Q of Green's function for 
the point P and surface S, then 
0(P,Q)^0{Q,P). 

Apply Green’s theorem in the 
form 6*6 (1) to the region bounded 
externally by 8 and internally by 
two spheres Sg of small radii 
ca having their centres at P, Q. 



8-2 
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Put V and V = 0' — , where r, r' denote 

T T 

distances from P and Q and O, O' denote Green’s functions 
for P and S and for Q and 8 respectively. Then we have 

.( 1 ). 


By hypothesis H, H' are solutions of Laplace’s equation, 
and since P, Q have been excluded from the region of volume 

integration ^ = 0 and ^ = 0, so that = 0 and = 0 

and the right-hand side of (1) is zero. 

Again, by definition, 0 and G* vanish on iS, so that the 
first integral on the left is zero. 

In the second integral on the left we put so that 


which clearly -> 0 as e 0, since H is assumed to be finite; and 

and as € ^ 0 this reduces to — inOp, where Op' denotes the 
value at P of Green’s function for Q and 8, i.e. G (Q, P). 

In the same way the third integral in (1) may be shewn to 
reduce to 4:rrGQ or iTrO (P, Q); so that (1) reduces to 


0{P,Q)^0{Q,P) 


.( 2 ). 


The theorem is also true for a region bounded internally 
by a closed surface 8, and externally by the infinite sphere, 
assuming that H and H' are at least of order R-^ at a great 
distance R. 
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1. Use Green’s theorem to prove that, if F, F' are solutions of 
Laplace’s equation such that F = F' at all points of a closed surface 
then F = F' throughout the interior of S, 

2. Prove that, if F, F' are solutions of Laplace’s equation such that 
dVldn = dV'/dn at all points of a closed surface S, then F—F' is 
constant at all points inside S. 

3. Shew that, if F^, V 2 are the potentials of two equipotential 
surfaces which surround bodies of mass M which produce a gravita¬ 
tional field, the contribution of the space between these equipotential 
surfaces to the total energy of the field is (F^ — Fj). 

4. Prove that there cannot be a function F which satisfies Laplace’s 
equation inside a closed surface S and makes dV/dn= —kV at all 
points of S, where dn denotes the normal drawn towards S from the 
inside and A; is a positive number. 

5. A number of bodies are formed of gravitating matter. Prove 
that the resultant force, resolved parallel to the axis of x, due to the 

attraction on one of the bodies, of imiform density pj j(x-’X)NdS; 

where the integral is taken over the surface of the body, x is the abscissa 
of the centre of gravity of the body, and N is the normal attractive 
force at each point of the surface. [C. 1890] 


6. Prove that, if dm, dm' denote elements of two different distribu¬ 
tions of matter which have the same external level surfaces, and 
any function which satisfies Laplace’s equation, then 

j<j>dm: j<l>dm' = M: M', 

where M, M' denote the total masses of the two distributions. 


7. Prove that, if two distributions of matter have the same external 
level surfaces, they have the same centre of gravity and the same 
principal axes of inertia at the centre of gravity. 

[In Ex. 6 put ^ = 37, and <f>==yz,] 


8. Prove that, if F is the potential of a given distribution of matter, 
^ a surface surroimding it, r the distance of an internal point from 
an element dS, dn an element of the outward drawn normal and 6 the 
€Wigle between r and 5n, then 


[IdV 
j r dn 


dS 


f VcosB 
J r* 


dS, 


where the integration is over the whole surfaoe. 
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9. Shew that, if F is a function which with its first derivative is 
finite and continuous throughout a region bounded by a surface S 
and r denotes distance from a point P within this region, then 


47rFp 



laF 
r dn 


dS- 


I 


^ V^Vdv, 

T 


where dn denotes the outward drawn normal and the integrations are 
over the surface and through the volxime. 


10. Two gravitating solids M and M' have potentials F and F'; 
the inward normal force across any element dS of the surface of M 
is the sum of N and R due separately to M and M'; and the inward 
normal force across any element dS' of the surface of M' is the sum 
of N' and R\ duo separately to M and M'. Shew that the potential 
energy of their relative position differs by a constant from either 

^^(RV-NV')dS or -N'V)dS'. [M. T. 1897] 


11. If V^Fi = 0, V^Fq = 0 and Fj = Fq = F at a closed surface, wliile 
Fq , Fi and their differential coefficients are everywhere finite without 
and within the surface respectively, and Fq vanishes at infinity, find 
the potential in all space of the surface distribution a where 


and dn is outwards. 
Shew that 


4 

dn dn * 


/''sC") 


dS is equal to 


-/ 


r dn 


dS when r is 


measured from an external point, and J ~ when r is measured 

from an internal point. [C. 1908] 


12. Prove that, if <l> denotes any function which with its first 
derivative is continuous, its mean value over a sphere of radius a 
exceeds its value at the centre of the sphere by 

integrated through the volume of the sphere, where r denotes distance 
from the centre. [Minchin.] 

13. Assuming Gauss’s theorem and Green’s theorem and that a 
distribution of matter consisting of volume density p, surface density a 
and separate particles has a potential function F which is infinite at 
eetoh particle and satisfies the relations V®F = — 47rp wherever there 

is volume density, and across each surface layer 
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and is otherwise continuous, prove that the potential must be of the 
form 

v=i:'^+ 

r J r Jr’ 

where m denotes the mass of a particle. 

14. Prove the theorem of the equivalent layer [5*4] by applying 
Green’s theorem to the region boimded externally by the infinite 
sphere and internally by the equipotential surface S and a small sphere 
round the point Q, and excluding the singularities of the potential 
function. 


15. Prove that, if M be the mass of a given distribution of matter, 
U the potential of a closed equipotential surface entirely surrounding 
it, and E the exhaustion of potential energy in collecting the matter 
from a state of diffusion at infinity into its actual distribution, then 
iMU<E. [M. T. 1895] 


16. A function V vanishes at all points outside a surface S, Shew 

that a distribution of matter p= — --- V^F within S and a = } upon 

^ 4:Tr 4:7rdn ^ 

S will have the function V for its potential. Here dVjdn means the 

rate of increase of V in the direction approaching the surface from the 

inside. [C. 1931] 

17. Prove that, if Vq is the potential at points outside a surface S 
due to masses within the surface, and is the potential at points 
inside the surface due to masses outside it, the expression 


47r 






yeF, 


/ 4wj 

\dn 

dn ) 


-dS 

r 


has the value Fq at external points and the value F^ at internal points, 
where r denotes distance from the point and dn an element of the 
outward normal to dS. [M. T. 1910] 


18. F is the potential and R the attraction at any point due to 
a given distribution of matter none of which lies between two level 
surfaces F = F^, F = Fg. The space between these surfaces is filled 
with matter of density R'^f (F). Prove that its potential at any external 
fVi 

point P is — 47rUj ^ f {V) dV, where U denotes the potential at P due 

to that part of the original distribution which lies inside the inner of 
the two level surfaces. 


19. If F is the gravitational potential of a distribution of matter 
inside a closed surface S, shew that the value of F at a point P outside S 
is given in terms of its values on S by the formula 

1 f^^dG 
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where the differentiation is along the outward normal to /S', O vanishes 
on Sy and, if r is distance measured from P, (3^ — 1/r is a harmonic 
function on and outside S, Give a physical interpretation of O and 
dOldn. 

Shew that, if V is the potential of a finite distribution of matter on 
one side of an infinite plane, and P is a point on the other side, 

v,=lfv,dn, 

where dfl is the solid angle subtended at P by an element of the plane 
at which the value of F is . [M. T. 1934] 


20. Use 5*61, 5*62 to shew that, if F be a function whose value is 
prescribed over a sphere S of radius a, and P a point inside the sphere 
at a distance f from the centre, then a solution of Laplace’s equation 
inside jS is given by 




4t7ra 



dSy 


where r denotes the distance of dS from P. 


21. If and <^2 ^^e any scalar functions whoso derivatives exist 
and are continuous at all points of a volume V and its boundary surface 
Sy and if 

^/k = J^(grad,/,j.grad^i)dV, (j,k=l,2), 
prove that /jj gt /u®. 

If is a given harmonic function in F and if <^2 = </>i ori S, shew that 

■f 11 = -f 12 • 

Hence shew that for such a <f >2 the necessary and sufficient condition 
that Ja 2 should be a minimum is that <f >2 should be harmonic in F. 
What is the physical interpretation of this result? [M. T. 1935] 



Chapter VI 

HARMONIC FUNCTIONS 
6"1. Any solution of Laplace’s equation 

027 027 027 

-1-1-— Q 

dy^ dz^ 

which is homogeneous in x, y, z is called a harmonic function 
or a spherical harmonic. The degree of homogeneity is called 
the degree of the function. We shall be concerned with the 
case in which the degree is an integer. 

The study of harmonic functions is an important branch 
of mathematical analysis. We shall only concern ourselves 
with so much of the theory as is useful for our immediate 
purpose and be content with arguments based upon physical 
ideas rather than attempt a rigorous treatment of the subject. 

6-11. IfV is a harmonic function of degree n, then 
dx^ dy^ 02 * 

is a harmonic function of degree n—p — q — t. 

For if we differentiate the equation V^F = 0 p times with 
regard to x, q times with regard to y and t times with regard 
to z, we get 


V2{—— i--Fl = o. 
\dx'’^dy^d^ I 


6*12. Surface and solid harmonics. In polar co-ordinates 
Laplace’s equation may be written 


_8 

dr 



1 3 
^sin030 



1 dW 
^ sin® 6 3^® 


[1*68 (1)] 
.(!)• 
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If we substitute F==r”^„, where is independent of r, 

0 / 0T^\ 

we get = + so that (1) reduces to 

J . ^ i \ cy yy /yyv 

'W sin20+" (”■ +1) = 0 • • -(2), 


or, if cos0 = /x, 


_0 

dfjL 



1-/^2 0^2 


+ 7i(7H-l)/S^ = 0 ...(3). 


A solution of equation (2) is called a Laplace’s function 
or a surface harmonic of order n. Since n{n + l) remains 
unchanged when we write for n, there are two solu¬ 

tions of (1) of which S',, is a factor, namely and 
and these are called solid harmonics of degree n and — (n + 1) 
respectively. 


6*13. It follows that if ?7 is a harmonic function of degree n, 
then is also a harmonic function. For we may write 

U = , so that , which is a harmonic function. 

For example, it is obvious that xyz is a solution of Laplace’s 
equation of the third degree; therefore xyzjr'^ is also a solution. 

Similarly, if C7 is a harmonic function of degree — (n 4-1), 
then r 2 / 1+1 [7 jg harmonic function. For we may write 

JJ = so that r 2 n+i [7 ^ r'^8^ , which is harmonic. 

These theorems combined with 6*11 make it easy to build 
up a great variety of harmonic functions. 

Thus x^ or 1 and d or tan~^ (yjoc) are both harmonic functions 


11 y 

of degree zero. If we divide them by r, we get - and ~ tan“^ - 

T T X 

as harmonics of degree — 1. Differentiating these with regard 

X %l z z ?/ 

to x, y or 2 , we find that ^and tan“^- are harmonics 

^ r3 r® r® ^3 ^ 

of degree — 2, and so on. 

Again, if dh denotes an element of length whose direction 

0 0 0 0 

cosines are l,m,n, so that +in-;^ +n-^,it follows that 

oh ox oy oz 
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6 * 2 ] 


is a harmonic function and h may be called its axis. 
Similarly, if ••• number of arbitrary direc- 

. r • n .330 

tions, we can form a harmonic function ... 

vh^ ^^8 




6*2. Legendre functions.* Let Jf be a fixed point (/, g, h) 
and P a variable point (x, y, z), then if MP = R, we know that 
the inverse distance 1 /i^ is a solution of Laplace’s equation, 
for it is the potential at P due to a particle of unit mass at M. 
Suppose that M lies on a sphere of radius a with its centre 
at the origin O, so thsbt= and let OP — r. Then 
we have ^ ^ 

^ = :B = {(^fr + {y- gf + {z -W}i .. 

Since x^ + y^-{-z^ = r^, the expression for V can be expanded 
in powers of 1/r, by the extension of Taylor’s theorem to 
functions of several variables, thus 


.1/^3 3 ^3 


+ 


dy dz 
(~1)- 




ni 




0 0^0 

+ + 


\ni 

) 


.( 2 ), 


the series converging if r > a, i.e, if the point P is outside the 
sphere. But, if we write fi for the cosine of the angle MOP, 
we have 

y ^ _I_ 

(r^ —2arju-4-a2)l’ 

and this can be expanded in powers of rja when r < a, and in 
powers of ajr when r > a. Thus if we denote the coefficient of 
yn/^n+i Qp Qf ^ have 


and 


n--Pi+.. 

a 

fti \ 

,.+-P„ + ...} r<a . 

.(3) 

1+-P, + .. 
r ^ 

«" T, \ 

+ »•>“ •• 

....(4). 


♦ Adrien Marie Legendre (1752-1833). French mathematician of 
great distinction. 
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The P’s are clearly polynomials in p, for they are the 
coefficients of successive powers of A in the expansion of 

\l-2pM+h^)-i or 

1 + +... (2M - +.(5)- 


By picking out coefficients of successive powers of h we 
find that 


= P2 = i(3/x2-l), P3 = |(5/x3-3/a), 

and 


2n\ 

2^n ! n ! 



n(n-l) 

2{2n-l)^ 


+ 


n{n—l){n — 2)(n — Z) 
2A{2n-l)(2n-Z) ^ 



....( 6 ), 


the series ending with a term in ja or a constant according as 
n is odd or even. For many purposes it is sufficient to know 
that P^ is a polynomial in /a of the form 

+ +..(7), 

where , A^_ 2 , ... are constants. 

The function P^ is called the Legendre function or Legendre 
coefficient of order n. It is written P^ (/x) when it is desired 
to specify the variable. We observe that ^ = cos0, where 6 
is an angle measured from an axis OM, called the axis of the 
function. 

In particular, for 0 = 0 or tt we have /a= ± 1, and 

1 ^ 

P^ (1) = coefficient of h^ in -—^ = 1, 

...( 8 ). 

and P„(-l) =.= 

If we write /x = cos 6, then 
{1 - 2^ + h^)-i = (1 - ( 1 - 

= ^1 + ...) (.1 + + + ... J . 
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It follows that is of the form 

P^ = J5^cos W0 +jB^_ 2 C 08 (r^~ 2)0 + 

where Pn~ 2 »• • • positive. Hence the numerical value of 

does not exceed the number obtained by putting 0 = 0 in 
the above expression. 

Thus |P,(/.)1 ^P,(l);i.e.lP,(/x)I ^1 .(9). 


6*21. and are harmonic functions. Since 

jx — cos MOP = (fx '\-gy-\- hz)jar, from 6*2 (7) we have 


r"P„ = ~~ {fx + gy+hz)”- 

CL 

+ {fx+gy + hzY'-’^ {x^+y^+z^) + .( 1 ), 


expressing r^^P^ as a homogeneous rational integral function 
of the co-ordinates of P. 

Further, we observe that in 6*2, (2) and (4) represent the 
same expansion of F, so that by comparing the terms in 
we have 


„ (-ir/. 9 2 T. 


( 2 ). 


But 1 /r is a harmonic function, therefore by 6-11 every 
term on the right-hand side of ( 2 ) is a harmonic function of 
degree — (n + 1 ). Hence is a harmonic function, and 

by 6*13 r^P^ is also a harmonic function. 


6-22. Rodrigues’s formula for P„.* To prove that 

This formula may bo verified by expanding (~ 1)” by the binomial 
theorem and differentiating the result n times; or we may proceed 
as follows: 

Let l~/iw = (l-2^A + ^*)* .(1), 

so that w = /X + — 1). 

♦ Correep^ sur V^Jcole polytechnique, vol. m (1814-16), pp. 361-85. 
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By a theorem of Lagrange* u can be expanded in powers of h in the 
form TF ^ / 7 V 

(^*- 1 )"+.( 2 ). 

But by differentiating (1) we have 

^^ = (l-2(uft + fe2)-i = 2:A"P„ .(3). 

Then by differentiating (2) and comparing the coefficients of in (3) 
and in the result of differentiating (2), we obtain the required formula. 


6-23, The equation (/>t) = 0. Since the equation (/u,® — 1)’* = 0 has n 

roots equal to 1 and n roots equal to — 1, the equation ^ (/i.* ~ 1)" = 0 

has n—l roots equal to — 1, n— 1 roots equal to 1, and a root f be- 
(P 

tween ~ 1 and 1. Then 1)" = 0 has n —2 roots equal to — 1, 

n — 2 roots equal to 1, a root between — 1 and f and a root between 
f and 1. Continuing the process it follows that the n roots of the 
equation P„ (;x) = 0 are all real and lie between ± 1. 


6*24. Differential equation for P^. Since is a har¬ 
monic function and P^ is a function of 6 only, satisfies 
Laplace’s equation in the form 



and as in 6*12 we have 

or putting cos 0 = /it, 

.<*>■ 

This is commonly known as Legendre's equation. It can be 
solved in series in the form 6*2 (6). 


and 


6*26. Integral theorems. To prove that 

J_^-Pn(p)-Pm(p)<^p = 0. n^m. 


271-1-1 


♦ Whittaker and Watson, Modem Analysis, 7*22, 
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From 6*24 we have 






Therefore 


{w (w +1) - m (m +1)} PnPmdfJ- 


and, integrating by parts, this 




dfjL dfjL dfM dfJL 


Again, to find J (fx) dfi, we have 


SA-P,(p)j =(l-2/xfe+A2)-i. 


Integrate between limits — 1 to 1 and use the fact that the 
integrals of all products (n4^m) vanish, and we get 

_j.2n 


J jF ^2 -- coeff. of in J ^ 


• 2/iA + 


1, l-h 
-h'^^TTh 


2 

'2^Ti’ 


6«26. From the expression for (ft) in powers of ft, it is clear that 
by solving a set of linear equations any positive integral power of ft 
can be expressed in terms of a finite number of F’s in the form 

ft** ~ Uq P "4* Ci>2 P n—2 "t ®4 Pn—4 "t • • • 
and consequently, from 6*25, if n < m, 
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And similarly, if F (ju) is a polynomial of degree less than m, 

f^^F(^)F„(^)d^ = 0. 

6*261. If a polynomial in /x is expressed in terms of Legendre 
coefficients, say in the form 

the coefficient A„ can be found by multiplying both sides by PniH') 
and integrating between limits — 1 and 1. For then by 6*25 all the terms 
on the right except the term in vanish, and 


6*27. Recurrence formulae. These are a set of formulae con 
necting successive P’s or their derivatives. 

We begin with the foimula 

+ + .(1). 

Differentiation with regard to h gives 

(fx-h){l^2fih + h^)-'i = Znh^-^P^ (/x) .(2). 

Multiplying by (1 — 2/x/i + h^), we get 

(^ - ft) |l +1 ft"P„ (^t)| = (1 - 2ij.h + ft2) ^ nh”-^P„ i^,), 

then equating the coefficients of gives 

(n + 1) Pn+i - (2n + 1) iLtP„ 4 nP„_^ = 0.[A]. 

AgSjin, differentiate (1) with regard to /x and we have 

ft (1 - 2^ft + ft2)-« = I ft"-/-5, 

1 

and therefore, from (2), 

00 00 r/p 

h Z nh^-^P^ ^U-h)Zh^ . 

1 1 

Whence, by equating coefficients of 
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and eliminating dPJdfj. from this result and [B] gives 

.[C]. 

Subtracting [B] from [C] gives 

. 

In [D] change n + 1 into n and subtract /x times [B] and we get 

= .[E]. 

Integrating [C] gives 

fp„Md^= .[*’J- 

Lastly, by adding successive formulae of the form [C] we get 
dP 

= (2n- 1) -f (2n - 5) P,_3 + (2n- 9) + ..., 

the series ending with 3Pi if n is even, and with Pq or 1 if n is odd 

.[C>]* 


6*271. Example. Use the theorems of 6*25 to prove 6*27 [G]. 

Since P„ is a polynomial of the nth degree in fi, therefore dPJdfjL 
is a polynomial of the (n — l)th degree and by solving linear equations 
is expressible in the form 

dP n ry , P , 

— ®n~l ^n-1 "r «n-8 ^ n-3 -T • • • • 

To find any coefficient, say multiply by P^ and integrate 
between limits — 1 and 1; then 


since by 6*25 every other term vanishes. 


2mH-1 * 


by 6*25 


( 1 ). 


But since m<n, dP^jd^ is of the form c„j_iP„j_i+ c„»_ 3 P,„_a + 
where the suffixes are all less than n. Therefore 

. 

By adding (1) and (2) we get 


2m + i 



since n — m is odd, 6*2 (8). 


RNA 


9 
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Therefore a„ = 2m + I, and 

-^” = (2n-l)P„_, + (2»»-6)P„_, + .... 

6*28. Definite integral forms for . There are two ex¬ 
pressions for as a definite integral, due to Laplace,* which 
may be obtained as follows: we have 

dcj} IT 

Jo a-h6cos</»~ — 

so that, by putting a—l—yJi and 6 = ± A — 1), we get 
1 1 
V(1 - ~ 

_ 1 P _ d(l> 

ttJo I — fjih ±h — 1) coB(f>* 

Expanding both expressions in powers of h, we get 

'Lh^P^ = - r 2A^{/Lt+ 1) cos <f)}^d<l>y 

0 '^■Jo 0 

and, by comparing the coefficients of 

-P« = -f {m+A/( p*-l-)co8.;6}"<i^6 .(1). 

1) is imaginary for values of /x between — 1 and 1, 
but if the binomial expression is expanded the integrals of all 
odd’^powers of cos</> vanish and a real result is obtained. 

By writing 

cos <f> — {fl cos i/j ± l))/{jLt± l)cos^}, 

the result may be transformed to 

Pn = -{ ± -1) cos dAjj .(2). 

6*29. After this lengthy discussion of the properties of the 
Legendre polynomials we revert to their introduction in 6*2 
as a special form of surface harmonic, namely a solution of 

♦ Micanique Celeste, liv. xi, ch, ii, p. 33 [1823h 
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6 ' 3 ] 

6*12 (2) independent of <j^. But we must point out that, if 
in 6*2 we take M as the point (a, d'y(f>*) and P as (r, 0,^), then 

F = {r® + — 2ar [cos 6 cos 6* 4- sin 6 sin 6' cos {<f> — 

.( 1 ) 

=:(r2 4-a2~2arA)-* (2), 

where X^iifi + \/(l — /x^) \/(l 
fi = cos 6, fji = cos O' 

This as before is expansible in powers of rja or a/r, and it 
may be shewn, as in 6*22, that the coefficient of or of 

anjrn+i is where A is given by (3). 

Expressed in this way the Legendre function (A) is 

a symmetrical function of 0, <}> and 0', (f>\ i.e. of the angular 
co-ordinates of P and M, and it is only when the polar axis 
of co-ordinates coincides with OM or OP that we have X — ^jl 
or A == ^' and the Legendre function takes the form of a poly¬ 
nomial independent of </> and </>'. 

6*3, Applications of the Legendre functions. Since the 
Legendre function (/x) is a function of fx or cos 0 only and 
r^Pn and r~^~^P^ are solutions of Laplace’s equation, it is 
evident that these functions may prove to be of use in connec¬ 
tion with the potentials of bodies with an axis of symmetry. 
In this connection there is a theorem of Legendre: If the 
potential of a body with an axis of symmetry is known at points 
on a finite length of the axis outside the body it can be calculated 
in terms of the Legendre functions at all points that can ber cached 
from this length of the axis without traversing matter. It is to be 
understood that there is symmetry both as regards the density 
and the shape of the body. 

Taking the axis of symmetry as the axis Oz, let Q be a point 
on the axis outside the body at a distance z from the origin. 
And suppose that the potential at Q is expressible in a series 
in the form 

P-<..+^»+a.2+^ + ...+a.*"+i +.(1) 



9"2 
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and that this series converges for the pertaining to a certain 

length of the axis. 

Also let F denote the potential at a point P, where OP = r 
and QOP = 6; and let 

F' = ao+^ + |air+~jjPi + ...+^a„r«+^jp„ + ... 

.( 2 ). 

Then F' is a harmonic function [6*21], so that F — F' is also 
a harmonic function. Now let P move up to then F becomes 
U and F' becomes the series on the right of (1) (since at Q 
T — z and d = 0), so that at Q, and other points on Oz for which 
(1) converges, F—F' vanishes. Therefore by 4*3 (v) the 
potential function F — F' vanishes at all points that can be 
reached in passing from Q and other neighbouring points on 
the axis without passing through the attracting matter or 
through a region in which (1) ceases to converge. We conclude 
therefore that the potential at any such point P is given by 

F = ao+^ + (air+^jPi + ... +(a„r"+^^P„ + ... 

.(3), 

and this series can clearly be written down, whenever (1) 
has been calculated, by writing r for z and introducing the 
appropriate P in each term. 

It is clear that (3) is absolutely convergent for all values 
of r which make (1) absolutely convergent when r is written 
for z, since | P^ | ^ 1. 


6*31. Potential of a circular ring. The potential of a 
circular ring of radius a and line density m at a point on its 
axis at a distance z from its centre is 27raml{a^+z^)^, This can 
be expanded in powers of zja or ajz according as z is less or 
greater than a, so that on the axis the potential is given by 


either 


rr ^ \ 


C4 = 27rm 


1^ 1.3a5 

2z^^ 2Az^ 



or 


z>a. 
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Therefore the potential at any point (r, 9) is either 


^\r 2 : 


:sP2 + 


1.3r‘ 1 

i 


, r<a, 

2 .4 0 * * J 


1.3a5 ^ 

[ 


, r>a. 

2.4r® * j 



6*32. Potential of a circular disc. With the same nota¬ 
tion the potential at a point on the positive axis of z is, by 
2*51, 27Tm{\/{z^-{-a^) — z}, where m is the mass of unit area. 
This can be expanded in the forms 

^ { 1^2 1.12^4 1.1.32:6 ) 

U^=^2TTm\a-z + ~ -+ ••• > 2;<a, 

[ 2a 2Aa^ 2.4.6a^ j 

, ^ (la 2 l.la4 1.1.3^6 ] 

and c 4 = 27rm{-- tt-t+ , z>a, 

* \2z 2.4z® 2.4.6Z® /’ 

Therefore the potential at (r, 0) is given by 


and 

27rm 


, Ir^ 


1 

.1 

„ 1. 

1.3 r« „ 1 




2 

74: 


T> 1 

4.'6a® * ’"J 

1, r <a 







....(1), 

l.la* 


1. 

1 

■ 3a6 „ 


....(2). 

2Aj^ 

^2 + 

Y, 

T 

_P _ 

... , r>a .. 


If r >a, the expression ( 2 ) is vaUd for since we can 

move from the point (r, 0 ) to the point (r, 0) keeping r>a 
without crossing the disc. 

If r<a, and 9 > we cannot move from (r, 0 ) to (r, 9) 
keeping r<a, without crossing the disc. Thus ( 1 ) only holds if 
0^9 ^ ^TT. 

If ^TT < 0 ^ TT, then, by symmetry, (r, 0 ) = 1 ^ (r, tt — 0). 

Now Pj, (tt — 0 )= -Fj^(9), and P 4 , ... are unaltered if 0 
is replaced by tt —9, 

Hence when r<a and ^tt< 0 ^ tt 

( 17*2 1 1 ) 

Thus, as in 2*51, the form of the potential is discontinuous 
in crossing the disc. 
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6*33. Example. Find, in the jorm of an infinite series, the potential 
due to a flai circular annulus (e,g. Saturn''s ring) at points moderately 
near the centre of the ring. 

If a small hut massive planet occupies the centre of the ring, shew that 
the maximum deflection of the plumb line on its surface {on account of the 
attraction of the ring) will occur in latitude 46°, approximately, 

[M. T. 1916] 

Let a, h be the inner and outer radii of the annulus. Its potential 
is the difference between the potentials of discs of radii a and 6. Taking 
the axis of the annulus as axis of z with the centre as origin, the potential 
at the point ( 0 ,0,2) is 

U = 27rm {\/ {b^ + z^) — \/{a^ + z^)], 

where m is the mass of unit area of the aimulus. 

For z<a, this gives 

Therefore, at points near enough to the centre for higher powers of 
rfa to be neglected, the potential is 

If the planet be of radius c, the deflecting force on the plumb line 
at a point on its siu-face at an angular distance 6 from Oz is dV/cdd, 
and if c is small compared with a, this is 

2rrmh c^ dP^ 27Tma c^ dl\_ /I \\dP 2 

W "de c“' 2o2 He ~ " a) He ' 

But P 2 = ^ (3 cos^ 1), so that the deflecting force is 
— f Trmc Q — ^ sin 26, 
and has its greatest value when ^ = 46°. 

6*4. Harmonic functions in general. We now turn back 
to the subject of harmonic functions in general. From 6*1 
will denote a harmonic function of degree n if it is homo¬ 
geneous of degree n in x, y, z and satisfies Laplace’s equation. 
If it is a rational integral function of degree n in x, y, z, it 
contains coefficients; but since V2jEir,^ = 0, and 

this relation is of degree n — 2, there are — relations 
between the coefficients. Therefore -¥\) {n -\-2) -- \n{n--1), 
or 2n-{-\, of the coefficients are independent, so that there 
are 2n-{-\ independent rational integral harmonic functions 
of degree n. 
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6*41. Integral theorem for surface harmonics. To prove 

that js^S^doj^Oy where 8^ are surface harmonics of 

different orders^ dw is an elementary solid angle^ and the integral 
extends to the whole surface of the unit sphere. 

Apply Green’s theorem 5*1 (1) to the region bounded by 
a sphere of radius a with its centre at the origin, putting 
V = r'^8^ and F' = , so that F, F' are harmonic. It follows 


^ 9F' 0F 
J dr dr 


dS==0. 


when the integration is over the sphere r == a, so that dS = a^da). 
The relation reduces to 

(n — m) 8^ 8^^ dco = 0, 

and, since w =t= m, therefore 

= 0 .(!)• 


Particular cases: (i) Since is also a surface harmonic, 
therefore c 


J«^F„dco = 0 .(2). 


(ii) Since therefore 




(iii) Since P^^ may be taken as special cases of 8^, 8^, 
and we may write dw = sin dd0d<f> = — d(id(f>, the integral over 

rl |*27r 

the sphere in this case becomes J J P^Pj^dpd<f>, or 

P^PndlL, so that (1) reduces in this case to the theorem 

of6-25,J P^P^diL = 0,m^n. 

6*42. Tesseral heirmonics. Reverting to the equation 
8’12 (2) for surface harmonics, viz. 

d^Sn .n^Sn 1 r> /i\ 

-^ + cot0-^+^in2^ g +«(TC+l)(Sf„-0 ...(1), 
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we proceed to examine solutions of the form 8^ = 0®, where 0 
is a function of 8 alone and O a function of alone. On multi¬ 
plying the equation by sin^0/0<I>, it becomes 


sin2 0(020 


0 (002 - 

We must therefore have 


-f-cot0-5^-f72,(?i4-l)0 = 0 ...(2). 


sin2 0(020 .30 ^ 


The latter equation has single-valued solutions when 

c—p, viz, ^ = Cp cos p(l>-\-Dj^ sin p<f> .(3); 

and © is then given in terms of 0 by 

or, in terms of p, by 


f 




0 = 0.(4). 


ror^ = 0 this is Legendre’s equation with a solution P^(/x), 
where ^ r dP \ 

. 

Differentiate (5) p times with regard to p and multiply the 
result by (1 This gives 

dp+2p 

dP P 

+ {n{n+l)-p(p+l)}{l-^^)ip^=0 .(6). 

dPP 

Let Z = (l-u2p^^; 

^ ' dpP ’ 

then 1=(1 
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Therefore 


A 

dyi. 


<‘-•'11 


/7P+2 P /fP+1 P 


d/iP+* 




P dPP 

-p{l-y?r +i>V (1 - - (7). 

By subtracting (6) from the right-hand side of (7), we get 


A 

dfi 




dz 

dfjL 


= I - »(n + 1 )+(i? + 1 ) - i? + (1 - 

p2 


= i)- 


i-p* 


so that ; 2 : is a solution of equation (4). It is convenient to 
denote solutions of this form by and, using Rodrigues’s 
formula for , write 


^ '2^'n\p\ dJI^ 


(/x2-l)-, 


^9 = 0, 1, ...n...(8). 


It is clear that is a rational integral function of sin 9 and 
cos 9, and that it vanishes if n+p> 2n, i.e. if 79 > 7 i. Also since 
00 is a surface harmonic, therefore 


Y 0 Y ico8^ X p^9^7)rk X 

■^n f -^n sin r> *** Bui-rr> •** sin 

are 2n-f 1 surface harmonics. They are called Tesseral har¬ 
monics. That they are linearly independent follows from the 
fact that if we had a relation 


AqX^^^ -h AiX^^ co8<f> + ... +A^Xn^coBn(/} 

4-j9iX„^sin^+ ... -h-B^X„^sinn^ = 0, 

it would be sufficient to multiply by and integrate from 
0to27r, toshewthat^^ — O* J5p = 0forp=l, 2, ...n; and then 
.40 = 0. 



138 


HARMONIC FUNCTIONS ' 


[643- 

6*43. expressed in tesseral harmonics. From 64, every 
surface harmonic of order n can be expressed as a linear com¬ 
bination of 2n^l linearly independent surface harmonics, 
so that 

S^ — AQX^ + A^X^GOS<i)-\- ... +-4^X,j,^cosn<ji 

+ ... +B^X,^^min<f> .( 1 ), 

where the Ji’s and J5’s are constants. 

Also since = therefore is expressible in terms of 
harmonics of the form r^X^^ slnP^y 
these are polynomials in x, y, z. 

Thus 642 (8) may be written 

“ p\ .' ’’ 

and taking for its expansion in powers of /x [ 6*2 (6)] we see 
that X^P is the sum of terms of the form sin^ 6 (cos ; and 

sinpcf) sin^ 6 (cos (r cos sin^ 6 sin^^. 

Here,the first two factors are polynomials in x, z\ and for 
the last factor, if ^ is odd say and =25+1, 

sin^</> = sin</>{sum of terms of the form (8in</>)^"^”*}, 

and 

sin^ 6 sin = r sin d sin <f> {sum of terms of the form 

(r sin 0)2"* (r sin d sin 

= y{mm of terms of the form (a:2 + y 2 )m^ 2 «- 2 mj 

Similarly if p is even; and thus f^X^ is expressed as 
a polynomial in x, y, 2 , and is expressed in terms of 2w + 1 
linearly independent polynomials. 


6*44. The integral 


^S^P^do), To prove that, if 8^ be any 

surface harmonic of order n of the angular co-ordinates 0, cf) and 
P^ the Legendre function of the same order having (0',^') as its 

r 477 

axis, then xS^P^do) — - - where the integration is over 

j Zn +1 

the whole surface of the unit sphere and 8^' is the same function 
of <f)' that 8^ is of 6, <f>. 
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Take the axis of as the axis from which d is measured 
so that fjL = cos 9. Then 

do) = sin dddd(f> = — dyid<f>, 

and 

/• rl r2rr ( n 

J = J J^ 1^0 + S {ApX^PC08p(f> 

+ BpX^P 8in^)| Pndpd4, 

where are constants. 

P2it P27r 

But cosp<f)d<f) = 0, sin^d^ = 0 and = P^, so that 
JO jo 

where i® constant. 

But from 6*42 (4), when = 1, vanishes save when p = 0. 

Hence, when /x= 1, 

But, when p=l, Pn~^^ so that the constant is the 
value of 8^ at the point on the sphere for which 9 is zero, 
i.e. on the axis Hence Aq = 8^' and 

. 


6*45. Expansion of a function in surface harmonics. A 

polynomial in cos 9, sin 9 cos (f), sin 9 sin ^ can be expressed in 
terms of surface harmonics. For this purpose we use the 
following theorem: If is a solid harmonic of degree n, then 
V2 ) = m (27^ + m +1) , 

where r^ — x^+y^ + z^. 


We have ^ = mr^-^xH^^ + 
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By adding the similar equations in y and 2 , we get 


[6-4S- 


a;^ + y-^+2-^j. 
But since is homogeneous and of degree n 


0Br afT dH^ „ 


therefore ) = m (2ii -f m H- 1) .(1). 

Now let/„ {x, y, z) be a homogeneous polynomial of degree n 
in X, y, z, and suppose that 


/» isc, y,z) = H^+r^n-i + +.(2). 

where are solid harmonics of the degrees in¬ 
dicated by the suffixes and the series ends with or 

r'^H^ according as n is odd or even. By using the result (1) we 
can shew how to determine the Af’s. 

Thus by applying (1) successively, we get 


Vy, = 2 (2n - 1) + 4 (2n - 3) 

+ 6(27i-~5)r^ir^^e-f... 

Vy^ = 2.4(2n-3)(2n-5)jy,_, 

-[-4.6(271—5)(2w—7) -f... 


(n— 1) (7z,-f2) (77 — 3)71 ... 2.5i?'i (77 odd) 

or V"/^ = 77 (77 + 1 ) (77 - 2) (77 — 1 )... 2 . (77 even). 


From the last of equations (3) we can find or according 
as 77 is odd or even, then from the preceding equation we can 
find or and so on until the up to determined 

and then is given by (2). 

By changing into polar co-ordinates, writing rco&O, 
rsinScos^, rain6sin^ for x, y, z, and then dividing by r^, 
it follows that any pol 3 momial in cos0, sin 0 cos sin 0 sin 
can in this way be expressed in a series of surface harmonics. 
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6 *6. Surface density in terms of surface harmonics. If 

as in 6*2 we find the potential at P due to a number of particles 
situated on the sphere of radius a, we get results of the form 

00 

Fx = S—»•<« .(1) 

and = .(2). 

where denotes the sum of a finite number of surface 
harmonics (one for each particle) and therefore itself a surface 
harmonic. 

We infer that the same result would be true for any arbitrary 
distribution of surface density on the sphere, since every 
particle of the matter gives rise to potential functions of the 
form stated. 

Conversely, we can assume (1) and (2) to give the potential 
of a certain distribution of matter and proceed to find it. By 
hypothesis harmonic, so that = 0 and = the 
matter therefore resides on the surface of the sphere, and its 
surface density u is given by 



It follows that if we accept the physical argument that an 
arbitrary distribution of surface density on the surface of 
a sphere produces the same kind of field of potential as an 
aggregate of particles distributed over the sphere, viz. the 
potential given by (1) and (2), then the arbitrary surface 
density is expressible in the form (3). This implies that an 
arbitrary function of the two variables 0, (f> on the surface of 
a sphere is expressible in a series of surface harmonics. The 
validity of this development in series has been discussed by 
many writers; a bibliography is given in Encychpedie des 
Sciences Mathimatiques, * 

♦ Tome n, vol. v, p. 176. F. also H. Poincar^, Figures d^equilibref 
p. 62 [1902]. 
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6*61. Reverting to the case in which the field is due to a single 
unit particle at a point M on the sphere r = a (6*2), it is of interest 
to verify that the formula 6*5 (3) for the surface density gives zero 
surface density except at the point M. 

Taking the values of V from 6*2 (3) and (4), the formula 6*5 (3) is 
easily seen to give 

47ra2cT = S(2n-f 1) 

0 


and it is convenient to take for the formula 6-21 (2) putting r = a, 
so that 


47ra^cr 


-00 

= 2 
_0 


(-l)M2n-f 1) 


-(/ 


d 

dx 






where/, h are the co-ordinates of M. 
This is equivalent to 


oo 

4t:7Tacr = 2 
0 


[ 


-2n(-l )»-^ 
n! 




/S . s ,I a 


1 

r 


n! 


^ dx'^^ Sy'^^dz 



But 


0 



7 r 0 n! 



) 


”1 

r 


Therefore 

4’“'[(-K4+4"'‘5)+‘}((r-/)‘+...)iL 



6*6, Potential of a thin spherical shell. In accordance 
with 6*5 we assume that any distribution of surface density 
on the surface of a sphere may be represented by a series of 
surface harmonics. We may therefore represent the surface 
density at any point M {a, by 

o=^ S^-]r 


(1). 
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If P is the point (r, 6 , <f>), the potential at Pdue to an element 
adS of the shell is (as in 6*2 and 6*29) given by 




and F. = - 


l+-Pi+^P2 + ... + —P„ + ... , r>( 


And the potential at P due to the whole sphere is got by 
substituting the series (1) for a and integrating (2) over the 
sphere, and by writing a^dco for dS and using the theorem 
6-44, wo find that 


>2 = ----- \ + 5 -'^! + s :: 2'^2 +••• + ? 


1 o" 
2 n +1 


r>a 


where Sq\ 82 , etc. denote the values of the surface har¬ 
monics ^ 0 , Sly 82 , etc. at the point on the sphere where the 
radius through P intersects it. 

It is easy to verify that the formula of 3*7 

' -Hf-fL. 

leads to the expression (1) for the surface density. 

6*61. ExEUtnples. (i) The density of a thin spherical shell of radius 
a is given by a=^Xz^y where z denotes distance from a diametral plane. 
Find the potential. 

We must first express the density in terms of harmonics. We have 
(T = = Aa* cos^ 6 

= 4AaMPo+2P,) .(1), 

since Po=l and 6—1). We can now either write down 

the potential from 6*6 (3), or we may assume a series of harmonics 
for V and find the coefficients by using 6*6 (4). Thus, let 

Fi = CoPo + c,^F,, r<o| 


Fg = C 0 ^ Pq 4 - Cj “ Pg, r > o 
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Then substitute for Fj, Fj and a in the relation 


\dT dr Jr^a 


and we get 


P 4- 


5co 


Pa = i7rAaMn + 2Pa) 


[6*61 


(3). 


Since this relation must hold good at all points on the sphere, 
i.e. for all values of 6, the coefficients of Pq and Pg must be the same 
on both sides, so that 

Cq = J nXa^ and Cg = |. fi ttA a® 
and ^1 = ^TrAa® -f|^P 2 ^» r<a 

F. = i|.Aa»r>a 

We remark that in adopting this latter method we only insert in 
the series (2) terms to correspond to the harmonics that occur in (1). 
If any additional terms were inserted, the relation (3) would serve to 
shew that their coefficients were all zero. 

Further, the total meiss M of the shell is given by 

M= f Aa®cos®0.27Ta®sin^d0 = fj7rAa*, 

J 0 

and since 2 r®P 2 = r® (3 cos® ^ — 1) = 2^® — a;® — y®. 



therefore the expressions for the potential may be written 




r<a 

r>a 


( 6 ). 


(ii) The density of a solid sphere of radius a is given by p:=:Xxyz. 
Find the potential at an external point. 

Consider a thin shell of radius a' and thickness da'. We may regard 
it as a layer of surface density 

a = pda'= Xxyzda' .(1). 

Now xyz is a harmonic function of the third degree, so that by analogy 
from 6-6 this and its associated harmonic xyzfr"^ [6«13] are the only 
harmonics which can occur in the expressions for the potential. We 
therefore assume that for the potentials of the shell 


Vi^Cxyz, r<a' 


and F,= Co'»^, r>a', 

arranging the constants so that Fj = Fj when r = a', Fj is finite at the 
origin and F 2 0 as r oo. 
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Then substituting in 


(SJ^ 


^ )r-a' 


= 4iro. 


we get 




SO that G = i^nXa'da'. 

Hence the external potential of this shell is iirXa'^da' . And by 
integrating with respect to a' between limits 0 and a, we get for the 
potential of the sphere rriTrAa®—?. 

(hi) A mass M is distributed over a spherical surface of radium a so 
that the surface density at any point is proportional to yV + + x^y^t 

where the origin is at the centre of the sphere. Find the potential inside 
and outside the sphere. 

It is necessary in the first place to express the surface density in 
a series of harmonic functions. Proceeding as in 6*45, w© put 


yV + 4- x^y^-fix, y, 2) = H4 + 4- .(1), 

where , Hq are harmonics of degree 4, 2, 0. 

Then V2/=2.7H2 4-4.5r2//o .(2) 

and VV‘=4.5.2.3Ho (3). 

But, by differentiation, 

V2/=4(a;2 + 2/^4-32) 
and VV=24. 

Therefore, from (3), (2) and (1), 


and 


^2 = 0 

H 4 = yV 4 - z^x^ 4 - x'^y^ — Jr*. 


We may therefore express the surface density in the form 
(T = A {y^z^ 4- 2 V 4- x'^y^ — Jr*) 4- J Aa*, 

and assume expressions for the potential containing the same har¬ 


monics, VIZ. 


Vx=A + B {yH^ 4- z^x^ 4- x^y^ — Jr*), 


r<a 


and 


F, = A “ 4- B {y^z^ 4- zV 4- x^y^ — Jr*), r > a. 




r r'* 

Then, by substituting in the relation 


[sr- 


RN A 
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A dJS 

we get - + —- (3/®z® + -h x‘^^ — Ja^) 

= 4trA (yh^ 4- 4- a;*!/® — Ja*) 4- firAa^, 

so that A = f ttAo® and J5 = ^rrAa. 

But ^ — j integrated over the sphere, 

= J jxa^dSf since the integral of a harmonic is zero, 
= f 7rAa«. 

Hence it follows that 

F, = M + g® , (yV + + a*2/* - Jr*)! 


[ 6 - 61 - 


and 


F, = Af { J +1 (2/«z* + z»x* + x^y^ - Ir*)}. 


6*7* Nearly spherical bodies. Sux^pose that the strata of 
equal density in a body are nearly spherical, so that the 
equation of a surface of equal density is 

r=:a(l+SC'„/SfJ, 

where the coefficients are so small that their squares and 
products can be neglected. 

The volume contained by such a surface is [1-42], 

which f 

= iJ(a^4-3a3SO,6fJda> 

= since ^S^doj = 0, [6*41], 

== vol. of a sphere of radius a. 


We note as special cases: 

(i) The surface r=a -j- eP ^, 

or r=a4-€Cos^. 

When € is small enough, this is a sphere of 
radius a with its centre <7 at a small distance e 
from the origin. 

(ii) The surface r = a4-€P2, 
r® = a® + 2 a€P2 • 



or 
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To the first order in powers of «, this is the same as 

r>=a2 + 2-€P,, 
a * 

or r2 = a*4-—€(3 cos* 0—1). 

a ' ' 

In Cartesian co-ordinates this can be reduced to 






SO that the surface is a spheroid of semi-axes a — Je, o -f c. 

We also notice that the normal to such a surface makes with the 
radius vector r an angle x which is at most of order €. So that for 
differentiation along the normal we have 


dn~' dr ^ ^rdd^r8mdd<l>* 

where A = cos x = 1 to the first order, while /x, v are of higher order than 
the first in e. Consequently, to the first order in €, we have 

eVjdn = dVldr, 


6*71 • Potential of a nearly spherical body. Consider a 
body stratified as in 6*7. The surface 


r=a(l+SC„5„) .(1) 


may be regarded as the boundary of a solid sphere of radius a 
and unit density surrounded by a layer of surface density 

alLCA- 

The potentials inside and outside the sphere are, from 6 * 6 , 


U, 


C„ 


; = f^(3a*>-r=>) + 4^S^— r<a 


4 rrCl^ C 

and ?72 = o — + ^ ^n', 

^ 3 r 2?i + lr"+t "" 


r>a 


...( 2 ). 


If we regard a as variable, the differentials of will 

give the potentials inside and outside the shell whose para¬ 
meters are a and a^da. The density p of the corresponding 
stratum of the given body is by hjrpothesis a function of a, 
and if c, c' are the values of a on the inner and outer boundaries 


10-2 
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[6-71~ 


of the body, the potentials inside and outside the shell are 
given by 


and 


r, = 4,rj%[a 


da 4-1 a'^~A 


2 . 

da 271 -f 1 r'^ 


I daj 


.( 3 ). 


The potential in the substance of the body, on the surface 
of parameter b say, will be found by adding together the first 
integral between limits b and c' and the second taken from 
c to 6. 


6*72. Example. The form of a nearly spherical uniform gravitating 
solid of mass M is given by 

r = a(l + €Zmn), 


where I, m, n are the direction cosines of r. Prove that neglecting the 
exterior potential is 


(1 3 


aHmn] 

~^r 


If the solid is surrounded by a thin layer of water just sufficient every¬ 
where to cover its surface, shew that, neglecting the self-gravitation of the 

water, the maximum depth of the water is * [C. 1931] 

Since c is small we may regard the solid «as formed of a sphere of 
radius a and imiform density p surrounded by a thin layer of surface 
density a = a€plmn; and by [6*7] ^ 7 ra^p = M, the whole mass. 

Also, since xyz or rHmn is a solution of Laplace’s equation, therefore 
[6*13] xyzjr"^ or Imnjr* is also a solution. 

The potentials inside and outside the sphere due to the surface 
layer may be taken to be of the form 

^3 ^4 

Fi=^ and Vo = A ~:lmn, 

^ ^ y4 » 


so that Fj is finite at the centre, Fg vanishes at infinity and Fi=: Fg 
at the surface. Then from 3*7 and 6*7, 

= 47rya= 4:7rya€plmn, 

Therefore 7 A = 4:7rya^€p = 3yAf e /a, 

and the external potential due to the whole body is 


{ev,_8V,\ 

\ 8r dr } 


( 1 ). 
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The free surface of the water must be an equipotential surface 
because in equilibrium the surfaces of equal pressure are at right angles 
to the resultant force. 

Hence the free surface is given by an equation 

~+|e^inOT=C7 .(2), 

where C is a constant to be determined from the fact that the water 
only just covers the solid. The surface of the water will touch the 
solid on its protuberancas, i.e. where the surface r = a(l4-€^mn) has 
its ma^^imum radii; and the water will be deepest over the hollows 
on the surface of the solid, i.e. where r = a (1 + €lmn) has its minima. 
It is easy to shew that Imn has maxima and minima values subject 

to Z® -f + n® = 1, where Z = ± , m = ± , n = ± , maxima 

VO Y o Y o 

when the signs are chosen to make Imn positive, and minima when 
negative. 

Hence O is to be found by substituting a -f 

for r; so that 

C-~(l _^ ^ i ^ \ 

a\ 3V3;^7v'3a"'a\ 21VV* 

and the free surface is given by 

= .(3). 


A minimum radius of the body is obtained by putting Imn = — « 


so that the minimum radius is r = a 


3^3)’ 


and the radius of the 


free surface in this position is given by 
1 1 €a8_l/ 

To the first power of €, this gives 


r a\ 2Iv'3^7v'3/ a\ 21^3/ 


So the greatest depth of the water 

^ 2173) " 373) 

=8*0/2173. 
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6*8. Clairaut’s Theorem.* Variation of gravity. As the 

result of the measurement of meridian arcs there is reason 
to suppose that the earth is a slightly oblate spheroid. If 
a, c are the major and minor semi-axes of such a spheroid 
and c = a(l —e), € is called the ellipticity of the spheroid. 
Neglecting the square of the ellipticity the equation of the 
surface is o . o o 

+ y ^ 1 

a* ''■a*(l-2e)"" 

or + -f-2e) = a2 

or (1 -h 2€ cos^ d) = 

or r = (x(l — ccos^^) .(1). 


For the earth, € is about 1/300. 

The earth rotates with a small angular velocity a> = 27r/24.60^. 
Bodies regarded as at rest relative to the earth are therefore 
subjected to a ‘centrifugal force’ as well as to the earth’s 
attraction; so that if V denotes the potential of the earth’s 
mass and X, F, Z denote the whole force on a unit particle 
on the surface, 


ox 






dV 

dz 


.( 2 ). 


As by far the greater part of the earth’s surface is fluid, 
it is reasonable to argue that the surface must be one of con¬ 
stant pressure, and from hydrostatical considerations this 


requires that Xdx+Ydy +Zdz = Q .(3), 

SO that F-f {x^-\-y^) — C .(4) 

over the surface, where (7 is a constant. 


If E denotes the mass of the earth, its potential at a great 
distance is j&/r, as a first approximation. We have to correct 
this because of the term ecos^^ in (1), and it will suffice for 
our purpose to add a solution of Laplace’s equation containing 
cos^0, i.e. a term APg/r®, where A is a small constant of 
order €. 


♦ Alexis Claude Clairaut (1713-1765). French mathematician. 
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-p Ap 

Thus we take V — —I—~ .(5) 

T 


and determine A so that (4) is satisfied by (5) when r is given 
by (1); i.e. 

-{l+€ cos2 e) 4- ™ (1 - cos2 e)^C 

Cl ji 

.( 6 ). 


r = a is a sufficient approximation in the second and third 
terms since A and are of order e. 

Since ( 6 ) is to be satisfied for all values of 6, by equating 
to zero the coefficient of cos^ 0 , we get 


A 




and the terms independent of cos® 6 serve to determine 


•(7), 

C. 


Hence F = (cos®0-i) .( 8 ). 

The radius from the earth’s centre to the point on its surface 
makes an angle d with the axis Oz, The force of gravity g 
at the earth’s surface is the resultant of the attraction and the 
centrifugal force and if it makes an angle v with the radius 
vector the radial force is — cos v. But v is so small that its 
square is negligible, so we have 

0 

g= — 0 ^ (F 4 - sin^ 0 ) at the surface, 


= ™ 4- ”j (cos^ ^ — J) coh sin^ 0, 

where r = a (1 — € cos^ 0 ) and is of order c. 

This gives 

S'^ . 

If 0 denotes the value of g at the equator 0 = ^tt, and m 
denotes the ratio of the centrifugal force to gravity at the 
equator so that = mO, we have 

(l + |m)^=(1+6)1 


( 10 ), 
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SO that jE/=(l+fm-€)C?a2.(H) 

and (9) becomes g^ = G^{H-(fm“€)cos2 0}.(12). 


We observe that if O is found by means of a pendulum, 
(11) determines the mass of the earth. 

The numerical value of m or oi^ajO is about 1/289, so that 
it is of the same order as c, and in the small term in (8) we may 
substitute o}H=^mO — mEja^, so that the potential of the 
earth at an external point is given by 

E Ea^ 

F = - + (im-c)~-^(co820-^) .(13). 


6*81. The Moon’s attraction. The forces exerted by the 
moon on the earth are equal and opposite to the forces exerted 
by the earth on the moon, and are approximately the same 
as if the moon’s mass M were collected at its centre of gravity. 
If r denotes the distance between the centres of the earth and 
moon, and 6 the moon’s north-polar distance, the earth’s 
attraction at the centre of the moon is composed of a radial 
dV IdV 

attraction P = — — and a force Q=z-— tending to increase 0. 

VT T Cu 

Hence from 6-8 (13) 


Jp 2 

P = 5 + 3(-^m-e)-^(co8i“0-i) 
Ed^ 

Q = — 2 (im — e) sin 6 cos 0 

r* 


( 1 ). 


The forces exerted on the moon are MP radially and MQ 
through the moon’s centre. The force system exerted on the 
earth by the moon consists of equal and opposite forces 
MPj MQ acting through the earth’s centre and a couple 
of moment MQr tending to cause rotation about an equatorial 
diameter, and since c — is positive the sense of the couple 
is such as to decrease 0. The moment of this couple is 

Eq2 

2(e — |w) sin d cos 6 


( 2 ). 
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6*82. We can obtain some information by comparing the 
results of 6*8 and 6*81 with MacCuUagh’s expression for the 
potential in 4*6 and the couple deduced therefrom in 4*62. We 
notice that the potential as given by 6*8 (13) is independent 
of the longitude, therefore when the attracting body is the 
earth MacCullagh’s expression 

E 


should be independent of the longitude. This is satisfied by 
taking a principal axis of the earth as the axis of rotation and 
writing .4 = ^, so that I— A sin^ 0 + (7 cos^ 0, and 



3(C-~^) 

2r^ 


(cos^e-i). 


A comparison with 6*8 (13) now shews that 


2/ 1 ^ 


The same result follows from comj>aring the couples of 

4*62 and 6*81. 


6*9. General solution of Laplace's equation. It has been 
shewn by Whittaker* that there is a general solution of 
Laplace’s equation of the form 


7 = ^^ f (z-h ix cos u + sin u,u)du, 

provided that / is such a function that differentiations with 
regard to x, y, z under the sign of integration are permissible; 
and it may be shewn that V can be expressed as a series of 
expressions of the types 



+ ix cos u + iy sin w)^ cos mn du, 
+ ix cos u -f iy sin u)^ sin mudu, 


* Math. Ann. Lvii (1902), p. 333. See also Whittaker and Watson, 
Modem Analysist ch. xvm. 
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where n and m are integers such that 0 < m < n; also that by 
changing to polar co-ordinates these may be reduced to solu¬ 
tions of the forms Xj* cosp(f>, Xj* &inp(f> in the notation of 

6 - 42 , 6 * 43 . 

EXAMPLES 

1. Find, in terms of /x, (/x), Pg (/x), Pg (/x). 

2. Prove that 




^ (*)+ sA-, + T-*- 


■2n~l"^2n-|-1^2n + 3‘ 


3. Prove that = 

4. Prove that P^^ P„ can be expressed in the form 

^0 ^m~n "t ^TO-n+l + • • • + Can ^m+n > 

where m>n and the c’s are constants. 

dPn dPn , d^Pft • 4.U X* 

5. Express , z and j” m the form 

A "f‘AiPl-f”A2P2 + A3P3-}“ ... . 

6. Prove that j d^ = n(n-M). 

7. Prove that 1 “ M*)= 0. 

2n(n-f 1) 
2n-f-1 

8. Shew that, if m and n are integers, the value of 


[M. T. 1912] 


[C. 1908] 


and 




is either 0, 2 or 


2n 


/ 


1 ^ dP^ , 


[M. T. 1908] 


2n-f-1' 

CO 

9. Prove that, if/(^x) is expressible in the form L Pn (i^)» then 




10. Prove that f P^ (/x) P„ (/x) dfi is equal to zero if m and n are 

Jo 

unequal and both even or both odd integers, while if m is even and n 
is odd, its value is 

i.3...(m~l) 3.5...n 

(m — n)(m + nH-1) 2.4...m 2.4... (n—1)* ^ 
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11. Explain why, if F (cos $) satisfies the equation 

^ I? ^ 

SO also does i^(Acos0 + /isin0cos<^+ vsin^sin^), 

where + [M. T. 1921] 

12. Shew that 

13. If F„ is a homogeneous function of x, y, z, of degree n satisfying 
V»F = 0, then 

+ Br* + F„ satisfies V^V*^ = 0. [C. 1891 ] 

14. Find the potentials of a surface distribution on a sphere in 

which the density at any point is proportional to the distance from 
a given plane. [C. 1904] 

15. A spherical surface of radius a is coated with matter of surface 

density proportional at each point to the cosine of the angle the radius 
vector from the centre to the point makes with the positive direction 
of the axis of x. Prove that the potential due to this matter at any 
external point is f^nara^xlr^, and at any internal point f^wcrx; where the 
origin is at the centre and ± or are the densities respectively at the 
points where the axis of x cuts the surface. [C. 1890] 

16. The surface density on a sphere of radius a is key. Find the 
potential inside and outside the sphere. 

17. The surface of a sphere of radius a whose centre is at the origin 
is covered by a thin layer of attracting matter, the surface density at 
any point P being proportional to ein* POZ. Find the potential at 
any point inside or outside the sphere. 

Prove that inside the sphere the lines of force are given by 
ylx = constant; zx^ = constant. [London Univ. 1925] 

18. Find the surface and volume distribution of matter for which 

F = A;r*^(3cos2 I) when a, 

and V^ka (3 cos^ ^ — 1), when r > o, 

fc, p and a being constants. [London Univ. 1936] 

19. If the mass per unit area of a thin heterogeneous stratum of 
attracting matter placed on a sphere of radius a is pz^ja^^ prove that 
the potential at a point within the sphere is 

g^7r/> [5z^ja^ — ZrHja^ + 7z), 

where r* = .r* + y® + 2*, and find the potential at a point without the 
sphere. [M.T. 1903] 
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20. Prove that Laplace’s operator V® is invariant for any change of 
rectangular axes, and that a solid harmonic (or a surface harmonic) 
of degree n transforms into a solid harmonic (or a surface harmonic) 
of degree n in the new co-ordinates. 


21. Prove that, if x be real and positive, 

ri, ^ 21 _ 

' ' J ~_l (cosh2.T — /x)i (2n4-1) 

(2) J ^(co8h2a; —/x)ip„(/i)d/x 

2i 2* 

”" (2n + 1) (2n -f 3) * {2n + 1) (2n - 1) »‘ 

[C. 1908] 


22. If Pn is Legendre’s coefficient of order n, prove that 

dap 71/2 

= n (n- 1) -f S^(2n~ 4r -f 1) {r (2n - 2r -f 1) - 2}P„_2r 

when n is even, and that when n is odd the summation is from 1 to 
i(n-l). [M. T. 1904] 

23. Express xyz (x^ + + s®) in the form 

where F^, F^, Fj, Fq are solid harmonics of the orders indicated by 
the suffixes. ^London Univ. 1931] 


24. Express in the form Fg F4 + F2, where F„ denotes a spherical 

solid harmonic of degree n, a function which satisfies Laplace’s equation 
within a sphere r = a, is finite at the origin, and is equal to x^y^z on 
the sphere r = a. [London Univ. 1926] 

25. Obtain a rational integral harmonic which has the value Az* 
at points on the surface of a unit sphere with its centre at the origin. 

26. A thin layer of matter is placed on a sphere of radius a whose 

centre is O, the density at any point P being proportional to cos* AGP 9 
where OA is a fixed radius. Determine the potential at all points 
inside and outside the sphere. [London Univ. 1927] 

27. Matter is distributed in a thin layer over the surface of a sphere 
of radius a and centre at the origin, the density being proportional 
to (x + y)® and the whole mass being m. Shew that the potential at an 
external point is given by 

y ^[1 + +6ary-28»)J 

and find the value at an internal point. [London Univ, 1938] 
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28v. If the surface density at the point {x, y, z) on a sphere of radius a 
is x^y — xy^, prove that the potential is 

at an internal point, and 

47ra® (x-^y) 

\ 9 ^ wi-r# / 

at an external point. [London Univ.] 


1 \ 
7r7 


29. The sphere x^y^ + is covered by a thin layer of attracting 

matter, the total amount of matter being M and the surface density 
at any point (X, F, Z) of the sphere being proportional to 

y2^2 4_^2X2-fX2F2, 


which may be written as 

j (a4 - X 4 _ y 4 _ ^4 ^ 3 y2^2 ^ 3^2x2 + ZX^ Y^) . 


Find the potential of the layer at any external point. 

Shew that the attraction on a small mass m, placed outside the 
sphere on the line a; = y = 2 at a distance r from the origin, is 


^ • [London Univ. 1925] 


30. If F„ is a rational integral harmonic function of degree n and 
=i y^ 2*, and U iso. homogeneous function of a:, y^ z of degree 2n, 

n 

shew (by expressing XJ in the form X op otherwise) that 

p-==0 


the integral being taken over the surface of the sphere r = a. 

[M. T. 1939] 


31. The density of a solid sphere of radius a, referred to axes through 
its centre, is x^yz. Shew that the potential of the sphere at an external 
point {x, 2/, 2) is given by 

jHt (ay /o^»l 

'^■■■63“ 111 \rj 5 -\r) j ’ 

whore = yz, H4 = %x^yz — y^z — yz^, [M. T. 1908] 

32. The potential over a sphere of radius a is equal to Axyz^^ and 

is due to some surface distribution a. Find the potential at an internal 
point, at an external point and the value of a. [C. 1913] 

33. Find the potential in all space when the potential at any point 

of the surface of a sphere is yH^. [C. 1906] 
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34. The density of a solid sphere of radius a at the point 
is where ik is a constant. Prove that the potential at any point 
inside the sphere is 

[| (6z» - Sf-*) (9a2 - 7r“) + 7o‘ - 3r*], 


where = x* + y® + 2 *. 


[London Univ. 1926] 


35. The density of a solid sphere of radius a is Po + Pi^ 2 > where 
Po» Pi constants and S2 is a spherical surface harmonic of the second 
degree with origin at the centre of the sphere. Prove that the potentials 
at internal and external points are 

f7ry/>o (3a®-r*) + fn-ypir®/S', + ij*-7rypir®/S', 


and 


Jwypoy + in'^ypi ^ 

[London Univ. 1931] 


36. The density of a solid sphere of radius a is proportional to a?® 2 ®, 
the centre being the origin. Prove that the potential at an external 
point is 


y_H + --- 


la^ 

99r« 


6 

99 




"^33 r8 r 


[London Univ. 1926] 


37. The matter of a heterogeneous solid sphere is arranged 
symmetrically about a diameter. It is found that the normal com¬ 
ponent of attraction due to the sphere at a point P on its surface where 
the radius makes angle 6 with the diameter A-\-B cos* 6, Prove 
that the tangential component of attraction at P is 

^Psin^cos^{2 + 7cos2^). [C. 1907] 

38. A given mass is distributed in a thin layer over a spherical 
surface, symmetrically with respect to a diameter AB and to the 
corresponding diametral plane, so os to produce inside a field, such 
that the force at any point towards AB is proportional to the distance 
from it. Shew that the force parallel to A P is proportional to the dis¬ 
tance from the corresponding diametral plane. Shew also that if the 
surface density in the diametral plane is zero, that at A is three times 
what it would be if the mass were uniformly distributed. [C. 1901] 

39. Shew that the potential of a imiform circular ring of radius a 
and mass M is 

■w j 0 

where the origin is at the centre and the axis of z along the axis of 
the ring. 
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40. V©rify, by using the formula© of 6*28 for , that the expressions 
Fi, F 2 of 6*32, for the potential of a circular disc, are equal when r = a, 
by proving them both equal to 

2m [ 1 + {/X + — 1) cos dfj) — 27nnafx. 


41. Prove that, if the potential of a body symmetrical about an 
axis at a point on the axis at a distance z' from an origin is 
then the potential at any point outside the body is 


- f {2 + + r2) cos«//} 


where r is the distance of the point from the origin, and z the projection 
of r on the axis. [M. T. 1887] 


42. Matter of mass M is distributed on the surface of a sphere whose 
centre is O and radius a, so that the density at any point is proportional 
to the square of its distance from an external point (7, where 0(7 = &. 
Prove that the potential at an external point P is 


M 


2 a^b x] 
3(o* + 6»)r»j’ 


where OP = r and a; = r cos POG. 


[C. 1897] 


43. The density of a solid sphere at any point P varies inversely 
as OP, where 0 is an external point distant / from the centre. Shew 
that the potential at any external point is 

' r oVfr) {2n+l){2n+Sy 

where M is the mass of the sphere, a its radius, and r, B polar co¬ 
ordinates whose origin is the centre and whose axis is the line from 
the centre to O. [M. T, 1909] 


44. The density of a thin spherical shell, centre (7, radius a and 
mass M, varies inversely as the distance from an external point 0. 
Shew that the potential at any external point P is 

IM faVf( _ f _ 

2 a Jo \x(r^ —2rxcosB-hx^) 

and at any internal point P is 

IM fr^/ff / _ 

2 a jo \x{r^--2rxcosB + x^) 

where r = OP, /= 00 and B is the angle POO. [M. T. 1895] 

46. A disc of radius a is of a mass M and its surface density is 
proportional to the mth power of the distance from its centre (m + 2 > 0). 
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Shew that the potential at points whose distance (r) from the centre 
is greater than a is 


yJf|-4-(m + 2) 2 (-1)" 


1.3.6...2n~l a2« 

! (m 4- 2n + 2) 



[London Univ.] 


46. A thin circular disc of radius a is such that the mass per unit 
area at distance r from the centre is — Find the attraction 
of the disc at any point on its axis. 

Hence shew that the component attraction parallel to the axis of 
the disc at any point whose co-ordinates are r, d, referred to the centre 
of the disc as origin and its axis as polar axis can be expressed in the 
form ^ 


if r>a. 


27rifca22(-l)«- 

1 


1.3.5...2n- 
(n-l-1)! 




_1 (cos d)j 

[M. T. 1920] 


47. Prove that at a point in free space, where r>a, the potential 
of a homogeneous gravitating hemisphere of radius a and density p is 


V = 27rypa^ 




(2n-4)! 


2^”-‘n!(n-2)! r* 




where r is the distance of the point from the centra of the base and the 
axis of harmonics is the axis of symmetry of the hemisphere drawn 
from the centre into the interior. [London Univ. 1931] 


48. Show that the potential of an oblate spheroid of density p 
and major and minor semi-axes a, c at large distances is approximately 


4:7rypa^c 


r 1 Po 

_i.3 r 


3.5 5.7 A 


[C. 1929] 


49. A uniform plane annulus of mass M is bounded by circles of 
radii a, 6, large compared with the dimensions of a spheroid whose 
centre of gravity coincides with the centres of the circles. Shew that, 
if the normal to the plane of the annulus makes an angle 6 with the 
axis of the spheroid, they exert on one another a couple of magnitude 


SyM{A-G) 
ah (a 4- h) 


sin 6 cos 6, 


tending to increase 0, where A, A, C are the principal moments of 
inertia of the spheroid, and y is the gravitation constant. 

[M. T. 1931] 


50. A imiform nearly spherical solid of density p has the surface 
r = a (14- cPg) ^ if® boundary. It is surrotmded by liquid of volume 
47 t( 6® —a®)/3 and uniform density a. Shew that, provided the solid 
is completely covered with liquid, the equation of the free surface is 
r = 5 (14- rjP 2 )t where 


3(p~or)a*c 

^"'62{6(/)-a)a» + 2a6»}’ 


[M. T. 1936] 



Chapter VII 

ATTRACTION OF ELLIPSOIDS 


7*1, Homoeoids and Focaloids. A homoeoid is a shell 
bounded by two similar and similarly situated concentric 
ellipsoids. A focaloid is a shell bounded by two confocal 
ellipsoids. 

Let a, 6, c be the semi-axes of the internal surface of a thin 
homoeoid and a -}- da, 6 + dft, c -f- dc those of the external surface. 
Let 0 be the ^common centre and let a radius OPQ cut the 
surfaces in P, Q, where OP — r and OQ ^V'Vdr. 

Since the figures are similar 


d'T_^da^db^dc^ 


( 1 ). 


Also the tangent planes at P, Q are parallel, and if p-^dp 
denote their distances from 0 and I, m, n the direction cosines 
of the perpendicular to them, we have 

^ ^2^2 ^ 52^2 ^ ^2^2 

and {p -i- dpY = (<x + daY + (6 + dbY + (c -f dcY 

= (1 -f {aH^ -f bhn^ -f c^n^), from (1), 
=:{l+kYp^. 

It follows that dp Ip = k .(2); 

and dp, the distance between two parallel tangent planes, 
is the thickness of the homoeoid at P. 

The volume of the internal eUipsoid is ^nobc and that of 
the external ellipsoid is f tt (a+da) {b-\-db) (c -hdc) or, from (1) 
^7rabc{l-^kY, Therefore, neglecting higher powers of k, the 
volume of the thin homoeoid is ^nkabc, when k — daja — ^tQ. 

Thin homoeoids are said to be confocal when their inner 
boundaries are confocal. 

For a focaloid for which the semi-axes of the internal and 
external ellipsoids are a, b, c and a', b\ d we have 

a'2=a2 + A, c'^^c^ + A.(3). 


RN A 


II 
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If the focaloid is thin and we take a-{-da, 6 -fd 6 , c + dc as 
the semi-axes of the outer ellipsoid, we have from (3) 

ada = bdb = cdc==^X .(4). 


And if p,p-\-dp are the distances from the centre of parallel 
tangent planes, and Z, m, n direction cosines of the normal, 
we have ^2 _ ^ 2^2 ^ 52^2 ^ ^ 2^2 


and {p + dpY = + A) Z^ + (6^ (c^ + A) 

so that pdp = \\ ..(6). 


The volume of the thin focaloid is 


( 6 ^ + A) (c^-fA)} —f 7 ra 6 c 

-|,A«6c(A+^ + i.) ..,(6). 

If p is the density of the matter and the thin shells are 
regarded as strata of matter condensed on the inner boundary, 
the mass per unit area is in either case pdp, and in the case 
of the homoeoid this is kpp, varying directly as p \ but in the 
case of the focaloid it is ^pA/p, varying inversely as p. 


7*2. Attraction of a homoeoid at an internal point. With 
any point 0 inside the homoeoid as vertex, draw a cone of small 
solid angle dw intersecting 
the homoeoid in frusta 
PQQ'P', BSS'P\ If pis the 
density, the mass of an ele¬ 
mentary slice of the cone of 
thickness dr at a distance r 
from O is pr^dwdr, and its 
attraction at 0 is pdcodr. 

Consequently the attrac¬ 
tions of the two intercepted 
frusta of the cone are p dwPQ 
and pdoiRS in opposite 
directions. But since the ellipsoids are similar, they make 
equal intercepts on any chord so that PQ===^It8 and the 
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attractions of the frusta are equal and opposite. By taking 
cones in all directions round O, we see that the resultant 
attraction of the homoeoid at O is zero. 

We appear to have neglected wedge-shaped elements at 
P and Q in integrating pdcjdr with regard to r, but the further 
integration with respect to dco implies that dcu tends to zero 
when we regard the integral as the limit of a sum, so that 
these elements contribute nothing. 


7*21. Potential of a homoeoid at an internal point. Since 
the attraction at an internal point is zero, therefore the 
potential is constant and equal to its value at the centre of 
the boundary surfaces. If we take a cone of small solid angle 
doj with its vertex at the centre, the contribution to the 
potential of the frustum of the cone intercepted by the 

homoeoid is easily seen to be J prdojdr=^\p{r^ — r^)da), 

where are the radii of the inner and outer surfaces in the 
direction of the cone. If we take (a, 6, c) and [ma, mb, me) as 
the semi-axes of the outer and inner surfaces = mr^ and the 
potential is 


Ip (l—m^) jr^dco .(1) 


integrated over the surface of the outer boundary. 

From the equation of the ellipsoid in polar co-ordinates, 
we have 

1 8mV\ CO8®0 


and putting (ia> = sin OdOcUf), we have 


|r^dct> = 8 


sm9ddd<f) 

. /cosset sin^cbX coB^O 

BinOdOdt 

sin*0 cos*0 ^ (&m^0 ^ cos^fl'^ 

0 + 


^ smddddt 

sin*0 cos*0 /&m^0 cob^0\ 


0 Jo ^ 


11-2 
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where ^ = But 


f" * 


^ 7T 


J A] 

0 ~2^y(ABY 


therefore 



sinddO 


cos^ 6\ I sin^ 6 

~W~ 






cos^ d\\ 

C 2 /) 


Put tan^ d, and we find that 

r c°° du 

jr doj-^iTobcj^ y'{(a*+«)(6* + M)(c*+«)} 

r 00 

We shall denote -Trr ^—rTTa- h 

J 0 \/{(«^ + u) (b^ + u) (c2 -^-u)} ^ 

substituting in (1), we find that the potential of a homoeoid 

at an internal point is 

V = rrpabc (1 — m^) /. 

But the mass M=^^7Tpabc (1 —m^)y so that 

F = |if/f-^3 .(3). 


If the homoeoid be thin, then m is nearly equal to 1, say 1 — e, 


«, that and 

l-m3 3€ 3’ 


V=iMI .(4). 


7*22. Attraction of a thin homoeoid at an external point. 

We shall prove first that this is directed along, the normal to 
the confocal elhpsoid through the external point. 

Let QP, QT be the normal and the tangent plane to the 
confocal at Q. Let the normal QP meet the plane of contact 
PB' of the enveloping cone of the homoeoid from Q in. P. 
Through P as vertex take a cone of small solid angle do> 
cutting the shell in elements MN, M'N', and let MM' meet 
the tangent plane to the confocal at T. Then the pole with 
regard to the homoeoid of the tangent plane QT to the con- 
focal lies on QP, but it also lies on the polar planeof Q, i.e. RR' ; 
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therefore P is the pole of the tangent plane QT. Therefore 
(MPM'T) is a harmonic range. But PQT is a right angle, 
therefore QP bisects the angle MQM\ 

Now if fjL, denote the masses of the elements MN, 
of the shell, fju: ix' — MP^: PM'^. But since QP bisects the 
angle MQM\ therefore MP : PM' = MQ : QM', 

Therefore or the attractions at Q of the 

elements MN, M'N' are equal; and they are equally inclined 
to QP, therefore their resultant acts along QP, Since the 
same holds good for all like pairs of elements into which the 
homoeoid can be divided, it follows that the resultant attrac¬ 
tion at Q is along QP. 


Q 



It is also evident that the plane of contact of the enveloping 
cone from Q divides the homoeoid into two parts which exert 
equal attractions at Q. 

Since the resultant force at a point external to the thin 
homoeoid is always directed along the normal to the confocal 
ellipsoid through the point, it follows that the external equi- 
potential surfaces are the confocal ellipsoids. 

It follows that any two thin confocal homoeoids have the 
same external equipotential surfaces, viz. the ellipsoids con¬ 
focal with them and external to both, and therefore by 
4*31 (iii) their attractions at any point are proportional to their 
masses. 

Now let a, 6, c be the semi-axes of a thin homoeoid of mass M 
and density p, and let a', 6', c' be the semi-axes of the confocal 
thin homoeoid having the point Q on its outer surface and of 
the same density. If p be the central perpendicular on the 
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tangent plane at Q we may take dp as the thickness and pdp 
as the mass per unit area of this confocal homoeoid. Since 
normal attraction increases by 4^7rpdp in crossing from inside 
to outside the surface (3*7) and the attraction is zero inside, 
therefore the attraction at Q of the outer homoeoid is ^rrpdp, 
but dpip — da'ja' (7-1), and the mass M' — 4:7Tpb'c'da' (7*1), so 
that the attraction of the outer homoeoid is M'pja'b'c', 

But the attractions of the two homoeoids are proportional 
to their masses, therefore the attraction at Q of the original 
homoeoid is Mpja'b'c^ where a\ b\ c! are semi-axes of the 
confocal through Q and p is the central perpendicular on the 
tangent plane to the confocal through Q. 

We may now deduce the potential of the homoeoid at an 
external point Q, for since p is normal to the equipotential 
surface through if F is the potential 

^ = force in direction p 
- ^Mpla'bW 

or dV= —Mpdp|\/{(a^-h^)(b^-h^)(c^-h^)}y 

where xy(a^ -f A) + t/^/(b^ + A) + z^/(c^ -f A) = 1.(1) 

is the confocal through Q, But 

^ (^2 ^ ^ (c^ 4- A) n^ 

= dH^ + b^m^ + c^n^ -f A, 
so that 2pdp — dX 

and dV = — \MdXj\/{{a^ ^X) {b'^-{-X) {c^ -^X)]. 

Hence F= JJlf '^{{a^+u)(b‘^+u){c^+u)} . 

where the lower limit A is given in terms of the co-ordinates 
X, y, z of as the greatest root of the cubic in A, (1), and the 
upper limit is determined by the consideration that F -> 0 
as A->oo. 

7*3* Attraction of a homogeneous solid ellipsoid at an 
internal point. Let jP be a point inside the ellipsoid. If we 
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draw a similar, similarly situated concentric ellipsoid through 
P, the matter outside this surface constitutes a homoeoid and 
exerts no attraction at P, so the 

problem is reduced to finding the - 

attraction of an ellipsoid at a point ^ 

P on its surface, ( II \ 

Let be the ( nr j 

ellipsoid through P and /, g, h the V j J 

co-ordinates of P^ Let r be the length 
of any chord PQ and Z, m, n the 
direction cosines of QP, so that the 

co-ordinates of Q SiTef--lr, g — mr, h^nr. Then since Q and 
P both lie on the ellipsoid, we have 


72 ^ 

_ O "t" 10 




fl gm hri 


_ 4. —J— 

2 ^ 62 ^ C2 


Taking a cone of small solid angle dco with vertex at P and 
axis along PQ, the attraction at P of the matter in this cone 

is ^pdo}dr=^prdoj, where p is the density. Hence the com¬ 
ponents of the resultant attraction at P are given by 

X=—pjzrdoi, Y=—p^mrdw, — pjnrdo) ,..{2). 

We must now substitute for r from (1) and integrate with 
regard to dw for cones in all directions on one side of the 
tangent plane at P. But since every straight line through P 
has a real intersection with the ellipsoid, if we integrate in all 
directions round P we should take the ellipsoid twice over. 


Therefore X = — 


Z2 m2 


where Z, m, n may now be regarded as co-ordinates of a point 
on a unit sphere, dco an element of area of this sphere and 
the integration extends to the whole surface. The terms in 
Im and In contribute nothing to the integral because terms Im 
and Z (— m) occur equally and cancel one another. 
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Therefore —pf 


Similarly Y — —pg 


and 


— p/i 



?* 

a* 


¥ 

m* 

n2 

a* 

+ 6^ 

-f-o 

C2 


m* 

F 


¥ 

m* 


o* 


“^“2 

C2 


»* 

F 


¥ 


7i2 

— 

-f- "Tir 

H-ST 

a2 

62 

C2 


d!a)= -Apf, say 


do) — — Bpg 


doj = — Cph 


.(3). 


We notice that A, B, C are functions of the ratios of the 
axes, so that it is immaterial whether a, b, c are semi-axes of 
the given ellipsoid or of the similar one through P. 

7*31. To shew that 

du 


A = 2'7Tabc 


poo 

Jo 


We have 




^ -h u)^ (6^ + u) (c^ -f u)} ‘ 
P 


b^^c‘ 


dojy 


and if we use polar co-ordinates this may be written 

p^TT piTT COS ^0 


A^S 


J 0 %j 

piTT p 


coB^d . «./cos^<i sin*<i\ 


sin6ddd<l> 


= 8 


cos^e . ^ 

—^smadOdt 


0 J 0 

piTT 


cos^ 0 sin^ 0 


lco&^0 ^ sin^^"^ 


where <=tan 


= 4ir 


62 


—5—smflap 


«2 


Ij / COB* 6 sin* 0\ I cos* 0 sin* d\\ ‘ 

0 A/il"^ “^A“^ 1^/) 
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7-32] 

Put o® tan® 6==u, and on reduction we get 

4 « I r* du 

A-2vabcj^ ^{(a® + «)®(6®+tf)(c® + tt)}. 

and similar expressions for B and C. 


7'32, Spheroids. We see from 7*3 (3) that 
j4 + jB + (7= 

and Aa^ + Bb^ + Cc^ 


doj 


.( 1 ). 


Two of the axes of a spheroid are equal, and if a = 6 then 
A=zB, X=--Apf, Y=--Apg, Z^--Cph, 
and (1) become 2^ + O = 47r, 

dw 


2-4a2+Cc2 = 




'TT /* 27 r gin 


Jo J 

= 277 


sin^ B cos^ B 
du 


where u = cos B, 

There are two cases: 

(i) Oblate spheroid, a>c. 


1 . / 1 2 ' 

ov 




= 4^a y(| . 


V(l-e*) 

where e is the eccentricity of the generating ellipse. 

(ii) Prolate spheroid, oa. 

2Aa^ + Ce* - loc 

2Aa +CC -y'(c«_o*)^°8c-V(c*-o») 

= 2wc»—— log, 


•( 2 ). 


‘ 1 _< 


•(3). 


where a* = c* (1 — e*), so that e is the eccentricity of the generating 
ellipse. 

A and O can then be found in either case, since 2^ + C = 477. 
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7*33. Attraction of an oblate spheroid of small ellipticity. 

When the ellipticity € is small, since c = a(l —c) and c* = a’®(l — e®), 
therefore e* = 2c approximately. 

Hence from 7*32 (i). 


2^ + C (1 - e») = ^ tan-t 


e 

V(l-e') 


“ I 31 -e»'^ 6 (l-e»)« •" 

Whence, since 24 + C = in. 

But 2A + C = 47r, so that ^ (1 — fc). 

Hence in this case 

X= ~|wp(l““|€)/'| 

r=~^,rp(l-.fc)^h.(1). 

JS=-^irrp(l + i€)h) 


These results are of importance in connection with figures of equi¬ 
librium of rotating liquid. 


7*34. Attractions of a nearly spherical ellipsoid. Let 

b = a(l — 7j) and c = a(l —c). Neglecting higher powers of r) and c 
than the first, we may put -4 = + Ftj -f F'c. 

If the ellipsoid were turned through a right angle about Ox, b and c 
would be interchanged, so we must have F:=^F'. Also an oblate 
spheroid for which 6 = a and ?? = 0 is a special caae, for which, from 7*83, 


therefore 

Hence 

Similarly 


and 

Therefore 


A=|ir{l-|c), 


Aq^^TT and — ^TT, 

4=Jjr(l-fi7-|«). 

<7=111 (l-fTj + Je). 

X=-inp(l-ir,-i,)A 

Y~ -irrp(l + iri-ie)gt .... 
Z=-iirp(l-h + te)Aj 


( 1 ). 
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We may obtain similar results in a symmetrical form, if we take k 


as the mean radius J (a + 6 + c), then 



11 

22a~6- 



^6 a 

) 



6 a —AjV 
^ 5 ifc ' / 


Similarly, 

B=u{ 

66-A;\ 
5 k ) 

• .. 

and 


6c-l;\ 

5 





7*35. Potential of a homogeneous solid ellipsoid at an 
internal point. From 7*21 (3) it follows that, since a solid 
ellipsoid can be divided into similar homoeoids, the innermost 
being of zero dimensions, by putting m = 0, the j)otential at the 
centre of a homogeneous solid ellipsoid is 


To find the potential at an internal point (x, y, z), we make 
use of the fact that the attraction components are given by 

X,Y,Z=-Apx,-Bpy,-Cpz (7*3) .(2). 

Hence, if V denotes the potential, 


7T7 

dV=^ — dx + -;^dy + -7^dz 
ax ay dz 

= Xdx-h Ydy-hZdz 

= — p (Axdx + Bydy 4- Czdz). 

Therefore V = \p{D — Ax^ — By ^~ Cz^) .(3), 

where the constant of integration \pD must be the potential 
at the centre. Hence by substituting for D from (1) above, 
and for A, B^ C from 7*31, we get 

/•co I 

V 'irpabcj ^ ^ ^^ (^2 ^ 

.( 4 ). 

\ o®+« c^+u/ 
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It follows from (3) that the equipotential surfaces inside 
the ellipsoid are ellipsoids 

-f Cz^ = const.(6), 


similar to one another though not similar to the boundary 
of the solid. 

7*4, Ivory’s Theorem.* Attraction of a homogeneous 
solid ellipsoid at an external point. If x, y, z and y\ z' 
are the co-ordinates of points P and P' one on each of two 


confocal ellipsoids 

x^jd^ -f y^jb^ -h =1.(1), 

and -f y^jb'^ 4- z^jd^ =1 .(2), 

then P and P' are called corresponding points if 

xja — x'la', ylb=^y'lb\ zlc = z'j& .(3). 


Let R8 be an elementary strip of the first ellipsoid of cross- 
section dydz, parallel to the 
x-axis; and let B'S' be the 
corresponding strip of the 
second ellipsoid of cross- 
section dy*dz\ so that 

dydzjdy'dz' = bcjb'c'. 

Then if /' (r) denotes the 
law of force at distance r, 
and p the density of either 
eUipsoid, the component 
parallel to the a;-axis of the 
attraction at P' due to the strip B8 

= —pdydz {r)co^P'Q8dXy 
when Q is the position of an element dxdydz. This 
, = -pdydz^r{r){^-^dx=pdydz{f{P'8)-f{P'R)}. 

* Sir James Ivory (1765-1842). Scottish mathematician. 
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7 - 4 ] 

In the same way the a:-component of the attraction at P 
due to the strip B'S' 

=^pdy'dz'{f{P8')-f(PR')}. 

But it is easy to deduce from (3) that for any pairs of 
corresponding points PS' = P'S and PR' = P'R. 

Hence the ratio of the a:-component of the attraction of 
RS at P' to the similar component of R'S' at P is 

dydz be 

dy'dz' b'c' . ^ 

If in like manner we take all such strips as JKS of the first 
ellipsoid and all such strips as R'8' of the second, we get 

X: J:' = 6c:6V .(5), 

where X is the ^c-component of the attraction of the first 
ellipsoid at P' and X' is the o^-component of the attraction 
of the second ellipsoid at P, 

This result is Ivory’s theorem and is true for any law of 
force. 

Now, from 7*3, X' = — A* px, where Pis the point {x, y, z) and 
A' is the same function of a', 6', c' that A is of a, 6, c. 

Therefore 




abc X _ 


ohc . f , 

7-.Apx, 


from (3), 


where x\ y\ %' are the co-ordinates of P'. 

Similar relations hold for the y and z components of attrac¬ 
tion, so that the attraction at P', {x\ y\ z'), of the ellipsoid 
(a, 6, c) has components 


X=-^A'pz' 

abc 




( 6 ), 


where a', b', c' are the semi-axes of the confocal ellipsoid which 
passes through P', and A', B', C' are the same functions of 
o', b', c' that A, B, C are of a, b, c. 
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7*41. MacLauriii’s Theorem. ♦ The aUrmtions of two con- 
focal ellipsoids at a point external to both are proportional to 
their masses and in the same direction. 

In 7*22 this was proved to be true of confocal thin homoeoids; 
that it is true for homogeneous solid confocal ellipsoids follows 
as a direct corollary from 7-4 (6), if we regard P' as a fixed 
point and consider the attractions at P' of two different 
confocal ellipsoids, of semi-axes a, b, c and density p and of 
semi-axes and density the x, y, and z components 

are aU in the ratio pabc : p^dib^^c^ and are therefore propor¬ 
tional to the masses and have resultants in the same direction* 


7*42. Potential of a solid homogeneous ellipsoid at an 
external point. Let a, 6, c be the semi-axes of the given 
ellipsoid and x, y, z the co-ordinates of the external point P. 
Let a\ 6', c' be the semi-axes of a confocal ellipsoid which 
passes through P, and F, F' the potentials at P of the first 
and second ellipsoids regarded as of the same density. 

Then by MacLaurin’s theorem (7*41) 

F: F' = a6c:aW. 

But, from 7*35 (4), 


^'^npa'b'c' 


Therefore 


0 ■\/{{a'^+u){b'^->ru){c’^+'u)] 

( x^ z^ \ 

1 


p 00 

V - TTpabcJ ^ ^(,'2 ^ + uj) 

\ a’^+u b'^+u c'^+uj ' ’ 

but since the ellipsoids are confocal, we have 
a'*, 6'*, c'2 = a* + A, 6® + A, c*+A, 
where A is the largest root of the cubic 


, y® . 2* 

: + r 5 - . ~x + ."^ = 1 


a*+A 6*+A ca + A 


.( 2 ). 


* Colin MacLaurin (1698-1746). Scottish mathematician. 
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7 * 6 ] 

Substituting in (1) for a'^, 6'^, c'^ and writing v for A + 
we get 

f" 1 

V - TTfMbcj ^ ^ ^ + J,) (c 2 ^ 

/, a:® 2/® z® \ _, 

where the lower limit A is determined by (2) as stated* 


7*5. Attraction of an infinitely long homogeneous solid 
elliptic cylinder. This can be deduced from 7*3 and 7*4 by 
making c oo. 

(i) Internal point. The components of attraction at a point 
{x, y, z) are given by 

X——Apx, Y——Bpyy Z = 0; 


where, from 7‘31, 

A = 27rab 


poo 

Jo v'{(a® 


du 


{b^^u)y 


By substituting for + u, we get 


A = 47Ta6 


rl/a 

Jo A 


vdv 


y/{l^v^{a^-b^)] 


or 


A=: 


4:7rb 
a -{-b 


Similarly 




4trra 


and 


X=-> 


Anpab x 
a + 6 a’ 


7 = 


^TTpab y 
a-\-b 'b 


.(1). 


(ii) External point. In like manner from 7*4 we find for 
the attraction components at an external point (a;, y) 

^ ^irpab X ^ ^irpab y 

- ^ f 1 J _ — , , , . T ^ . 

a -f 0 a a +o o 

where a', b' are the semi-axes of an ellipse confocal with the 
cross-section of the cylinder and passing through {x, y). 
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7*6. Equilibrium of rotating liquid. When a mass of gravi¬ 
tating Hquid rotates uniformly about an axis, it is concjeivable 
that for a certain form of the free surface the liquid particles 
may be in relative equihbrium. Since the resultant attraction 
on a particle depends on the form of the boundary surface, 
which is unknown, the problem does not admit of a complete 
solution. But assuming the liquid to be homogeneous, it can 
be shevrn that certain forms are possible forms of relative 
equihbrium. 

Thus if p be the density and a> the angular velocity, whether 
spheroidal or eUipsoidal forms are possible forms of relative 
equihbrium is found to depend on the numerical value of 
(jD^j2'iTp,* We shall illustrate the theory by considering two 
simple cases. 

7*61. To shew that an oblate spheroid of small ellipticity 

is a possible form of relative equilibrium of a mass of uniform 
liquid rotating about an axis with a small uniform angular 
velocity, (MacLaurin’s theorem.) 

The problem can be treated as a statical one if we compound 
with the force of attraction the reversed effective force, wh per 
unit mass, where r denotes distance from the axis of rotation. 

The equation for the pressure is then 

— = Xdx4- Fdy + Zdz 4- o)^ (xdx-hydy) 

P 

dn 

or -Jt- = — (xdx ^-ydy) — Cpzdz + (xdx 4 ydy), 

P 

assuming that the boundary is an oblate spheroid. 

But the pressure must be constant over the boundary, so 
that its equation must be of the form 

(o}^-^Ap) (x^ y^) -- Cpz^ = const .(1). 

If a, a, a (1 — c) are the semi-axes of the oblate spheroid and 
the elhpticity e is small, its equation is 

x^+y^+j^^ = a’^ .( 2 ). 

♦ Besant and Bamsey, Treatise on Hydrornechanics, Part I, ch. viu, 
where a discussion of the subject with full references will be found. 
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By comparing ( 1 ) and ( 2 ), we get 

to^-Ap=-Cp(l-2e) .(3). 

But for an oblate spheroid of small eUipticity, from 7’33, 
^=|7r(l-fe) and- C7 = |7r (1+te); 

15 

and, by substituting these values in (3), we find that e = —, 

loirp 

and with this relation between € and oj the condition is 
satisfied. 

7-62. Jacobi’s ellipsoid.* An ellipsoid with three unequal 
axes is also a possible form. 

For an ellipsoid the pressure equation is 

— = — Apxdx — Bpydy— Cpzdz + {xdx + ydy), 

P 

where A, B, C are given by 7*31. 

The surfaces of constant pressure are therefore 

{(A)^ — Ap)x^-\‘ (a>^ — Bp) y^ — Cpz^ = const.( 1 ). 

And if we assume the free surface to be 

x^ 

. 

a^ 

we get by comparing ( 1 ) and ( 2 ) 

a^{w^--Ap)^b^{io^-Bp)=^-c^Cp .(3). 

By ehminating ci>, we get 

a262(B-.A) = {a2-62)c2C.(4). 

If we put D = {(a2 4 . u) (b^ 4 - u) (c^ 4 u)]^, then, from 7*31, 

A = 277060 I 7-5 - ™=r, J5= 277060 I TTS-;- 

Jo (o24w)Z)’ Jo (b^-¥u)D^ 

^ ^ .L du 

and (4) becomes 

(a^-b^)r^{ _ ^ _^ 1 = 0 . 

^ ^Jo i) 1(0®+«)(&*+«) c»+uj 

♦ Carl Gustav Jacob Jacobi (1804-1851). German mathematician. 
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Hence, if a+ 6, the lengths of the axes must be such as to 
satisfy the equation 


"^^udu /1 1 1 ^ \ /V 

0 


( 6 ). 


If a and b are given, then (6) is an equation for determining 
c, and since the left-hand member is negative when c is small, 
and becomes positive as c increases, there must be a real value 
of c which satisfies the equation. 

Further, is given from (3) in the form 

where M is the mass of liquid. So that when a + 6, 


O) 


2 = 



udu 

(a^-f w) (b^-\-u)D' 


and since this expression is positive there is a real value of to, 
thus establishing that the ellipsoid is a possible form. 

Further, from (3), 


a2^2 ^ ( J[qj2 _ (7^2J p 


:|Jf (a^ —c^) 


i 


udu 


0 {a^-^u)(c^-^u)D' 


so that for w to be real we must have c < a. Similarly c < 6, so 
that the eUipsoid must rotate about its least axis. 


7*63. Bxample. A rigid gravitating sphere of radius a and uniform 
density p is surrounded by a layer of gravitating liquid of volume 
^TT (6* — a®) and uniform density a and placed midway between Vwo spheres 
of 7nas8 M whose distance apart (2/) is so great that powers of bjf above 
ibe third may be neglected. Prove that the boundary of the liquid is a prolate 
spheroid of small elUpticity 

45M6» 

47r/®{5a*(p-<7) + 2a6*}' 

and determine the small angular velocity of the whole mass wMch, to this 
approximation, would make the boundary a sphere. [P. 1933] 

Let A, Bhe> the centres pf the two spheres of mass M and 0 that of 
the sphere of radius a. Let (? be a point in the liquid at distance r 
from 0, where the angle Q0B=x$, and let OB be the axis of z. 
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It is easy to shew that the potential at Q due to the two spheres of 

M M 2M / r* \ 

mass ikf, being + gg»to [1 4-J2 ^zj * neglecting higher 

powers than 6 ®//^, This represents a small disturbing field of force 
which distorts the boundary of the fluid from what would otherwise 
be a sphere. We assume therefore that the boxmdary surface is 
r= 6 (l + €P2), where € is small, introducing the same Legendre coef¬ 
ficient as occurs in the potential of the disturbing field. 



We now regard the potential at Q as due to the spheres at A and B, 
a solid sphere of radius a and density /o — <r, on which is superposed 
another solid sphere of radius b and density a, together with a surface 
distribution of density abeP^* 

If we assume that this layer produces potentials 

= {r<b) and = (r>6), 

WO have ^^rrabePo = --5— o ) = — » 

^ \ dr dr b 

so that A = f 7r€(762, and = iTreor^Pj. 

The whole potential at Q is therefore 

F==j7r(/D~<7)^ + t7rcx(362-r2)4-f7r€CTr2P2 + y^(l+^^ ......(1). 

Let us also suppose that the whole hS/S a slow angular velocity a> 
about AB, then the equation for the pressure in the liquid is 

-2 = dF4*a>^ (xdx-\-ydy)y 
P 

so that F + Joi* (x* -f 2 /*) = const. 

represents surfaces of constant pressure, and must include the boun¬ 
dary surface. 

But a;® + y® = r® sin® ^ (1 — Pj), so that we must have 

Jtr (p - or) ~ (36* - r*) -f ^Ttear'^P^ 

+ -j(l p.) + (I - p,) = const. 

at all points on the surface r = 6 (1 + eP*). 
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A comparison of the equation of the boimdary surface with the 
expression for the disturbing field of force suggests that c is of order 
(6//)®, Hence, by substituting for r, we find that the coefficient of 
P, vanishes if ai® is also of order (6//)® and 

finable + + 2Mp - = 0 . 

Hence, when there ia no rotation, we have 


16 ibr 6 » 

* ~ 2ff/» {6 (p - v) a* + 2(76®} 


( 2 ), 


and € would be zero if there were an angular velocity given by 

But we have assumed the boundary surface to be r = 6(l + cP 2 ), 
or r® (1 — 2€pa) = 6®, which reduces to 




z® 

(l + €)® 


= 6®, 


and represents a prolate spheroid of elliptioity fc, and hence the 
required result. 


EXAMPLES 

1. Prove the formulae 

A = jf,do>. B = C = lpo>. 

where the integrations are over the whole surface of the ellipsoid. 

2. Shew that, if a>b>c, then Aa^> Bb^> (7c®; and that at points 

on the surface of such an ellipsoid the potential is greatest at the end 
of the axis c. [London XJniv. 1931] 

3. The bounding surfaces of a thin homoeoid are prolate spheroids 
and a cone is drawn having its vertex on a directrix plane. Shew that 
the portions of the homoeoid cut out by this cone exert equal attractions 
at the focus which corresponds to the directrix plane. [M. T. 1898] 

4. Prove that, if an ellipsoid is a level surface for any distribution 
of matter within it, the external level surfaces are confocal ellipsoids. 

5. Find the distribution of matter which will produce zero potential 
outside the surface a;®/a* + 2/®/&* + 2 ®/c* = 1 and potential 

J(l~aT*/o«~^®/6®-z®/c®) 

throughout the spaoe enclosed by it. 

Deduce MacLaurin’s theorem on the attractions of confocal ellip¬ 
soids. 
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6. Prove that for a prolate spheroid, [a=6 = e (1 — e*)*], 

{ ex> gSSn \ 

*“\?i(^+l)(2n + 3))- 

[London Univ. 1925] 

7. Prove that a spheroid of uniform density cannot have its 
boundary surface as one of its level surfaces. 

8. Prove that the level surfaces inside a solid homogeneous ellipsoid 
are similar and similarly situated ellipsoids. 

If from this ellipsoid a co-axial ellipsoid be removed, shew that the 
level surfaces within the cavity are hyperboloids whose asymptotic 
cones have each three mutually perpendicular generators. [C. 1907] 

9. Find the amoimt by which the gravitational potential energy 
of a uniform solid ellipsoid exceeds that of a uniform sphere of 
equal volume and mass. [M. T. 1912] 

10. A solid homogeneous ellipsoid is divided by a plane perpen¬ 
dicular to an axis. Prove that the mutual attraction of the parts for 
varying positions of the plane varies as the square of the area of the 
section. [C. 1881] 

11. Prove that the resultant attraction between the octant of the 
ellipsoid x^/a^ -h 2/^/5* + = 1 in which x, y, z are all positive and the 

rest of the body is a single force of magnitude 

*7rp*a6c (A*a2 -f 

acting along the line 

Ibnx-^ 16a__ Ibiry —165 _ ISttz — 16c 

35 “ “ Cc ' 

where p is the density assumed to be uniform. 

12. A uniform oblate spheroid is divided into two parts by its 
equatorial plane; show that, at a point of the equator, the component 
attraction of either half perpendicular to the equatorial plane is 

4ypa(l -e ) (taQjj-i g_ 

where a is the equatorial rculius and e the eccentricity of a principal 
elliptic section. [London Univ.] 

13. Shew that the resultftnt attraction of one half of a solid homo¬ 
geneous oblate spheroid, cut off by an equatorial plane, at a point on 
the rim of the base is inclined to the plane of the base at an angle whose 
tangent is 

4c (tanh”^ e — e) /(a sin“^ e — cc) tt, 

where a and c are the semi-axes and e the eccentricity of the meridian 
section. [C. 1903] 
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14. Shew that for a nearly spherical ellipsoid for which 

a = A;(l-fA), 6 = A;(l + ft), c = A;(l + >'), 

where A, fi, v are small, k^^abc and A4-jLt+v=0 approximately, the 
components of attraction at an internal point are 

•“^7rp(l-gA)a:, ~^7rp(l - Jv)2!. 

15. Prove that, if a solid uniform ellipsoid of mass M is nearly 
spherical and has semi-axes /x, {fjL^ — h)^, (p,® —A;)*, the potential at an 
external point is 

~+^^{xHh + k) + y^-2h + k) + zH-2k + h)} 
to the first order of small quantities. [M. T. 1895] 


16. An ellipsoid is cut into two equal portions, by a plane whose 
equation referred to the principal axes is lx + my -^nz — 0. Shew that 
the attraction of one half of the ellipsoid on the other reduces to a single 
force, whose x component is ^npAa^bcllp, where = 

[M. T. 1899] 


17. A solid ellipsoid of imiform density p is bounded by 
x^ja^ -f ~h z^jc^ = 1, 

and is divided by the plane y = x tan a. ProVe that the force parallel 
to this plane on either part required to prevent it from sliding over the 
other part is 

-6°) sing COS g " udu _ 

* ^ {a“sm*a + 6*cos*a)*y 0 {(a®+ m)*(6“ + m)®(c®+ «)}*' 


18. Prove that the potential of an ellipsoid of uniform density p 
at an internal point (x,y,z) exceeds the potential of a thin shell of the 
same mass M coincident with its surface and bounded by confocal 
ellipsoids by 


where a, b, c are the semi-axes of the ellipsoid and its principal axes 
are the axes of co-ordinates. [M. T. 1894] 


19. Prove that, if a sector of a homogeneous solid elliptic cylinder 
of infinite length is cut out by planes through the axis and through 
conjugate radii CP, CP' of the cross-section, then the resultant force 
on this sector is constant for all positions of the cutting planes, and 
its direction bisects the angle QCQ', where CQ, CQ' are the corre¬ 
sponding radii of the auxiliary circle of the cross-section. 

[M. T. 1904] 

20. Prove by direct integration or otherwise that the \5omponents 
of attraction at an internal point {x, y, z) of a solid homogeneous cylinder 
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of density p, whose cross-section is the ellipse -|- 3 /^/ 6 * = 1 , and 
whose length is infinite in both directions, are 


— 47ryp 


ab X 
a + ba* 


— 4Tryp 


ab y 
0 + 66* 


0 . 


If (a;, 3 /, z) is an external point, and a', b' are the semi-axes of the 
cylinder through the point which is confocal with the given cylinder, 
verify that, if X', Y\ denote respectively 


then 


— 4t'nyp 


ab X 


dX' 

dy 


dY' 

dx 


and 


ab y 

dX' dY'_ 
dx dy 


Hence, or otherwise, shew that X', Y' are the components of attraction 
at the external point. [London Univ.] 


21 . Prove that, if the surface density of an elliptic disc of semi-axes 
a, b is inversely proportional to (\-‘X^jd^ — y^lb^)^t it has the same 
attraction at points external to both, as a thin homoeoid whose focal 
ellipse is the boundary of the Klisc. 


22. Shew that any plane divides a homogeneous solid ellipsoid 

into two parts such that the attraction between them reduces to a single 
force. [C. 1891] 

23. The space between a solid homogeneous ellipsoid and a con¬ 
centric ellipsoidal envelope is filled with a homogeneous liquid; prove 
that the resultant attraction between the solid and liquid reduces to 
a couple, which vanishes only when the ellipsoids are co-axial. 

[C. 1904] 


24. Shew, by the method of 7-63, that an oblate spheroid of small 
ellipticity is a possible form of relative equilibrium of a mass of 
homogeneous gravitating liquid rotating slowly. 

26. A particle of meiss M is held at a fixed distance r(>a) from the 
centre C of a uniform circular disc, and the line joining the particle 
to G makes an angle B with the axis of the disc. Shew that, neglecting 
terms containing fifth or higher power of (a/r), the couple about G 
tending to increase B is given by 


where p and a are the surface density and radius of the disc and Q is 
the constant of gravitation. [P. 1935] 

26. If the law of attraction were ymmWy shew that a finite mass of 
fluid, which is rotating with uniform angular velocity o> about a fixed 
axis, could only be in relative equilibrium if cn* were less than yAf, M 
being the total mass of the fluid. Shew that the excess of the pressure 
at the origin ^over that at the surface is for a given mass of fluid pro¬ 
portional to (yAf — ’ [C. 1915] 
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27. A uniform gravitating liquid sphere of radius a and density p 
is made to rotate with small uniform angular velocity its surface 
being free. Find the form of the free surface, and shew that the pressure 
at the point ( Ja, 6) is 

^trOp^d^ — ^pcj^a^ — ^ffpto^d^ cos® 6, 

where G denotes the constant of gravitation. [M. T. 1928] 

28. A homogeneous gravitating solid is in the form of a prolate 

spheroid of small ellipticity €, where € = (a — c)la, and a, c are the 
semi-axes of a meridian section. The spheroid is rotating about its 
axis with angular velocity a>. Prove that the direction of a plumb line 
at any point on the spheroid will pass through the centre provided 
ai® = 6€flf/6a, where g is the approximately constant value of gravity 
at the surface. [M. T. 1927] 

29. A planet is at the centre of a imiform plane annulus of mass M 

and inner and outer radii a and 6, large compared with the planet’s 
radius. Shew that to a first approximation the effect of the ring on the 
figure of eqxiilibrium of the planet is equivalent to that of a rotational 
velocity a>, where = ZQM jah (a + 6), G being the constant of gravita¬ 
tion. [M. T. 1926] 


30. The surface of a uniform solid of mass Af is an oblate spheroid 
of semi-axes a(l —c), a, a, where c is small; find, correct to the first 
power of €, the gravitational potential of the body at an external point. 

The body rotates slowly with constant angular velocity w about its 
axis of symmetry A'OA^ O being the centre of the spheroid, and a 

particle rests on the surface at a point P such that AOP = B. Find an 
approximate expression for the force exerted by the particle on the 
body, and shew that its direction is inchned to OP at an angle 



sin 2^, 


approximately, where y denotes the constant of gravitation. 

[M. T. 1933] 


31- Two masses ikf are placed at distances c on opposite sides of the 
centre of a gravitating sphere of liquid of radius a and total mass M\ 
Shew that, neglecting powers of ale above the third, the liquid is 
deformed into a prolate spheroid, the ratio of the min or axis to the 
major being / 16Ma»\ 

V 2Af'cV‘ 

If two more particles of the same mass M are placed also at distances 
c from the centre on an axis at right angles to the fine joining the former 
pair, shew that the surface becomes an oblate spheroid, the ratio of 
the axes being the same as in the previous case. [C. 1980] 
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